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Antiferromagnetic order and spin dynamics in iron-based superconductors 
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High-transition temperature (high-T c ) superconductivity in the iron pnictides/chalcogenides 
emerges from the suppression of the static antiferromagnetic order in their parent compounds, 
similar to copper oxides superconductors. This raises a fundamental question concerning the role 
of magnetism in the superconductivity of these materials. Neutron scattering, a powerful probe to 
study the magnetic order and spin dynamics, plays an essential role in determining the relationship 
between magnetism and superconductivity in high-T c superconductors. The rapid development of 
modern neutron time-of-flight spectrometers allows a direct determination of the spin dynamical 
properties of iron-based superconductors throughout the entire Brillouin zone. In this review, we 
present an overview of the neutron scattering results on iron-based superconductors, focusing on 
the evolution of spin excitation spectra as a function of electron/hole-doping and isoelectronic 
substitution. We compare spin dynamical properties of iron-based superconductors with those 
of copper oxide and heavy fermion superconductors, and discuss the common features of spin 
excitations in these three families of unconventional superconductors and their relationship with 
superconductivity. 
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I. INTRODUCTION 


The interplay between magnetism and superconduc¬ 
tivity has a long history. For example, it is well known 
that superconductivity in conventional Bardeen-Cooper- 
Schrieffer (BCS) superconductors (Bardeen et al. 1957) 
such as the element lanthanum can be suppressed by as 
little as 1% gadolinium substitution as the magnetic im¬ 
purity (Mathias et al. 19581. Within the BCS frame¬ 


work, magnetic impurities can act as pairing break¬ 


ing agents rapidly suppressing superconductivity (Bal- 
atsky et al. 2006). In the case of unconventional 


su¬ 


perconductors such as heavy fermions (Gegenwart et al. 


2008 Steglich et al. 1 

979 Stewart 2001), copper oxides 

(Bednorz and Muller 

1986 Lee et al., 2006 

, and iron 

pnictides/chalcogenic 

3S (Canfield and Bud’ko 

2010 Hsu 

et al., 2008 Johnston 

2010 Kamihara et al., 2006 2008 

Paglione and Greene 

2010j Stewart 2011; Wen et al. 


2011), the observation that superconductivity always ap¬ 


pears near the static antiferromagnetic (AF) order (Ue- 


mura 


2009) suggests that magnetism may be a common 


thread for understanding the microscopic origin of un¬ 
conventional superconductors and high-transition tem¬ 


perature (high-T c ) superconductivity (Scalapino 2012). 


Based on this premise, much work has focused on study¬ 
ing the interplay between magnetism and superconduc¬ 
tivity, particularly the high-T c copper oxide supercon- 


ductors (cuprates) since its discovery in 1986 1 

Armitage 

et al, 2010; Fujita et al., 2012; Kastner et al., 1998; Kivel- 

son et al. | 2003; Tranquada et al. 

2014 

). Although under- 


standing the magnetism and its relationship with super¬ 
conductivity in cuprates is still an area of active research 
(Fujita et al. 2012[ Tranquada. et al. 20141, the discovery 
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of AF order in the parent compounds of iron-based super- 2008). 


conductors in 2008 

de la Cruz et al. 

2008 Huang et al. 

2008a Zhao et al. 

2008b 

) provided a new opportunity 


to study the interplay between magnetism and supercon¬ 
ductivity. There are already several review articles sum¬ 
marizing the general progress in the field of iron-based 


superconductors (Canfield and Bud’ko 2010 Chubukov 

2012 

Dagotto, 2013 

Hirschfeld et al, 2011; 

lohnston 

2010 

Mazin 2010 Paglione and Greene, 2010 

Stewart 

2011 

Wang and Lee 

2011 Wen et al. 2011). The inter- 


play between magnetism and superconductivity studied 
b y neutron scattering has b een briefly discus sed as well 
(Dai et al. 2012 Lumsden and Christianson 2010). 


The purpose of this review article is to present a 
comprehensive account of the AF order, spin excita¬ 
tions, and their relationship with superconductivity in 
iron pnictides and chalcogenides. Since magnetism is 
generally believed to play an important role in the 
electron pairing mechanism of high-T c superconductiv¬ 
ity (Scalapino 2012), it is important to summarize the 


progress in the field over the past several years and com¬ 
pare the outcome with previous work on high-T c copper 
oxides and heavy fermion superconductors. Neutrons- 
with their wavelengths comparable to the atomic spac¬ 
ing and their spins directly probing the unpaired elec¬ 
trons in solids-have played a unique role in determin¬ 
ing the magnetic properties of high-T c superconductors. 
Soon after the discovery of iron pnictide superconduc¬ 
tor LaFeAsOi-zFj, with T c = 26 K (Kamihara et al. 


2008), neutron and X-ray scattering experiments have 


discovered that its parent compound LaFeAsO exhibits a 
tetragonal-to-orthorhombic structural distortion followed 


by a collinear AF order (de la Cruz et al. 2008 No¬ 


mura et al. 2008). Since the presence of a collinear AF 


structure in LaFeAsO was predicted earlier in the band 
structure calculations as due to a spin-density-wave or¬ 
der arising from nesting of the hole and electron Fermi 
surfaces (Dong et al. 2008), its confirmation by neutron 


scattering and the semi-metallic nature of these materials 
(Kamihara et al. 2008) provided strong evidence for the 


itinerant origin of the magnetism in the iron based super¬ 
conductors (Hirschfeld et al. 2011 Mazin 20101. This 


is fundamentally different from the parent compounds of 
copper oxide superconductors, which are Mott insulators 
because the Coulomb repulsive energy cost U of having 
two electrons (or holes) on the same site is much larger 


than the electron hopping energy t (Lee et al. 

For a Mott insulator, the static AF order arises from a 
saving in energy of 4t 2 /U via virtual hopping and is due 
to electron correlation effects. Since the Mott insulating 
state of copper oxides is believed to play an essential role 
in the pseudogap physics and mechanism of high-T c su¬ 
perconductivity (Lee et al. 2006), it is interesting to ask 


whether iron-based superconductors are also close to a 
Mott insulator and determine the effect of electron cor¬ 
relations and local moments to their electronic properties 


and spin dynamics ( 

Fang et al. 

2008 

Haule et al. 

2008 

Qazilbash et al., 2009J Si and Abrahams 2008 Xu et al. 


From the experimental point of view, a systematic 
determination of the magnetic structures and the dop¬ 
ing evolution of spin excitations in different classes of 
iron-based superconductors and their associated materi¬ 
als will form the basis to establish whether magnetism 
is responsible for high-T c superconductivity. For cop¬ 
per oxides, superconductivity can be induced by charge 
carrier doping (electrons or holes) into the Cu0 2 , re¬ 
sulting complicated phase diagrams with many incipi¬ 
ent states competing with superconductivity (Armitage 


et al.\ 2010; Fujita et al., 2012f Kastner et a/.[|1998||Kivel- 

son et al. j2003 

Tranquada et al. 

2014). The undoped 

copper oxides such as La 2 CuC>4 

Vaknin et al. 

1987) 


and YBa 2 Cu 3 0 6+a; (Tranquada et al. 1988) are simple 


antiferromagnets with the neighboring spins oppositely 
aligned. When holes are doped into the parent com¬ 
pounds, the static AF order is gradually suppressed, but 
spin excitations (or short-range spin fluctuations) sur¬ 
vive across the entire superconducting dome and couple 
with superconductivity via a collective magnetic excita¬ 
tion termed neutron spin resonance (Eschrig 2006 Fujita 
et al. 2012 Tranquada et al. 2014). However, inelas¬ 


tic neutron scattering experiments designed to study the 
doping evolution of spin excitations were only carried out 
on the La 2 _ x Sra;Cu04 family of cuprates across the en¬ 
tire phase diagram from the undoped parent compounds 
to heavily overdoped non-superconducting samples (Lip- 


scombe et al. 

2007 

Wakimoto et al. 

2007) 


There are 


phase diagram on other cuprates due to material limita¬ 
tions (for example, YBa 2 Cu 3 06 +a; cannot be hole over¬ 
doped to completely suppress superconductivity) or the 
difficulty in growing large single crystals suitable for in¬ 
elastic neutron scattering experiments. 

In the case of iron-based superconductors, there are 
two major classes of materials, the iron pnictides and iron 
chalcogenides (Johnston |2010[ Stewart| 2011). Com¬ 


pared with the hole-doped La 2 _ x Sr ;c Cu04 copper oxide 
superconductors, where superconductivity can only be 
induced via substituting the trivalent La by the diva¬ 
lent element Ba or Sr, superconductivity in iron pnictide 
such as BaFe 2 As 2 (Rotter et al. [2008b) can be induced 


by ionic substitution at any element site. These include 
Ba by K/Na to form hole-doped Bai_ a: K a .Fe 2 As 2 ( |Pra-| 
manik et al. 2011 Rotter et al. 2008a), Fe by Co or Ni 


20061. to have electron-doped BaFe 2 _ a; r a ;As 2 (T =Co, Ni) (TLJ 


et al. 

2009a 

Sefat et al. 

2008 

), and As by P in the isova- 

lent ( 

or isoelectronic) doped compounds BaFe 2 As 2 _ a ,P a; 

(Jiang et al. 

2009 

). While K or Na doping to BaFe 2 As 2 


induces the same numbers of holes to the FeAs layer, 
the effect of Ni-doping is expected to introduce twice 
the number of electrons into the FeAs layer as that of 
Co-doping from naive electron counting point of view 
(Johnston 2010| Stewart 2011). Since large sized sin¬ 
gle crystals can be grown by self-flux method in many of 
these cases (Canfield and Bud’ko 2010), doped BaFe 2 As 2 


materials provide a unique opportunity to study the evo- 
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lution of the static AF order and spin excitations as a 
function of hole, electron, and isovalent doping through¬ 
out the entire AF order to superconductivity phase dia¬ 
gram, and determine their connection with superconduc¬ 
tivity. These experiments, together with neutron scat- 


tering studies of related materials (Inosov et al. 2013 

Ishikado et al.\ 2009 

; Johnston et al 2011 [ Kim et al, 

2011; Lamsal et al. 

2013 Marty et al. 2011 Shamoto 

et al. 

2010; Simonson et al., 2012| 

Singh et al. 2009 

Taylo 

r et al.\ 2013 

Wakimoto et al. 

2010), will estab- 


excitations in different families of iron-based supercon¬ 
ductors. The outcome, together with the results from 
high-T c copper oxide and heavy fermion superconduc¬ 
tors, and can form a basis to determine if magnetism is 
indeed responsible for superconductivity in these materi¬ 


als (Scalapino 2012). 


Compared with other techniques suitable to study 
magnetism in solids including Muon spin rotation (/xSR) 


(Uemura 


loul et al. 


2009), nuclear magnetic resonance (NMR) (Al- 


2009), and resonant inelastic X-ray scattering 


(RIXS) (Ament et al. 2011), neutron scattering has sev¬ 
eral unique advantages: (1) The neutron itself is a charge 
0 fermion with a spin S = 1/2, resulting a magnetic 
dipole moment which can interact with unpaired elec¬ 
trons and magnetism in solids; (2) The charge neutral¬ 
ity of the neutron renders it a weakly interacting probe 
with well-known scattering cross sections; (3) The avail¬ 
able wavelength and energies of neutrons as a scattering 
probe are ideally suited to study static magnetic order 
and spin excitations in solids. The general scattering 
principle involved is simply to measure the number of 
neutrons scattered into a given solid angle at a known en¬ 
ergy (E = hui, where h is the reduced Planck’s constant 
and u! is the angular frequency) and momentum transfer 
(Q). The laws of conservation of momentum and energy 
are satisfied via Q = k, — kf and E = hoj = Ei — Ef , 
where ki , Ej = h 2 k 2 /2m, kf , and Ef = h 2 k 2 /2m are the 
incident neutron wave vector, energy, outgoing neutron 
wave vector, and energy, respectively, and m is the mass 
of a neutron. The coherent magnetic scattering cross 
section from a system with a single species of magnetic 
atoms is then (Xu et q£j|2013a), 


GT<7 

dfldE 


N kf 2 -2 w 
hh p e 

IL /vj 


£(*a,0-QaQ/OS a/, (Q,w). 


a ,/3 


Here N is the number of unit cells, p = (-^) 2 g 2 f(Q) 2 
(where = 0.2695x 10” 12 cm, g k, 2 is the electron spin 
^-factor, and /(Q) is the magnetic form factor), e~ 2W is 
the Debye-Waller factor, a, /3 are the Cartesian coordi¬ 
nates x, y, and z, and Q a , Qp are the projections of the 
unit wave vector Q onto the Cartesian axes. S a P (Q, w) is 
the dynamic spin correlation function, and is associated 
with the imaginary part of the dynamic susceptibility 
x"^(Q,w) via the fluctuation-dissipation theorem: 

X'/PQ,^ = g 2 y 2 J(1 - e-WBT)^ Q, w ). 


For a paramagnet with isotropic spin excitations, 
S xx (Q,lu) = S yy (Q,co ) = S zz (Q,u). Since neutron 
scattering is only sensitive to spin (fluctuations) direc¬ 
tion perpendicular to the wave-vector transfer Q, the 
5(Q,w) of an isotropic paramagnet measured by unpo¬ 
larized neutron scattering experiments (see Section III. F 
for neutron polarization analysis) is related to S' S2 (Q, u>) 
via S(Q,u) = 2S’ 2Z (Q,w). By measuring S( Q,w) in 
absolute units via phonon or vanadium normalization 


(Xu et al., 2013a), one can estimate the energy depen¬ 


dence of the Q-averaged or the local dynamic susceptibil¬ 
ity X'"(w) = fp,? x" (Q’ u )dQ/J BZ d Q within a Brillouin 
Zone (BZ) ( Lester et al.\ 2010). The overall strength of 
the magnetic excitations, corresponding to the loc al fluc- 
tuating moment (m 2 ), can then be computed via (Lester 
et al. 1120101) 


<™ 2 > = “ / 


X"M du 


-oo 1 - exp {-Tiio/k B T) 


One of the central purposes of inelastic neutron scattering 
experiments is to determine the energy and wave vector 
dependence of x w (Q, w) in absolute units for various iron 
pnictides, and compare the outcome with those in copper 
oxide and heavy fermion superconductors. 

In this article, we present a comprehensive review of 
recent neutron scattering results on iron-based supercon¬ 
ductors, mainly focusing on the evolution of the static 
AF order and spin dynamics of iron pnictides and chalco- 
genides. In Section II, we summarize the static AF or¬ 
der for various iron pnictides/chalcogenides and its dop¬ 
ing evolution. This includes the effects of electron and 
hole-doping on the static AF order and tetragonal-to- 
orthorhombic structural transitions (Section II. A.); how 
impurity (Section II. B.) and isoelectronic substitution 
(Section II. C.) affect the magnetic and structural phase 
transitions. Section III summarizes spin excitations and 
their relationship with superconductivity, including spin 
waves in the parent compounds of iron-based supercon¬ 
ductors (Section III. A.); as well as neutron spin reso¬ 
nance and its relationship with superconductivity (Sec¬ 
tion III. B.); the electron and hole-doping evolution of 
the spin excitations in the BaFe 2 As 2 family of iron pnic¬ 
tides (Section III. C.); evolution of spin excitations in 
iron chalcogenides and alkali iron selcnides (Section III. 
D.); impurity effects on spin excitations of iron pnictide 
and chalcogenide superconductors (Section III. E.); neu¬ 
tron polarization analysis of spin excitation anisotropy in 
iron pnictides (Section III. F.); electronic nematic phase 
and neutron scattering experiments under uniaxial strain 
(Section III. G.); comparison with /rSR, NMR, and RIXS 
measurements (Section III. H.); and comparison of spin 
excitations in iron-based superconductors with those in 
copper oxide and heavy fermion superconductors (Sec¬ 
tion III. I.). Section IV provides a brief account of cur¬ 
rent theoretical understanding of spin excitations in iron- 
based superconductors. Finally, we summarize the re¬ 
sults and discuss possible future directions for the field. 
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FIG. 1 (Color online) Crystal structure and magnetic order 
in different families of pnictides. (a) The crystal and magnetic 
structures of the BaFe 2 As 2 in the AF ordered phase. The yel¬ 
low, green, and blue balls indicate Ba, As, and Fe positions, 
respectively. The red arrows mark the ordered moment di¬ 
rections of Fe in the AF ordered state (the C-type). The 
ay, &t, and c show the Cartesian coordinate system suitable 
for the paramagnetic tetragonal phase of BaFe 2 As 2 (Huang| 


et al., 2008a I. (b) The AF structure of BaMn 2 As 2 , where 


the ordered mo ments on Mn are along the c-axis direction 
(the G-type) (Singh et al. |2009 >. (c) The crystal structure 
of CaCo 2 As 2 , where the ordered moments on Co form the A- 
type AF structure ( |Quirinale et al] 1 2013 [ ) . (d) The collinear 
AF order in NaFeAs doubles the crystalline unit cell along 
the c-axis (Li et al. 2009c). (e) The collinear C-type AF 


structure in the Fe plane, where the green dashed box marks 
the tetragonal crystalline unit cell in the paramagnetic state 
and the magenta dashed box indicates the orthorhombic mag¬ 
netic unit cell. The u 0 ,af and 6 0 ,af mark directions of the 
orthorhombic lattice, (f) The in-plane moment projections 
for the G-type antiferromagnets. 


II. STATIC ANTIFERROMAGNETIC ORDER AND ITS 
DOPING EVOLUTION 

A. Lattice and magnetic order in the parent compounds of 
iron-based superconductors 
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FIG. 2 (Color online) Crystal and magnetic structures of 
iron chalcogenides and alkali iron selenides. (a) The crystal 
structure of the FeTei-aSe^, iron chalcogenide, where Fe and 
Se/Te positions are marked as blue and green, respectively, 
(b) The tetragonal phase crystal structure of A r Fe 2 - :E Se 2 . 
The positions of A, Fe, and Se are marked as light blue, blue, 
and green, respectively, (c) The in-plane bi-collinear magnetic 
structure of FeTe, where the arrows indicate the moment di¬ 
rections (|Bao et al . | |2009| |Fruchart et al.\ |1975| |Li et al\ 


2009d|. (d) The in-plane crystal and magnetic structures of 
A 2 Fe 4 Se 5 in the vacancy ordered insulating phase. Only iron 
positions are plotted and the grey dashed lines mark the struc¬ 
tural and magnetic unit cells. The light blue (magenta), green 
(blue), and yellow (dashed yellow) lines represent the nearest 
(Ji, Ji), next nearest (J 2 , J 2 ), the next next nearest (J 3 , 
J 3 ) neighbor exchange interactions, respectively (|Bao et al. 


2011 |Wang et al. 201 lb|c[ Ye et al. 20111. (e) The crysta 


and magnetic structure s of A 2 Fe 3 Se 4 in t he vacancy ordered 
semiconducting phase (Zhao et al. 2012). The nearest and 


next nearest neighbor exchange couplings are clearly marked. 
(f,g) Other possible magnetic structures of Tl 2 Fe 4 Ses in the 
vacancy ordered phase (|May et~aT. 20121. 


From a crystal structure point of view, the parent com¬ 
pounds of iron-based superconductors can be classified 
into five different families: RFeAsO (R = La, Ce, Pr, 
Nd, Sm,..., the 1111 system), AFe 2 As 2 (A = Ba, Sr, Ca, 
K, the 122 system), AFeAs (A = Li, Na, the 111 system), 
Fei+yTei-^Se^ (the 11 system), and A a ,Fe 2 _ y Se 2 alkali 
iron selenides (A = K, Rb, Cs, Tl, ..., including the insu¬ 


phase A 2 Fe 3 Se 4 ) (Aswathy et al. 

2010} 

Dagotto, 

2013 

Johnston, 2010; Paglione and Greene, 2010( Sadovskii 

2008 Stewart 2011| Wang et al. 

2014 iZhao et al ., 

2012) 


ers in the unit cell and other systems have single FeAs(Se) 
layer. A recent development in the field is the synthe¬ 
sis of iron selenide superconductors via intercalation of 
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molecular complexes between layers of FeSe (Burrard- 
Lucas et al.\ 2013[ Krzton-Maziopa et al. 2012} |Yiiig 
et al. 2012). The crystal structures at room tempera 


axis [Fig. 2(d)] (Bao et al. 2011| |Wang et al. 2011b 


ture are all tetragonal except for the insulating 245 phase 
and some of them will become orthorhombic at low tem¬ 
perature below T s . Neutron diffraction measurements 
have established that the long range AF order in the 
iron pnictides including the 1111, 122, 111 families is 
collinear with moment aligned along the a 0 axis direc¬ 
tion of the orthorho mbic s tructure [Figs. 1(a), 1(d), and 
1(e)] ( Lynn and Dai| 2009), except for the stoichiometric 
LiFeAs system which is superconducting without a mag¬ 


netically ordered parent compound (Pitcher et al. 2008 


Tapp et al. 2008} Wang et al. 2008). Although the in¬ 


plane collinear AF structures for different classes of iron 
pnictides are identical [Fig. 1(e)], the ordering along the c 
axis is material dependent. In the 122 system, which has 
two magnetic irons per formula unit, the ordering is AF 
within a unit cell along the c axis [Fig. 1(a)]. For the 111 
system with one iron per unit cell, the magnetic struc¬ 
ture doubles the chemical unit cell along the c axis [Fig. 
1(d)]. While the collinear AF structure in iron pnictide 
is the so-called C-type antiferromagnet stemming from 
the original work of Wollan and Koehler on perovskite 
manganese oxides (Wollan and Koehler 1955), the re¬ 


lated pnictide materials such as BaMn 2 As 2 (Singh et al. 
2009) and CaCo 2 As 2 (Quirinale et aZ. |2013 ) have the G 


[Figs. 1(b), 1(f)] and A-type [Fig. 1(c)] AF structures, 
respectively. Recently, another AF parent compound was 
found in electron-overdoped LaFeAsOi-^Ha, (x ~ 0.5) 
system in addition to the usual collinear AF structure at 
x = 0 (Hiraishi et al. 20141. 


For the iron chalcogenides (the 11 family) and alkali 
iron selenides, their crystal structures are shown in Figs. 
2(a) and 2(b), respectively. Instead of a collinear AF 
structure, the parent compound of the 11 family has a 
bi-collinear AF spin structure as shown in Fig. 2( c) ( |Bao 
et al . | 2009} Fruchart et al . | 1975 Li et al. 2009d). Com 
pared to the collinear spin structure of the 122 family 
in Fig. 1(e), the iron spins are rotated 45° within the 
a 0 b 0 -plane in the 11 system. The magnetic structure in 
the 11 family is sensitive to the excess iron population 
in the interstitial iron site (Rodriguez et al. 2011 20101. 


While the bi-collinear magnetic order is commensurate 
for Fei_|_ x Te with x < 9%, it exhibits incommensurate he¬ 
lical magnetic order that competes with the bi-collinear 
commensurate ordering close to TV for x > 12% (Ro- 


Ye et al. 


2011) and in the plane [Figs. 2(f) and 2(g)] 


(May et al. 2012), respectively. Compared with the par¬ 


ent compounds of the 122 system, the ordered moments 
of the 11 and 245 systems are much larger. In Table 
I, we summarize the lattice parameters, structure tran¬ 
sition temperature T s , the AF phase transition temper¬ 
ature TV, and the static ordered moments for the par¬ 
ent compounds of different iron-based superconductors. 
In the 1111, 111, and 245 systems, the structural tran¬ 
sition occurs at a temperature higher than that of the 
AF phase transition (Johnston 2010[ Stewart. 2011). In 
the 122 system, the structural and magnetic transitions 
almost occur simultaneously in the undoped parent com¬ 
pounds (Kim et al. 2011), but are well separated upon 


electron-doping (Canfield and Bud’ko 2010). 


B. The effect of electron, hole-doping, impurity, and 
isoelectronic substitution on the static AF order and 
tetragonal-to-orthorhombic structural transitions 

Before discussing the impact of electron/hole doping 
on the long range AF order, we shall define momen¬ 
tum transfer in reciprocal space and compare the sizes 
of the Brillouin zones for the parent compounds of dif¬ 
ferent families of high-T c superconductors. Figures 3(a), 
3(b), and 3(c) show spin arrangements within one layer 


of La 2 CuC >4 (Vaknin et al. 1987 

1, BaFe 2 As 2 (Huang 

et al., 

2008a), and FeTe (Bao et al.\ 

2009 

Fruchart et al. 


1975 


Li et al. 2009d), respectively. The chemical unit 


cells are marked as green dashed lines and the mag¬ 
netic unit cells are magenta shaded. The positions of 
Cu 2+ /Fe 2+ and 0 2_ /As 3_ /Te 2_ are also marked. The 
momentum transfer Q at ( q x , q y , q z ) in A^ 1 can be de¬ 
fined as ( H,K,L ) = (q x a/2n, qyb^n, q z c/2n) in recip¬ 
rocal lattice units (r.l.u.), where a (or a 0 ), b (or b a ), 
and c are lattice parameters of the orthorhombic unit 
cell. In this notation, the AF order in the parent com¬ 
pound of copper oxide superconductors occurs at Qaf = 
(iT, K) = (0.5 ± m, 0.5 ± n) where m,n = 0,1,2,--- 
and the first magnetic Brillouin zone is shown as the 
magenta shaded box in Fig. 3(d). Another equivalent 
Brillouin zone near T is marked by the magenta dashed 
line, while the first Brillouin zone of the chemical unit 
cell is the green dashed box. If the ordered moment is 
entirely on the iron site in BaFe 2 As 2 , the chemical unit 


driguez et al. 

2013 

2011 201( 

)). The alkali iron selenides cell is twice the size of the magnetic unit cell along the 

(the 245 family) ([Fang et al.\ 

>011 

Guo et al. 

2010) have b 0 axis direction due to out of plane positions of the As 


two different iron vacancy structures including the insu¬ 
lating a/ 5 x V5 iron vacancy ordered phase [Figs. 2 (d), 


2(f), and 2(g )] ( Bao et cd\ 2011} Wang et~cd\ 2011b Ye 


et al, 2011) and the semiconducting rhombus iron va- 


cancy ordered 234 phase [Fig. 2(e)] (Wang et al. 2014 


Zhao et al, 2012). While the 234 phase has a AF struc¬ 


ture similar to the parent compounds of the 122 family 
[Fig. 2(e)] (Zhao et al, 2012), the insulating 245 phase 
have the block AF structure with moments along the c 


atoms [Fig. 3(b)]. In a completely detwinned sample, the 
first magnetic Brillouin zone is the magenta shaded area 
around Qaf = (H,K,L) = (1 ± 2m, 0 ± 2 n,L) where 
L = ±1,3, 5, ••• r.l.u., larger in size than the chemical 
Brillouin zone in dashed green. Because the AF order 
in iron pnictides is always preceded by a tetragonal-to- 
orthorhombic lattice distortion, the twinning effect in the 
orthorhombic state means that AF Bragg peaks from the 
twinned domains appear at positions rotated by 90° [Blue 
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TABLE I Summary of the structure transition temperature T s , magnetic transition temperature Tn , and ordered magnetic 
moment per iron for the AF ordered parent compounds of the iron-based superconductors. The lattice parameters in the 
paramagnetic tetragonal state are also listed. 


Materials 

aT = bT (A) 

C (A) 

Ts (K) 

T n (K) 

Moment/Fe (/is) 

LaFeAsO ° 

4.0301 

8.7368 

155 

137 

0.36-0.6 

CeFeAsO b 

3.9959 

8.6522 

158 

140 

0.8 

PrFeAsO c 

3.997 

8.6057 

153 

127 

0.48 

NdFeAsO d 

3.9611 

8.5724 

150 

141 

0.25 

LaFeAsOo.sHo.g e 

3.975 

8.67 

95 

92 

1.21 

CaFe 2 As 2 s 

3.912 

11.667 

173 

173 

0.80 

SrFe 2 As 2 9 

3.920 

12.40 

220 

220 

0.94 

BaFe 2 As 2 h 

3.957 

12.968 

~140 

~140 

0.87 

Nao.gssFeAs 1 

3.9448 

6.9968 

49 

39 

0.09 

Fei.oesTe 9 

3.8123 

6.2517 

67 

67 

2.25 

K 2 Fe 4 Se 5 k 

8.7306 

14.113 

578 

559 

3.31 

Rb2Fe4Se5 1 

8.788 

14.597 

515 

502 

3.3 

Cs 2 Fe 4 Se 5 m 

8.865 

15.289 

500 

471 

3.4 

TlFei. 6 Se 2 n 

~8.71 

14.02 

463 

100 

~3 

a (|de la Cruz et al. 

12008| |Huang et al.\ 2008b |McGuire et al.\ 


i Kimber et al. ZUU8 Zhao et al. ZUU8c 


(Chen et al. 


(Hiraishi et aTTftZUITl 


Goldman et al. 2008 2009 Kreyssig et al., 2008). 

Jesche et al. Zl)U8 Kancko et at. Z008, Zha o et~ak| |2008d[ ). 


Sow 

i 


(Huang et al. 


(Li et aT. 


et al. 


ig - 


2008a 


2009c). 


Kim et al. 


2dnf" 


9 

Bao et al. 

2009 

Li et al. 

2009d 


*■ 

Bao et al. 

2013 

. 



1 ■ 

Wang et al. 

2011b 

Ye et al. 

2011 


Ye et at., 2U 
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n [ 

May et aC. ZUJ 

2). 




dots in Fig. 3(g)]. Therefore, to properly account for the 
twin domain effect, one needs to carry out wave vector 
integration within the region marked by solid black lines 
in Fig. 3(g) to obtain the local dynamic susceptibility 
x"(w). Figure 3(c), 3(f), and 3(h) summarizes the bi- 
collinear spin structure of FeTe, its associated magnetic 
Bragg peaks in reciprocal in detwinned, and twinned 
samples, respectively. Depending on the size of the unit 
cell, the same AF Bragg peak for collinear AF order in 
BaFe 2 As 2 and and bi-collinear AF order in FeTe can ap¬ 
pear with different Miller indices. For example, if we 
choose unit cells of BaFe 2 As 2 with one [half the magenta 
shaded area in Fig. 3(b) and ignoring As], two [similar to 
magenta shaded area in Fig. 3(a)], and four irons [dou¬ 
bling the magenta shaded area in Fig. 3(b) along the 
c-axis], the same AF Bragg peak would occur at in-plane 
wave vectors (0.5,0), (0.5, 0.5), and (1,0), respectively. 
For the bi-collinear AF order FeTe, one iron and two 
irons per unit cell would have the same AF Bragg peak 
at the in-plane wave vectors (0.25,0.25) and (0.5,0), re¬ 
spectively. 

Figure 4(a), 4(b), and 4(c) summarizes the effective 


nearest neighbor and next nearest neighbor magnetic ex¬ 
change couplings for La 2 CuC> 4 , BaFe 2 As 2 , and FeTe, re¬ 
spectively. Figure 4(d), 4(e), and 4(f) shows the cor¬ 
responding reciprocal space with Qaf positions marked 
as red and blue dots for the two different twin domains. 
While neutron scattering typically studies the magenta 
region of the reciprocal space within the first Brillouin 
zone near Qaf, RIXS can only probe spin excitations 
within the purple circles near the origin T due to kine¬ 
matic constraints from energies of the incident and out¬ 
going photons used to enhance the magnetic scattering 
at Cu £ 3 ~edg e and Fe T 3 -edge (Tacon et al. 2011 Zhou 
et al.. 2013aI. 


Although the field of iron-based superconductors 
started with the discovery of the 1111 family of mate¬ 
rials (Kamihara et al. 2008), majority of recent neu¬ 


tron scattering work has been focused on the 122 fam¬ 
ily due to the availability of high quality single crystals 
(Canfield and Bud’ko 2010). In the undoped state, the 


prototypical 122 such as BaFe 2 As 2 exhibits tetragonal- 
to-orthorhombic lattice distortion and AF order below 
T s ~ Tn ss 138 K (Huang et al. 2008a). Figure 5 sum- 






































































































7 


La 2 Cu0 4 

o ® o @ 

e o • 

b #-<>-♦ 
o ? ~ * 


(b) BaFe As 


(d) 


Detwinned 


o 


-1 0 1 
H (2t i/a) 


(e) 


4 - 4 


B 


(g) 


Cu 2+ / Fe 2+ 

0 2 7 As 3 7Te 2 ' 

AF unit cell/BZ 
Lattice unit cell / BZ 

Zone for x"(to) 
calculation 

Twinned 


(a) 


La CuO„ 

2 4 


(d) 


- 2-1 0 1 



® 


n ©■■■■# ® 

-4► a 


£ 0 


- 2-10 1 2 
H (27t/a) 


(b) BaFe 2 As 2 


(C) 


FeTe 

* ® 


4 *^ 

I ® 

o ® 





(C) 


FeTe 


(f) 


- 2-10 1 2 
H (2jt/aJ 


H UK/a M ) H (2n/aj 


94 


9 

^0.5- 

FIG. 3 (Color online) Comparison of the AF structures of 
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marked as dashed green and pink area, respectively (Fujita 
et aZ.| 2012]). (d) The reciprocal space where the solid red dots 


represent the AF ordering wave vectors. The magenta area 
and dashed green indicate the size of the in-plane magnetic 
and chemical Brillouin zone, respectively, (b) The in-plane 
AF structure of BaFe 2 As 2 , where the open and filled blue cir¬ 
cles indicate As positions below and above the iron plane, re¬ 


spectively (Johnston 2010). The magnetic and chemical unit 


cells are marked as magenta area and dashed green lines, re¬ 
spectively, and (e) the corresponding reciprocal space, where 
red dots indicate Qaf- (g) The effect of twin domains for 
AF order and Brillouin zones. The solid black lines mark the 
integration area in reciprocal space to obtain x" (w). ( c ) The 
in-plane AF structure of FeTe, and (f) the corresponding re¬ 
ciprocal space in a detwinned sample, (h) The effect of twin 
domain in reciprocal space. The magenta dashed boxes in 
(d), (e), and (f) indicate the AF Brillouin zones near F point 
probed by RIXS. 


marizes evolution of the structural and magnetic phase 
transitions for the electron and hole doped BaFe 2 As 2 - 
From transport, neutron diffraction, and X-ray diffrac- 


-1 -0.5 0 0.5 1 
H (27i/aJ 


FIG. 4 (Color online) Real and reciprocal space of La 2 CuC> 4 , 
BaFe 2 As 2 , and FeTe with magnetic exchange couplings and 
regions of reciprocal space probed by neutron scattering and 
RIXS. (a) The in-plane AF structure of La 2 Cu 04 , where near¬ 
est and next nearest neighbor magnetic exchange couplings 
are marked as Ji (solid blue line) and J 2 (solid green line), 
respectively, (d) The reciprocal space where the solid red 
dots represent the AF ordering wave vectors. The purple cir¬ 
cular area at T indicates the region of reciprocal space that 
can be probed by RIXS using Cu L-edge, while the AF Bril¬ 
louin zo ne probed by neutr on scattering is marked by magenta 
square (Tacon et al. 2011). (b) The in-plane AF structure 


tion measurements (Canfield and Bud’ko 

2010 

Chris- 

tianson et al.\ 2009 

Chu et a/.||2009; Lester et al. 

05 

O 

O 

CM 

Nandi et al. 2010; 

Pratt et al. 

2009 

1, the phase dia- 


gram or eiecirun-uopeu joa^rei-^oa;J 2 AS 2 as shown in 
Fig. 5(b) has been established. Upon electron-doping 
via Co substitution for Fe to suppress the static AF or¬ 
der and induce superconductivity, the structural and AF 
phase transitions are gradually separated with the struc¬ 
tural transition occurring at higher temperatures than 
the magnetic one. The collinear static AF order co¬ 
exists and competes with superconductivity in the un- 


of BaFe 2 As 2 with the nearest neighbors (Ji a ,Jib) and sec¬ 
ond nearest neighbor J 2 magnetic exchange couplings, (e) 
The corresponding reciprocal space, where the purple circle 
indicate reciprocal space area covered by RIXS using Fe L 3 - 
edge (Zho u et al.\ 2013a). (c) The in-plane AF structure of 
FeTe with the nearest neighbors (Ji a ,Jib) and second near¬ 
est neighbors (J 2 a, J 2 b) magnetic exchange couplings, (f) The 
corresponding reciprocal space with the purple circle showing 
the reciprocal space covered by RIXS. 


derdoped regime marked as green shaded area in Fig. 

|2009 ). 


5(b) (Christianson et al. 2009 Pratt et al. |2009 1. For 
electron-doping levels near optimal superconductivity, 
the orthorhombic lattice distortion S = (a — b)/(a + b ) 
initially increases with decreasing temperature below 
Tat, but then decreases dramatically below T c . For 
BaFei. 874 Coo.i 26 As 2 , the orthorhombic structure evolves 
smoothly back to a tetragonal structure below T c and 
the system is believed to enter into a “reentrant” tetrag- 
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(a) 3«i r 


(c) Ba, Na Fe.As, * \ BaFe 2 _ x Ni,As 2 

“ W: 



FIG. 5 (Color online) The structural and magnetic phase 
diagrams of electron and hole-doped BaFe 2 As 2 - (a) The 
coupled structural and AF phase transitions in hole-doped 
Bai_^Ka;Fe 2 As 2 as determined from neutron powder diffrac¬ 


tion experiments (Avci et al, 2012). The structural and AF 


phase transitions are denoted as T s and TV, respectively, (b) 
The phase diagram of Ba(Fei_ a ,Co a: ) 2 As 2 determined from X- 


ray and neutron diffraction experiments (Nandi et al. 2010). 


(c) The structural and magnetic phase diagram of hole-doped 
Bai_ ir Na :! ;Fe 2 As 2 from neutron diffraction experiments (Avci 
\etaJ] [20l3l [2014| ) . The green shaded region denotes the pres 
ence of a tetragonal AF phase, (d) Similar phase diagram for 
electro n-doped BaFe 2 - ;c Ni; C As 2 (Lu et al. |2013 


Luo et al. 


2012b |. Here the incommensurate (IC) AF order is a spin- 


glass phase coexisting and competing with the superconduct¬ 
ing phase ( |Lu et al ] 2014b). 


onal phase as shown in Fig. 5(b) (Nandi et al. 2010). 


Subsequent neutron diffraction experiments revealed that 
the static AF order in the underdoped regime changes 
from commensurate to transversely incommensurate for 


Ba(Fei_ a; Co a ;) 2 As 2 with 0.056 < x < 0.062 (Pratt et al. 


2011). These results have been hailed as direct evidence 


for spin-density-wave order in iron pinctides (Pratt et al. 


2011), where the static AF order arises from Fermi sur¬ 


face nesting between the hole and electron pockets at 
the r and M points of the reciprocal space, respectively 
(Dong et al. I 2008 Fink et al. 2009 Vorontsov et al. 


2009). 


Figure 5(d) shows the phase diagram of electron-doped 
BaFe 2 _ a; Ni 2 ,As 2 obtained from X-ray and neutron scat¬ 


2013 


Luo et al. 2012b). 


tering experiments (Lu et al. 

Here, the AF order decreases with increasing Ni-doping 
and disappears near optimal superconductivity in a first 
order like fashion with an avoided quantum critical point 
( Lu et ad. \ |2013 ). Similar to Ba(Fei a; Co a ;) 2 As 2 with 
0.056 < x < 0.062 ( |Pratt et al. 2011), there is short- 
range (60 A) incommensurate AF order for samples near 


Ba(Fe 1 . x Cr x ) 2 As 2 

SDW 
Structural Transitio 
G-typeAFNI 


PM 

V, (Tetragonal) / G-Typ 

^4 / AFM 

SDW A ■_ 

(Orthorhombic) 
u 0.0 0.1 0.2 0.3 0.4 0. 







FIG. 6 (Color online) The electronic phase diagrams of 
various other doped 122 family of iron pnictides. (a) 
The structural and magnetic phase diagram of Cr-doped 
Ba(Fei_ 2 ,Cr ;c ) 2 As 2 (Marty et al. 201 1| ). For x < 0.2, Cr 
doping suppresses the coupled structural and magnetic phase 
transition without inducing superconductivity. For x > 0.3, 
the system becomes a G-type antiferromagnet (see Fig. 1). 
(b) The structural and magnetic phase diagram for Mn- 
doped Ba(Fei_a,Mn x ) 2 As 2 (K im et al. | |201 0). (c) Similar 

phase diagram for Ba(Fei_ x Mn a ;) 2 As 2 obtained by another 
group (I nosov et al .112013 1. Here the system is believed to 
form a disordered spin glass (Griffiths) phase for x > 0.1. 
(d) The structural and magnetic phase diagram of isoelec- 
tronic doped BaFe 2 (Asi_ 2 ,P a; )2 determined from transport 
and NMR measurement ( |Shibauchi et al. 20141. Recent neu¬ 
tron powder diffraction experiments indicate a coupled struc¬ 
tural and magnetic phase transition (Allred et al. 2014). (e) 
The structural and magnetic phas e diagram of the is oelectron- 
ically doped Ba(Fei_ x Rua;) 2 As 2 {Kim et al. :2011). There is 
no evidence of a quantum critical point near optimal super¬ 
conductivity. (f) The pressure dependence of the structural 
and magnetic phase transitions in CaFe 2 As 2 . The system en¬ 
ters into a collapsed tetragonal (cT) p hase above ~ 0.4 GPa 
where magnetism disappears (Goldman et al. 2009). 


optimal superconductivity (Luo et al. 2012b). Although 


these results indicate an avoided quantum critical point 
in BaFe 2 - a; Ni a; As 2 , they are in direct conflict with a re¬ 
cent NMR work suggesting the presence of two quantum 
critical points (Zhou et al. 2013b). However, these NMR 


results are inconsistent with systematic NMR and neu¬ 
tron scattering results on nearly optimally Co and Ni- 
doped BaFe 2 As 2 samples revealing a cluster spin glass 


state for the magnetic order (Dioguardi et al. 2013 


Lu et al. 2014b). The spin glass picture of the mag¬ 


netic order near optimal superconductivity in electron- 
doped iron pnictides is clearly inconsistent with the spin- 
density-wave explanation of the transverse incommensu- 
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FIG. 7 (Color online) The doping evolution of the electronic 
phase diagrams for P-doped CeFeAsi-ajP^O and H-doped 
LaFeAsOi-^Ha;. (a) The structural and magnetic phase di¬ 
agram of CeFeAsi-^P^ O determined from neu tron powder 
diffraction experiments (de la Cruz et al., 2010). A quantum 


critical point is found near x = 0.4 without the presence of 
superconductivity. The inset shows the Tn from Ce magnetic 
ordering, (b) Phase diagram of CeFeAsi-^P^O determined 
by transport measurements. In the underdoped regime, the 
data is consistent with the results of neutron powder diffrac- 
2010|). The system becomes a Ce- 


tion (de la Cruz et al. 


ordered ferromagnetic metal for 0.4 < x < 0.9. For samples 


with x > 0.9, it becomes heavy fermion-like metal (Luo et al. 
[2010). (c) The structural, magnetic, and superconducting 

phase transitions in H-doped LaFeAsOi-a-H^ (Hiraishi et al. 


20141. There are two AF phases with different magnetic struc¬ 


tures near two superconducting domes. 


and superconducting coexistence region is also present. 
In particular, for Bai_ a ,Na. r Fe 2 As 2 near x = 0.28, 
a new magnetic ordered phase in the C 4 tetragonal 
symmetry of the underlying lattice has been found (Avci 
et al. 2014). In addition, superconductivity appears 


in heavily hole-doped regimes, much different from the 
electron-doped case. 

In copper oxide superconductors, superconductivity 
can only be induced by electron and hole doping into the 
nearly perfect Cu0 2 plane, and impurity substitution at 
the Cu sites by other elements dramatically suppresses 


superconductivity (Armitage et al, 

2010 

Fujita et al. 

2012; |Kastner et al., 1998; Kivelson et al. 

, |2003| |Tran- 

quada et al. 

2014 

). The situation is much different for 

iron pnictides. While impurities such as Cr and Mn sub- 


stituted for Fe in BaFe 2 As 2 suppress the static AF order 
in the parent compound without inducing superconduc- 
[Figs . 


tivity 


6 (a ), 6 (b) and 6 (c)] (Inosov et al. 2013 


Kim et al. 2010| [Marty et al. 2011), isoelectronic sub- 


stitution by replacing As wit h P |Fig. 6 (d)] (|Jiang et~al. 


2009 

Shibauchi et al. 

2014 

) or Fe with Ru [Fig. 6 (e)] in 

BaFe 2 As 2 ( 

Kim et al. 

2011 

) can induce superconductiv- 

ity. For Cr-doped Ba() 

H ^l— X 

Cr a ,) 2 As 2 , neutron diffraction 


experiments on single crystals have established the struc¬ 
tural and magnetic phase diagram, showing a suppres¬ 
sion of the collinear AF order for samples with x < 0.3. 
For x > 0.3, the system becomes a G-type antiferro- 
magnet with a tetragonal structure as shown in Fig. 
6 (a) (Marty et al. 2011). The situation in Mn-doped 


Ba(Fei_ x Mnj,) 2 As 2 is somewhat similar. With increas¬ 
ing Mn doping in BaFe 2 As 2 , the structural and AF phase 
transitions ar e gradually suppressed as shown in Figs. 
6 (b) and 6 (c) (Inosov et al. |2013 Kim et al. 2010). For 
Mn concentration x > 0.1, the system goes into a mixed 
phase, possibly in the Griffiths regime, with coexisting 
short-range spin excitations at AF wave vectors similar 
to those in BaFe 2 As 2 [Qaf = Qstripe] and BaMn 2 As 2 


[Q = QjveeZ rotated 45° from Qaf] (Inosov et al. 2013). 


In contrast to Cr and Mn doping, isoelectronic doping 
by replacing As with P in BaFe 2 (Asi_zIG ) 2 induces op- 


tirnal superconductivity near x = 0.3 (Jiang et al. 2009). 


From the systematic transport and London penetration 
depth measurements on BaFe 2 (Asi_ x P a ;) 2 , a quantum 
critical point has been identified near optimal supercon¬ 


rate magnetic order (Pratt et al. 2011). These results ductivity at x = 0.3 [Fig. 6 (d)] (Shibauchi et al. 2014). 


suggest that the incommensurate AF order in electron- For isoelectronic Ba^ei-aRu^^/optimal supercon- 


doped iron pnictides arises from a localized moment 


(Fang et al. 2008 Si and Abrahams 

2008 

Xu et al. 

2008 

), instead of spin-density-wave ore 

er induced from 


nested Fermi surfaces like incommensurate AF order in 


pure chromium metal (Fawcett et al. 1994). 

The electronic phase diagrams of 
Bai_ x K x Fe 2 As 2 and Ba 1 _ a: Na ;c Fe 2 As 2 


ductivity a gain appear s near Ru-doping level of x = 0.3 
[Fig. 6 (e)] (KimFFaT] 2011). However, there are no re¬ 
port for the presence of a quantum critical point in this 
system. Figure 6 (f) shows the pressure dependence of the 


hole-doped 
determined 

from neutron scattering experiments on powder samples 
are summarized in Figures 5 (a) and 5(c), respectively 
( |Avci et ~aL\ 2013 2014) 2012). Compared with electron- 
doped BaFe 2 As 2 , hole-doping does not separate the 
structural and magnetic phase transitions and the AF 


CaFe 2 As 2 phase diagram (Goldman et al. 2009). While 


superconductivity can be induced directly via ap plying 
hydrostatic pressure in BaFe 2 As 2 and SrFe 2 As 2 (John¬ 


ston 


2010 Stewart, 2011), external pressure exerted on 


CaFe 2 As 2 results in a nonmagnetic collapsed tetragonal 
(cT) phase, eliminating the static AF ordered moment 
and spin excitations without inducing superconductivity 
(|Goldman et al.\ 2009). 
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FIG. 8 (Color online) Spin waves in BaFe 2 As 2 and SrFe 2 As 2 
determined from neutron time-of-flight spectrometry, (a) 
Spin waves at E = 26±10 meV. The presence of peaks at wave 
vectors Qaf = (±1,0) and (0, ±1) is due to twinning effect. 
The white regions are detector gaps. Similar spin waves at (b) 


E = 81±10, (c) 157±10, and (d) 214±10 meV (Harriger et al. 


2011). The color bars indicate magnetic scattering in absolute 
units of mbarn-sr _1 -meV _1 -f.u. _1 . (e) Spin-wave dispersion 
curves and fits using a Heisenberg Hamiltonian with differ¬ 
ent exchange couplings along the [1, A'] direction, (f) Similar 


Heisenberg Hamiltonian fits along the [H, 0] direction (Har¬ 


riger et al. 20111. (g) Spin waves in the energy and wave 


vector cuts along the [0.5, K ] direction for SrFe 2 As 2 (|Ewings| 


et al. 2011). (h) The dashed line shows fit of a Heisenberg 
Hamiltonian assuming one set of exchange coupling constants, 
(i) A RPA calculation of y"(Q,u;) based on a 5-band model 
(Ewings et al. 2011). The reciprocal space notation in (g,h,i) 


is tetragonal where Qaf = (0.5, 0.5), different from those in 
(a-f). 


Although a majority of neutron scattering work has 
been focused on the 122 family of materials because of the 
availability of high-quality single crystals, there are also 
important phase diagram results in the 1111 family. For 
example, P-doping in the CeFeAsO family of materials 
suppresses static AF order, but doe s not induce supercon¬ 
ductivity (de la Cruz et al.. 2010). Systematic neutron 


scattering studies of the structural and magnetic phase 
transitions in powder samples of CeFeAsi-^P^O suggest 
that the pnictogen height [the average Fe-As(P) distance] 
and orthorhombicity of the CeFeAsi-^P^O unit cell crit¬ 
ically control the iron AF ordered moment and Neel tem¬ 
perature of the system. Figure 7(a) shows the P-doping 
dependence of the structural and AF phase transition 
temperatures in CeFeAsi-^P^O, suggesting the pres¬ 
ence of a magnetic quantum critical point near x = 0.4 
(de la Cruz et al. 2010| ). A complete mapping of the 
CeFeAsi-zP^O phase diagram shown in Fig. 7(b) was 
obtained via transport and susceptibility measurements, 
which reveal that superconductivity does not appear in 
the entire phase diagram, possibly due to heavy fermion 
properties of the rare earth element Ce (Luo et al. 2010). 


Another recent advance is the discovery of bipartite mag¬ 
netic parent phases in the H-doped LaFeAsOi-^Hj, fam¬ 
ily of materials (JHiraishi et al. 2014). In contrast to the 
general phase diagram of iron pnictides, superconductiv¬ 
ity in LaFeAsOi_ x Ha; appears in two domes adjacent to 
two different AF phases with different magnetic struc¬ 


tures and Neel temperatures [Fig. 7(c)] (Hiraishi et al. 


2014). These results again confirm the notion that super¬ 


conductivity in iron-based superconductors is intimately 
connected with the magnetic interactions. 


III. SPIN EXCITATIONS AND THEIR RELATIONSHIP 
WITH SUPERCONDUCTIVITY 


The rapid development of neutron time-of-flight chop¬ 
per spectrometers in recent years has allowed measure¬ 
ments of spin excitations in high-T c superconductors 
throughout the Brillouin zone for energy transfers up to 1 
eV. In the case of copper oxides, spin waves in La 2 CuC >4 


(Coidea et al. 

2001 

Headings et al. 

2010 


While the 

low energy spin excitations are well described by theory 
based on the Heisenberg Hamiltonian, high-energy spin 
waves are damped near the Q = (1/2,0) position in re¬ 
ciprocal space and merge into a momentum dependent 
continuum suggesting the decay of spin waves into other 


excitations (Coidea et al. 2001 Headings et al. 2010). 


The effective magnetic exchange couplings of La 2 Cu 04 
determined from the Heisenberg model are summarized 
in Table II. The doping evolution of spin excitations as a 
function of electron and hole doping and their coupling to 
superconductivity are reviewed recently dArmitage et al 


2010 Fujita et al.. 2012 Tranquada et al. 20141. In 


the case of iron-based superconductors, the situation is 
more complicated. Of the 11, 111, 122, 1111, 245 fami¬ 
lies of materials, spin waves in the AF parent compounds 


11 (I 

dpscombe et al. 

2011 Zaliznyak et al. 20T 

n, in 

(Zhang et al., 2014a 

, 122 (Diallo et al. 20091 1 

Dwings 

et al 

2011| Harriger et al, 2011| Zhao et al., 2009), and 

245 

Chi et al., 2013| Wang et al. 2011c Xiao et al. 

2013 

) families of materials due to the availability of large 


single crystals needed for inelastic neutron scattering ex¬ 
periments. Although single crystals of the 1111 family 
of materials are still not large enough to allow a deter¬ 
mination of the entire spin wave spectra, measurements 
of low-energy spin waves reveal that the system is highly 
two-dimensional with weak magnetic exchange coupling 
along the c axis (Ramazanoglu et al. 2013). In the Sec¬ 


tions III A and B, we describe spin wave measurements in 
the parent compounds of different families of iron-based 
superconductors and discuss their relationship with su¬ 
perconductivity. 
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A. Spin waves in the parent compounds of iron-based 
superconductors 

Inelastic neutron scattering studies of spin waves in 
the parent compounds of iron-based superconductors be¬ 
gan soon after the availability of single crystals of the 
122 family (Canfield and Bud’ko 2010). For a magnet¬ 


la) 1 5±3 meV 


ically ordered system, spin waves occur when the mag¬ 
netic moments precess around their ordered configura¬ 
tion. Regardless of the microscopic origin of the mag¬ 
netic order, spin waves of an ordered system should ex¬ 
hibit sharp excitations in the long wavelength limit (small 
Q) and can be described by a suitable Hamiltonian us¬ 
ing perturbation theory. For a spin Hamiltonian, one 
can starts with a Heisenberg model where the energy of 
spin waves depends only on the relative orientation of 
neighboring spins. In the initial neutron scattering ex- 


periments on low-energy spin waves in SrFe 2 As 2 

Zhao 

et al. 

2008a), CaFe 2 As 2 (McQueeney et al. 

2008 

, and 

BaFe 2 As 2 (Matan et al. 

2009), a spin gap due to 

mag- 


low-energy spin waves were described by either a lo¬ 


cal moment Heisenberg Hamiltonian (McQueeney et al. 


2008 


Zhao et al. 2008a I or the spin excitation contin¬ 


uum from itinerant electrons (Diallo et al. 2009 Matan 
et al. 2009). However, these measurements were unable 


to reach spin waves near the zone boundary and thus did 
not allow a conclusive determination of the effective near¬ 
est and next nearest neighbor magnetic exchange cou¬ 
plings denoted as Ji a /Jib and J 2 , respectively [Fig. 4(b)]. 
In the itinerant picture of the magnetism in iron pnictides 
(Hirschfeld et al. 2011[ Mazin 2010), spin waves from 
the AF ordered phase should arise from quasiparticle ex¬ 
citations between the electron and hole Fermi surfaces 
and form a spin excitation continuum at high-energies 
(Kaneshita and Tohyama 2010). In the initial neutron 


time-of-flight experiments on CaFe 2 As 2 , spin waves up to 
an energy of ~100 meV were measured and found to fit 


a Heisenberg Hamiltonian (Diallo et al. 2009). However, 


no spin wave signals were found for energies above 100 
meV consistent with ab initio calculations of the dynamic 
magnetic susceptibility, indicating that the high energy 
spin excitations are dominated by the damping of spin 


2009). 


waves by particle-hole excitations (Diallo et al. 

In subsequent neutron scattering experiments on 


CaFe 2 As 2 (Zhao et al. 2009), BaFe 2 As 2 (Harriger et al. 
2011), and SrFe 2 As 2 using more sample mass (Ewings 


et al. 2011 ), spin waves were mapped out throughout 


the Brillouin zone and the zone boundary energy scales 
were found to be around 220 meV. Figures 8(a)-8(d) show 
images of spin waves in BaFe 2 As 2 in the AF ordered state 
at energies of E = 26±10, 81±10, 157±10, and 214± 10 


meV, respectively (Harriger et al. 2011). With increas¬ 


ing energy, spin waves become diffusive but one can still 
see clear excitations near the zone boundary at E = 214 


(b) 45±3 meV 



0 ( 1 , 0 )( 1 , 1 ) ( 0 , 0 ) ( 1 , 0 )( 1 , 1 ) ( 0 , 0 ) 


FIG. 9 (Color online) Comparison of spin waves in NaFeAs 
and combined DFT and DMFT calculations, (a) Spin waves 
of NaFeAs at E = 15 ± 3 meV. Similar spin waves at (b) 
E = 45 ± 3; (c) 60 ± 10; (d) 100 ± 10 meV. The magnetic 

intensity are in absolute units of mbarn-sr -1 -meV -1 -f.u. -1 . 
(e) The solid circles show dispersion of spin waves, while the 
color plots are calculations from a combined DFT and DMFT 
theory, (f) Similar comparison for BaFe 2 As 2 (]Zhang et al .| 
2014a I. 


tions. Using a Heisenberg Hamiltonian with anisotropic 
spin wave damping, one can fit the entire spin wave spec¬ 
trum with a large in-plane nearest neighbor magnetic ex¬ 
change anisotropy (Ji Q > 0, Jp, < 0) and finite next 
nearest neighbor exchange co upling (J 2 > 0) (|Harriger 


et al. 2011 Zhao et al. 20091. The details of Heisenberg 


Hamiltonian for spin waves have been discussed in 

Diallo 

et al. 

2009 

Harriger et al. 

2011 

Zhao et al. | 2009 

). The 


exchanges are shown as solid lines in Fig. 8(e), while 
the dashed lines in the Figure are calculations assum¬ 
ing isotropic in-plane magnetic exchange couplings. The 
discovery of large in-plane exchange anisotropy is surpris¬ 
ing given the small orthorhombic lattice distortion in the 


AF ordered state (Wysocki et al. 2011). Only by prob¬ 


ing spin waves at high energies near the zone boundary, 
one can conclusively determine the effective magnetic ex¬ 
change couplings in the system. Different magnetic struc¬ 
tures and spin exchange couplings in iron-based materials 
has been studied using a localized moment model with 
different nearest and next nearest neighbor exchange cou- 


meV, different from the earlier experiment (Diallo et al. plings ( Hu et al. 2012). 


2009). Figures 8(e) and 8(f) show spin wave dispersions 
of BaFe 2 As 2 along the in-plane [1 ,K\ and [17,0] direc- 


Although spin waves in CaFe 2 As 2 (Zhao et al. |2009 ) 
and BaFe 2 As 2 (Harriger et al. 20111 can be modeled by 
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28.5 ±25 meV 



0 0.5 1 1.5 
h (r.I.u.) 

85 ± 15 meV 


115 ± 15 meV 


(h) a 


m 

■ 1.0 

■ 

_ 

:■ 

N 


:■ 


0 0.5 1 1.5 
h (r.I.u.) 


0 0.5 1 15 
h (r.I.u.) 


-1 -0.5 0 0.5 1 

K (r.I.u.) 

Q=(0,0.45,0-1.8) 
-o- 10K 
-o- 80K 
o 300K 




5 10 15 20 25 

E (meV) 


J la=-17.5; J1 b=-51.0; J2a=J2b=21.7; J 3=6.8 




FIG. 10 (Color online) Spin waves in Fei.osTe and FemTe. 
(a) Wave vector dependence of spin waves in Fei.osTe at 
E = 7.5 ±2.5 meV. Similar spin waves at (b) E = 28.5 ± 2.5, 
(c) 85 ± 15, (d) 115 ± 15 meV (Lipscombe et al. 2011| . The 
diffusive nature of spin waves is clearly seen at high ener¬ 
gies. (e,f) Dispersion curves of spin waves along the [H, 0] 
and [1, A'] directions and the Heisenberg Hamiltonian fits us¬ 
ing nearest, next nearest, and next next nearest neighbor ex¬ 
change couplings (Lipscombe et al. 20111. (g) Energy and 


wave vector cut of spin waves in Fei.iTe. (h) Energy depen¬ 
dence of spin waves at different temperatures, (i) Tempera¬ 
ture dependence of the integrated moments for Fei.iTe. The 
data suggests an increased total integrated moment on warm¬ 


ing from 10 K to 100 K across Tn and T s (Zaliznyak et al. 
2011 ). 


a local moment Heisenberg Hamiltonian, one still has 
to use anisotropic spin wave clamping characteristic of 
an itinerant electron system. In the neutron scattering 
work on spin waves of SrFe 2 As 2 [Fig. 8 (g)] (Ewings et al., 
2011), the authors report that a Heisenberg Hamiltonian 


that can fit the low-energy spin wave data fails to describe 
the spectrum near the zone boundary [Fig. 8 (h)]. The 
overall spin wave spectrum is instead best described by 
an itinerant model with large spin wave damping near 
the zone boundary [Fig. 8 (i)] (Ewings et al., 2011). 

Similar to the 122 family of materials, NaFeAs, the 
parent compound of the 111 family of iron pnictides, has 
the collinear AF structure albeit with a greatly reduced 


ordered moment size (Li et al. 2009c). Triple-axis neu¬ 


tron scattering experiments on single crystals of NaFeAs 
studied low-energy spin waves and found a small gap in 
the excitation spectrum (Park et al.l 2012 Song et al. 


2013b). Figure 9 summarizes the evolution of spin waves 


to the zone boundary as a function of increasing energy 
(Zhang et al. 2014a). Compared with the spin wave zone 
boundary energy of ~220 meV in BaFe 2 As 2 as shown in 
Fig. 8 , spin waves in NaFeAs reach the zone boundary 


at the in-plane wave vector Q = (1,1) around ~110 meV 
[Fig. 9(d)]. This means that the overall magnetic exci¬ 
tation bandwidth in the 111 family is considerably lower 
than that of the 122 family of iron pnictides. Figures 
9(e) and 9(f) compare the experimental and the com¬ 
bined density functional theory and dynamical mean field 
theory (DFT±DMFT) calculations of spin wave disper¬ 
sion of NaFeAs and BaFe 2 As 2 , respectively. The out¬ 
come suggests that pnictogen height is correlated with 
the strength of electron-electron correlations and conse¬ 
quently the effective bandwidth of magnetic excitations 
in iron pnictides (Yin et al . 2014 Zhang et al. 2014aI. 


Figure 10 summarizes spin wave measurements for the 


iron chalcogenide Fei +a; Te ( 

Bao et al. 

2009 Fruchart 

et al. 

1975 

Li et al., 2009d 

), the parent compound of 


the 11 fami 


y of iron-based superconductors (Lipscombe 


et al. 201 1[ Zaliznyak et al. 2011). The static magnetic 


order and spin excitations of Fei+^Te are sensitive to 


the excess iron in the interstitial sites (Rodriguez et al. 


2011 2010 


Stock et al. 2011 Wen et al. 2011). This 


is rather different from the iron pnictides, which cannot 
accommodate any excess iron in the crystal structure. 
For Fei.osTe and Fei.iTe, the AF structu re is commensu¬ 
rate bi-collinear (Rodriguez et al. 2011| ). Figures 10(a), 
10 (b), 10 (c), and 10 (d) show the two-dimensional images 
of spin waves in Fei.osTe at E = 7.5 ± 2.5, 28.5 ± 2.5, 


85±15, and 115± 15 meV, respectively (Lipscombe et al. 


2011). The dispersion of spin waves is different from that 


of the 122 and 111 families, and becomes diffusive for 
energies above 85 meV without well-defined spin waves 
[Figs. 10(c) and 10(d)]. The solid lines in Figs. 10(e) 
and 10(f) show fits of the dispersion using a Heisenberg 
Hamilton ian assuming exchange couplings Ji a , Jib, J 2 , 
and J 3 (Lipscombe et al. 2011). In a separate neutron 


scattering experiment (Za 


iznyak et al. 2011 ), the au¬ 


thors find that the low-energy spin excitations can be 
well-described by liquid-like spin plaquette correlations 
[Fig. 10(g)]. Furthermore, the integrated magnetic exci¬ 
tation intensity increases on warming [Fig. 10(h)]. The 
effective spin per Fe S « 1 at T ss 10 K in the AF ordered 
phase grows to S ~ 3/2 at T = 80 K in the paramagnetic 
phase, suggesting that the local magnetic moments are 
entangled with the itinerant electrons in the system [Fig. 
10(i)] ( |Zaliznyak et al. 2011). 

Of all the iron-based superconductors, alkali iron se- 


Fang et al. 

2011 

Guo et al. 

2010 ) 


icmuco e~l x ±. ^2—y' _ _ 

are unique in that superconductivity in this class of ma- 
terials always coexists with a static long-range AF order 


with a large moment and high Neel temperature (Bao 

| Wang et al.\ | 2011 b 


et al., 2011 

May et al. 

et al., 2011 


Zhao et al. 


2012 


Ye 


2 OI 2 J) . Although there is am¬ 
ple evidence indicating that the superconducting alkali 
iron selenides are mesoscopically phase separated from 
the insulating A 2 Fe 4 Se 5 phase with the y/5 x \/5 block 


AF structure as shown in Fig. 2(d) (jCarr et al., 2014 

Charnukha et al. 2012| Ksenofontov et al. \ 2011 Li et al. 

2011 bj Ricci et al. 

2011 Shermadini et al. 2012[ Shoe- 

maker et al. 2012 j 

Speller et al. 2012 Texier et al. 20121 
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[H 0 , H 0 ](r.l.u.) 


(a) 


(b) 



-1 

1 

0 * 

-1 

1 

0 

-1 




0.8 1 1.2 
H(r.l.u.) 


0.8 1 1.2 
H (r.I.u.) 


FIG. 11 (Color online) Spin waves in the insulating 245 and 
semiconducting 234 phases. (a,b,c) Spin-wave dispersions in 
the insulating Rbo.ggFei. 58 Se 2 (Wang et al 2011c). The solid 


circles are data from the cut directions marked by the ar¬ 
rows in the right panels of (a,b,c). The solid lines are fits 
from a Heisenberg Hamiltonian considering nearest, the next 
nearest, and next next nearest neighbor exchange couplings 
( |Wang et al. 201 lc|. (d) The energy dependence of the local 
dynamic susceptibility, where the solid line is the calculated 
value from the Heisenberg Hamiltonian. The total moment 
sum rule appears to be satisfied (Wang et al. 2011c). (e,f,g,h) 
The wave vector dependence of spin waves from the semicon¬ 


ducting Ko. 85 Fei. 54 Se 2 (Zhao et al. 2014). 


Wang et al. 2012cI, there is still no consensus on the 


chemical and magnetic structures of their parent com- 


pounds (May et al. 2012 

Wang et al.\ 

2011b 

Ye et al. 

2011 

Zhao et al. 2012 

). Assuming that the insulating 


H. 2 Fe 4 Se 5 phase is the parent compound of the supercon¬ 
ducting yt a: Fe 2 -ySe 2 , its spin waves have been mapped 
out by several groups (Chi et al. \ 2013] Wang et al. 2011c 


Xiao et al. 2013). Compared with spin waves in iron 


hole-doped 


0 1 

electron-doped 


Q = (l,0) o = (1/2, 1/2 ) t 
H (r.I.u.) 




V 


y+ 


0 10 20 30 40 50 60 70 

Temperature (K) 



FIG. 12 (Color online) Neutron spin resonance in electron- 
and hole-doped iron pnictides. (a) The schematic drawings 
of the wave vector dependence of the low-energy spin exci¬ 
tations in optimally hole- (upper panel) and electron-doped 
(lower panel) superconducting iron pnictides. (b) Tempera¬ 
ture dependence of the resonance at E = 16 meV, showing 
clear superconducting order parameter-like enhancement be 
low T c for a powder sample of Bao.6Ko.4Fe2As2 (Christianson 


et al. 2008). (c) Temperature dependence of the magnetic 


scattering at energies E = 3.0, 9.5, and 16 meV for optimally 
electron-doped superconducting BaFei 85 C 00 15 AS 2 (T c = 25 
K) ( |Inosov et ffll.| [2010). While magnetic intensity at the res- 
onance energy (E = 9.5 meV) shows a clear enhancement 
below T c at the expense of opening a spin gap at E = 3.0 
meV, the scattering at E = 16 meV is not sensitive to super¬ 
conductivity. (d) Temperature dependence of the resonance 
energy for BaFei. 85 C 00 . 15 AS 2 (Inosov et al. 2010). 


sion of the spin waves in insulating A 2 Fe 4 Ses has two op¬ 
tical branches at high energies and one acoustic branch at 
low energy, where the arrows are wave vector scales and 
the thin dashed line separates the vertical energy scale for 
the acoustic and low-energy optical spin waves from the 
high-energy optical spin waves [Figs. ll(a)-ll(c)] (Wang 


pnictides and iron clialcogenides (Figs. 7-10), the disper- 


et al. 2011c). By integrating the energy dependence of 
the local dynamic susceptibility in Rbo. 89 Fei. 5 gSe 2 [Fig. 
11(d)], it was found that the total moment sum rule is 
exhausted for magnetic scattering at energies below 250 
meV. Therefore, spin waves in insulating Rbo.g 9 Fei. 5 gSe 2 
can be regarded as a classic local moment system where 
a Heisenberg Hamiltonian is an appropriate description 
of the spin wave spectrum. 

On the other hand, if the semiconducting AF phase 
with rhombus iron vacancy order [Fig. 2(e)] is the par¬ 
ent compound ( Wang et al. 2014| Zhao et al . | 2012), one 
finds that spin waves of the system are rather close to 
those of iron pnictides. Figure 11(e), 11(f), 11(g), 11(h) 
shows the evolution of spin waves as a function of increas¬ 
ing energy for the semiconducting Ko. 8 sFei. 54 Se 2 with 












































































































14 



Temperature (K) 



FIG. 13 (Color online) Effect of electron-doping on mag¬ 
netism and superconductivity in electron underdoped iron 
pnictides. (a) Temperature dependence of the nuclear (2, 2,0) 
and (1/2,1/2,1) (in tetragonal notation) magnetic scattering 
in the electron underdoped Ba(Feo. 953 Coo.o 47 ) 2 As 2 (T c = 17 
K). The structural, magnetic, and superconducting transi¬ 
tions are clearly marked, (b) A weak resonance appears below 
T c at E = 4 meV (Pratt et al. 2009). 


collinear AF order and TV = 280 K (Zhao et al. 20141. 


The data agrees well with calculations using a Heisen¬ 
berg Hamiltonian. A comparison of the observed spin 
wave spectrum in this system with those of the CaFe 2 As 2 


single crystals (Zhao et al. 2009) reveals remarkable sim¬ 
ilarity, and thus suggesting simi ar effective magnetic ex¬ 
change couplings in these systems (jZhao et al. 2014). 

Table II summarizes the effective magnetic exchange 
couplings for the parent compounds of known iron-based 
superconductors. We also list the effective magnetic ex¬ 
change couplings for La 2 Cu 04 , the parent compound of 
copper oxide superconductors. They are dominated by 
the large nearest neighbor and weak next nearest neigh¬ 
bor magnetic exchange couplings (Fig. 4). For the par¬ 
ent compounds of iron-based superconductors, it is in¬ 
structive to compare their effective magnetic exchange 
couplings. In spite of their dramatically different AF 
structures summarized in Figs. 1-3 , they all appear to 
have similar next nearest neighbor magnetic exchange 
couplings (see Table II). This is consistent with the idea 
that the next nearest neighbor coupling J 2 is mainly de¬ 
termined by a local superexchange mechanism mediated 
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FIG. 14 (Color online) The evolution of Fermi surfaces in 
electron and hole-doped BaFe 2 As 2 . (a) Schematics of Fermi 
surfaces corresponding to BaFe 2 As 2 with possible nesting vec¬ 
tors marked with arrows ([Dai et al. ][2012|. The d xz , d yz , and 
d xy orbitals for different Fermi surfaces are colored as red, 
green and blue, respectively, (b) Fermi surfaces when 10% 
electrons are doped into BaFe 2 As 2 to form optimal supercon¬ 
ductivity. (c) Fermi surfaces with 30% electron doping when 
superconductivity is suppressed ( |Wang et al. | [2013b|. (d) 

Schematics of Fermi surfaces for hole-doped Bai_a;K x Fe 2 As 2 
with increasing K-doping to x = 0.23 (upper left panel), 0.40 
(upper right panel), 0.65 (lower left), and 0.86 (lower right 
panel) ( |Richard et al . | |2011 1. (e) A comparison of the Fermi 
surfaces for x = 0.4 and 1 in the folded Brillouin zone (|Richard| 

TTai] [20ll l. 


by As or Se/Te, regardless of their metallic or insulat¬ 
ing ground states (|Abrahams and Si[ |2011| |Hu and Ding 
2012 ]). 


B. Neutron spin resonance and its relationship with 
superconductivity 

The neutron spin resonance is a collective magnetic 
excitation occurring below T c with a temperature de¬ 
pendence similar to the superconducting order param¬ 
eter (Eschrig 20061. First discovered in hole doped 


YBaaC^Oe+a; copper oxide superconductors (Rossat- 


Mignod et al. 1991), the resonance is located near the 


AF ordering wave vector Qaf of the nonsuperconduct¬ 
ing parent compound and occurs at an energy related to 


the superconducting T c (Dai et al. 2000 Wilson et al. 


2006a) or gap energy (Yu et al. 2009). It has been ar¬ 


gued that the mode is a signature of the d-wave pairing 
as a result of quasiparticle excitations between the sign 


reversed d -wave superconducting gaps (Eschrig 2006). 


Soon after the discovery of iron pnictide superconduc¬ 


tors (Kamihara et al. 2008), a neutron spin resonance 


was found in powder samples of Bao. 6 Ko. 4 Fe 2 As 2 (Chris¬ 


tianson et al. 2008). Since the resonance occurs below 
T c at the momentum transfer (Q = 1.15 A ) close to 


Qaf i n BaF e 2 As 2 [Fig. 12 (a) and 12(b)] ( |Christianson 
et al. 20081, the mode is believed to arise from the sign 


reversed quasiparticle excitations between the hole and 
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TABLE II Comparison of the effective magnetic exchange couplings for parent compounds of copper-based and iron-based 
superconductors. Here the nearest, next nearest, next next nearest neighbor, and c axis exchange couplings are SJi a (SJib), 
SJ2a(SJ2b), SJ3, and SJ C , respectively, where S is the spin of the system. 


Materials 

SJia (meV) 

SJu (meV) 

SJ 2a (meV) 

SJ 2b (meV) 

SJ 3 (meV) 

SJ C (meV) 

La 2 Cu0 4 0 

55.9 ±2 

55.9 ±2 

-5.7 ± 1.5 

-5.7 ± 1.5 

0 

0 

NaFeAs b 

40 ± 0.8 

16 ±0.6 

19 ±0.4 

19 ±0.4 

0 

1.8 ±0.1 

CaFe 2 As 2 c 

49.9 ±9.9 

-5.7 ±4.5 

18.9 ±3.4 

18.9 ±3.4 

0 

5.3 ±1.3 

BaFe 2 As 2 d 

59.2 ±2.0 

-9.2 ± 1.2 

13.6 ± 1 

13.6 ± 1 

0 

1.8 ±0.3 

SrFe 2 As 2 (L) c 

30.8 ± 1 

-5 ±4.5 

21.7 ±0.4 

21.7 ±0.4 

0 

2.3 ±0.1 

SrFe 2 As 2 (H) f 

38.7 ±2 

-5 ±5 

27.3 ±0.3 

27.3 ±0.3 

0 

2.3 ±0.1 

Fei. 05 Te 9 

-17.5 ±5.7 

-51.0 ±3.4 

21.7 ±3.5 

21.7 ±3.5 

6.8 ±2.8 

~1 

Rbo.8gFei.58Se 2 h 

-36 ±2 

15 ±8 

12 ±2 

16 ±5 

9 ± 5 

1.4 ±0.2 

(Tl,Rb) 2 Fe 4 Se 5 ; 

-30 ± 1 

31 ±13 

10 ± 2 

29 ±6 

0 

0.8 ±1 

Ko.85Fei.5 4 Se 2 9 

-37.9 ±7.3 

-11.2 ±4.8 

19.0 ±2.4 

19.0 ±2.4 

0 

0.29 ±0.06 


Coidea et al. 


ZKaiig~~et al . 


Zhao et al. 


I] [2001 

:. -iliila . 

W 


I- 


(Harrier et al. ^Dllb. 
(iiiwings et a/.||zUll[). 


^'i'he L and H are fits using low and high-energy spin 

sp ectively. __ 

g ( Lipscombe et al. 2011}. 


I Wang et al. 2011c). 


“(Chi et al. 2013). 


h’Zhao et al. 2014 1 . 


waves, re¬ 


electron Fermi surfaces near the T and M points in recip- 2010). 


rocal space, respectively (Figs. 12-14) (Hirschfeld et al. 
2011 Mazin 2010|). In subsequent inelastic neutron scat 


tering experiments on single crystals of electron doped 


2009 


Ba(Fei_ x Coa,) 2 As 2 (Inosov et al. 2010 Lumsden et al. 


2009 ) and BaFe 2 -a;Ni a ,As 2 superconductors (Chi et al. 
|Li et al. \ |2009b|), the resonance was indeed founc 


at the in-plane AF ordering wave vector Qaf = (1, 0) 


of LaFeAsOi x F x 1 

Ishikado et al., 

2009 

Shamoto et al. 

2010 

Wakimoto et al. 

2010) and molecular-intercalated 

FeSe 1 

Taylor et al. 2011 

l) also revealed resonance like spin 


pendence of the imaginary part of the dynamic suscepti¬ 
bility ^"(Qaf,^) for energies below (E = Tito = 3 meV), 
at (E = 9.5 meV), and above (E = 16 meV) the res¬ 
onance in superconducting BaFei.gsCoo.isAsg (T c = 25 
K). It is clear that the intensity gain of the resonance be¬ 
low T c is at the expense of opening a spin gap at energies 
below it. By carefully monitoring the temperature de¬ 


pendence of the resonance, the authors of (Inosov et al. 


2010) suggest that the energy of the mode decreases with 


increasing temperature and may be directly correlated 
with the temperature dependence of the superconducting 
gap energy [Fig. 12(d)]. However, recent experiments on 
the nearly optimally doped BaFei.go 4 Nio.o 96 As 2 super¬ 
conductor found that the resonance energy is essentially 
temperature independent on warming (Luo et al. 2013aI, 
different from that of BaFei.ssCop.isAsg Qlnosov et al. 


In the electron underdoped regime where static AF or¬ 
der coexists and competes with superconductivity [Figs. 
5(b) and 5(d)], the static AF order occurs at a lower 
temperature than T s . Figure 13(a) shows the tempera¬ 
ture dependence of the nuclear peak intensity at ( 2 , 2 , 0 ) 
and magnetic Bragg scattering at Qaf for underdoped 
Ba(Feo. 953 Coo. 047 ) 2 As 2 (T c = 17 K) (Pratt et al. 2009). 
In the high temperature tetragonal state, the observed 
neutron scattering intensity from a strong nuclear Bragg 
peak ( 2 , 2 , 0 ) is lower than that expected from the struc¬ 
ture factor calculation due to multiple scattering ef¬ 


fect, termed as the neutron extinction effect (Hamilton 


1957). When the symmetry of the system is reduced from 


tegragonal to orthorhombic, there is a dramatic intensity 
gain below ~60 K arising from the release of the neutron 
extinction effect and the AF order occurs below Tjv ~ 48 
K. Upon entering into the superconducting state, the in¬ 
tensity of the static AF order decreases. Simultaneously, 
a weak neutron spin resonance appears at E = 4 meV 
[Fig. 13(b)], suggesting that the intensity gain of the 
mode arises from suppression of the static AF order and 


spin excitations at energies below the resonance (Chris¬ 


tianson et al. 2009 Pratt et al. 20091. Application of 


a magnetic field which partially suppresses superconduc¬ 
tivity will enhance the intensity of the static AF order 
and suppress the resonance ( Wang et al] 2011d). These 
results further suggest that the static AF order coexists 
and competes with superconductivity in electron under- 
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H (r.I.u.) H (r.I.u.) H (r.I.u.) H (r.I.u.) 


FIG. 15 (Color online) Comparison of wave vector evo¬ 
lution of the low-energy spin excitations in electron-doped 
BaFe 2 -^Ni ;r Aso with the RPA calculation based on a rigid 
band shift. (a,b,c,d) RPA calculation results obtained for an 
energy E = 8 meV for electron dopings of x = 0,0.07,0.1, 
and 0.15, respectively. As the doping increases from (a) to 
(d), one clearly sees an enhancement of the anisotropy in spin 
excitations (transverse elongation). (e,f,g,h) In-plane wave 
vector dependence of the spin excitations at E = 8 meV for 
x = 0,0.065,0.1, and 0.15, respectively (Luo et al. 2012b). 
For the electron overdoped x = 0.15 sample, two transverse 
incommensurate peaks are expected from the RPA calcula¬ 
tion. This is indeed observed in neutron scattering experi¬ 
ments (Luo et al. 2013a). 


doped iron pnictides. 


From density functional theory calculations (Kuroki 


et al. 


2008 Mazin et all 2008) and angle resolved pho¬ 


toemission spectroscopy (ARPES) experiments on elec¬ 


tron/hole doped iron pnitides (Chen et al. 2014 | Richard 


et al. 2011), we know that Fermi surfaces in most of these 


materials are composed of hole-like pockets near F and 
electron-like pockets near M point at Qaf = (1, 0). The 
neutron spin resonance in iron pnictides at Qaf = (1, 0) 
can arise from the sign reversed quasiparticle excitations 
between the hole and electron Fermi surfaces in an s" 1 - 


wave superconductor as shown in Fig. 14 (Hirschfeld 


et al. 2011 Mazin 2010), exhibiting the same signature 


as the sign changed su perco ndu cting gap function in the 
d -wave copper oxides ( |Eschri g| 2006). With increasing 
electron-doping, the hole and electron Fermi surfaces de¬ 
crease and increase in size, respectively [Figure 14(a)- 
14(c)]. Similarly, the hole Fermi pockets at the T point 
increase in size with increasing hole-doping, while the 
electron Fermi surfaces exhibit a Lifshitz transition at 
M point before becoming hole overdoped KFe 2 As 2 [Fig. 
14(d)-14(f)] ( |Chen et al.\fM4§ . 

Using the random phase approximation (RPA) based 
on a three-dimensional tight-binding model in the lo¬ 
cal density approximation (LDA) (Graser et al.\ 2010), 
calculations can predict the momentum anisotropy of 


the low-energy spin excitations and the resonance (Park 


et al. 2010). For electron-doped BaFe 2 - a; T' a ;As 2 , low- 


energy spin excitations become progressively elongated 
ellipses along the transverse direction relative to the spin 
waves in BaFe 2 As 2 due to the enhancement of the intra- 



FIG. 16 (Color online) The K-doping evolution of the 
neutron spin resonance and low-energy spin excitations in 
Bai_ 2 ,K a; Fe 2 As 2 . (a) Inelastic neutron scattering experiments 
on superconductivity-induced low-energy spin excitations of 
powder samples in Bai_ I ,K a! Fe 2 As 2 at x = 0.3. A clear res¬ 
onance is seen around 12 meV and Q = 1.25 A^ 1 as shown 
in (c). (b,d) The resonance becomes broader in Q at x = 0.5, 
and splits into two peaks at different wave vectors at (e,g) 
x = 0.7 and (f,h) x = 0.9 due to the changing Fermi sur¬ 


faces (Castellan et al. 201 1[). (i) Longitudinal scans along 


the [H, H, 1.3J direction above T c at E = 8 meV for single 
crystals of KFe 2 As 2 - Two incommensurate peaks are seen. 
The inset shows similar scan below T c at E = 3.5 meV. Lon¬ 


gi tudinal (j ) and transverse (k) scans at E = 6 meV (Lee 
et al. 2011). 


orbital, but inter-band, pair scattering process between 
the d xv orbitals [Figs. 14(a) and 14(b)] (Zhang et al. 


2010). Figure 15 shows the comparison of the RPA cal- 


culations and experimentally measured in-plane spin ex¬ 
citation anisotropy in BaFe 2 - a; Ni a ;As 2 superconductors 
(Luo et al. 2012a), confirming that the quasiparticle ex¬ 
citations between the hole and electron Fermi surfaces 
are consistent with the wave vector evolution of the low- 
energy spin excitations (Luo et al. 2013a). 

In the case of hole-doped materials, RPA calculations 
have predicted that spin excitations should be longitudi¬ 
nally elongated, and thus rotated 9(P from tho se of the 
electron-doped BaFe 2 _ a; T a ;As 2 (P ark etaTj 2010). Inelas¬ 
tic neutron scattering experiments on hole-doped single 
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FIG. 17 (Color online) The wave vector and energy depen¬ 
dence of the neutron spin resonance for BaFe 2 (Asi_ a: P ;r )2 
and NaFei-^CozAs. (a) The energy dependence of the dy¬ 
namic susceptibility at the in-plane AF wave vector position 
and different L values above T c for BaFe 2 (Asi_ :E P ; r )2 with 
x = 0.34. (b) Identical scans at temperature T = 5 K well 
below T c . The resonance shows clear dispersion for different 
values of L (Lee et al. 2013aI. (c) The energy dependence 


of the superconductivity-induced double resonance for an un¬ 
derdoped (UD) NaFei-xCozAs with x = 0.015. There are 
two peaks in the energy scan at E = 3.5 and 6 meV. (d) 
The double resonance in the underdoped sample becomes a 
single resonance for electron overdoped NaFei-xCo^As with 
x = 0.045 ( |Zhang et aZ.[[2013c[ ). 


crystals of superconducting Ba 0 . 67 Kao. 33 Fe 2 As 2 (T c = 38 
K) reveal longitudinally elongated spin excitations for 
energies near the resonance, consistent with RPA cal¬ 
culations (Zhang et al. |2011a ). Figure 16(a)-16(h) plots 
the hole-doping dependence of the resonance obtained 
using high-quality powder samples of Bai_ ;E K; C Fe 2 As 2 
(Avci et al. 2011 Castellan et al. 2011). Although these 


measurements do not provide precise information con¬ 
cerning the wave vector dependence of the spin excita¬ 
tions, they do give the hole-doping evolution of the total 
momentum transfer of the mode. With increasing hole¬ 
doping, the sharp resonance centered at Q ss 1.25 A -1 
for x = 0.3 [Fig. 16(a) and 16(c)] becomes broader in 
Q and splits into two peaks for x = 0.7 and 0.9 [Fig. 
16(e)-16(h)] ( |Castellan et al.\ 2011). This is consistent 
with the RPA result that hole-doping induces longitu- 


tution is much less effective in suppressing AF order and 
inducing superconductivity. Inelastic neutron scatter¬ 
ing experiments on BaFe 2 _ ;i: Ru :r As 2 near optimal super¬ 
conductivity reveal a neutron spin resonance similar to 
electron-doped BaFe^xT^As^ but with greatly damped 
intensity, possibly due to the weakening of the electron- 


electron correlations by Ru doping (Zhao et al. 2013). 
In the case of BaFe 2 (Asi_ :r P :E ) 2 , initial neutron scat¬ 
tering experiments on powder samples with T c = 30 K 
have revealed the pre s ence of a resonance at E ss 12 

Figure 17(a) and 17(b) 


meV (Ishikado et al. 2011). 


shows the energy dependence of x"(Q- UJ ) above and 
below T c , respectively, obtained for single crystals of 


BaFe 2 As 1 . 32 P 0.68 (T c = 29.5 K) (Lee et al. 2013a). In 
the normal state, x"(Q, w) is featureless and changes only 
slightly at different momentum transfers along the c axis 
(L = 0,0.25,0.5,0.75,1). Upon entering into the super¬ 
conducting state, a neutron spin resonance is formed and 
its energy is significantly dispersive along the c axis [Fig. 
17(b)] (Lee et al. 2013a). Since the bandwidth of the dis¬ 
persion becomes larger on approaching the AF ordered 
phase, the dispersive feature may arise from the three- 
dimensional AF spin correlations in the undoped parent 


(Lee et al. 2013a). 


So far, most of the neutron scattering work has been fo¬ 
cused on single crystals of electron/hole-doped BaFe 2 As 2 
family of materials. For the NaFei-^Co^As family of 
materials (Parker et al. ] |2010| |Tan et alj |2013| ) , the air 
sensitive nature of these materials makes it very diffi¬ 
cult to perform inelastic neutron scattering experiments 
(Tanatar et al. 2012). By using hydrogen free glue to 


coat the samples, neutron scattering experiments can be 
carried out to study the evolution of spin excitations in 


NaFei-^Co^As ( 

Park et al. 

2012 

Song et al. 

2013b). 

From ARPES experiments ( 

Ge et al. 

2013 

Liu et al. 


2011 ), it was found that the superconducting gap in the 


electron Fermi pockets of the underdoped regime near 
x = 0.0175 has a large anisotropy, which is absent in the 
hole Fermi pocket. The superconducting gap anisotropy 
disappears upon increasing x to 0.045. Figure 17(c) and 
17(d) shows the intensity gain of the resonance below 
T c for underdoped x = 0.015 (Zhang et al.. 2013c) and 


overdoped x = 0.045 (Zhang et al. 2013a), respectively. 


Instead of a single resonance peak, superconductivity in¬ 
duces a sharp resonance at E r 1 = 3.25 meV and a broad 


dinal incommensurate spin excitations {Castellan et al. resonance at E r2 = 6 meV ( |Zhang et al.\ |2013 C p . Similar 
20nl). Indeed, neutron scattering experiments on hole measurements on electron overdoped x = 0 045 reveal 


overdoped KFe 2 As 2 found two incommensurate spin exci¬ 
tation peaks located longitudinally away from Qaf [Fig. 
16(i)-16(k)], again confirming the notion that low-energy 
spin excitations in hole and electron-doped iron pnictides 
are controlled by quasiparticle excitations between the 


only one sharp resonance (Zhang et al. 2013a). The ap¬ 


pearance of the double resonance and the superconduct¬ 
ing gap anisotropy in the underdoped sample was inter¬ 
preted as originating from cither the orbital dependence 


of the superconducting pairing (Yu et al. 2014 Zhang 


hole and electron Fermi surfaces (Lee et al. 2011). 


et al. 2013c) or superconductivity coexisting with static 


In addition to electron or hole-doping to BaFe 2 As 2 , 
isoelectronic substitution to BaFe 2 As 2 by replacing Fe 


AF order in the iron pnictides (Lv et al.. 2014 Rowe 


et al. 2012). 


with Ru (Thaler et al. 2010) or As with P (Jiang et al. 


2009) can also induce superconductivity. Compared with 


the electron-doped BaFe 2 _ x T a; As 2 , isoelectronic substi¬ 
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FIG. 18 (Color online) The evolution of spin excitation dis¬ 
persions for hole and electron-doped BaFe 2 As 2 - (a) The elec¬ 
tronic phase diagram of electron and hole-doped BaFe 2 As 2 , 
where the arrows indicate the doping levels of inelastic neu¬ 
tron scattering experiments. The right inset shows the crystal 
and AF spin structures of BaFe 2 As 2 . The inset above x e = 0.1 
shows the transversely elongated ellipse representing the low- 
energy spin excitations in electron-doped BaFe 2 - a; Ni a ;As 2 in 
the ( H , K) plane of reciprocal space. The left insets show 
the evolution of low-energy spin excitations in hole-doped 
Bai_*K l; Fe 2 As 2 in the (H, K) plane. C-SF and IC-SF in¬ 
dicate commensurate and incommensurate spin fluctuations, 
respectively, (b-e) The solid lines in the figure are spin wave 
dispersions of the undoped BaFe 2 As 2 along the two high- 
symmetry directions. The symbols in (b), (c), (d), and 
(e) are dispersions of spin excitations for BaFe 1 . 9 Nio. 1 As 2 , 
BaFe 1 . 7 Nio. 3 As 2 , Bao.67Ko.33Fe2As2, and KFe 2 As 2 , respec¬ 
tively. The shaded areas indicate vanishing spin excitations 
(Wang et al. 2013bI. 


C. The electron and hole-doping evolution of the spin 
excitations in the BaFe2As 2 family of iron pnictides 

To understand the interplay between magnetism and 
superconductivity in iron pnictides, one must first deter¬ 
mine the electron and hole-doping evolution of the spin 
excitation spectra throughout the Brillouin zone. Since 
single crystals of electron and hole-doped BaFe 2 As 2 are 
available, one can systematically map out the evolution 
of spin excitations at different electron/hole-doping lev¬ 
els marked with arrows in the phase diagram [Fig. 18(a) 


(Chen et al. 

2014 

Harriger et al. 2011} Lee et al.\ 2011 

Liu et al. 2012a[ Luo et al, 2013a) Tucker et al. | 2012a 

Wang et al. 

2013b 

). The solid lines in Figs. 18(b)-18(e] 


show the dispersion of spin waves in BaFe 2 As 2 along 


the [1 ,K] and [H, 0] directions (Harriger et al. 2011). 
Upon electron-doping to induce optimal superconduc¬ 
tivity, spin excitations become broader at low-energies 
(E < 80 meV) and couple to superconductivity via the 
resonance while remaining almost unchanged at high en¬ 
ergies ( E > 80 meV) ( Liu et cd\ 2012a Tucker et al. 


2012a[). The red circles and yellow upper triangles in 


Fig. 18(b) show spin excitation dispersions of the opti¬ 
mally electron doped BaFe 1 .gNio. 1 As 2 at T = 5 K and 
150 K, respectively (Liu et al. 2012a). Figure 18(c) 


shows the dispersions of spin excitations of the electron- 
overdoped nonsuperconducting BaFe 1 . 7 Nio. 3 As 2 , where 
a large spin gap forms for energies below ^50 meV 
( |Wang et al. 2013b). Figure 18(d) and 18(e) shows the 
dispersions of spin excitations for optimally hole-doped 
Bao. 67 Ko. 33 Fe 2 As 2 and hole-overdoped KFe 2 As 2 , respec¬ 


tively (Wang et al. 


2013b). While electron doping does 


not much affect the high-energy spin excitations and dis¬ 
persion, hole-doping suppresses the high-energy spin ex¬ 
citations. 

Figure 19 reveals the evolution of the two-dimensional 
constant-energy images of spin excitations in the ( H,K ) 
plane at different energies as a function of electron-doping 


for BaFe 2 _ rE Nl r As 2 (Harriger et al. 

2011 } 

Luo et al., 

2013a 

Wang et al. 

2013b). In undoped BaFe 2 As 2 , there 


is an anisotropy spin gap below ~15 meV, thus there 
is essentially no signal at E = 9 ± 3 meV [Fig. 19(a)] 
(Matan et al. 20091. For nearly optimally electron doped 


x = 0.096, the spin gap is suppressed and low-energy spin 
excitations are dominated by the resonance [Fig . 19(f) 
( Chi et al. 2009| Inosov et al. 2010 Li et al] |2009b 


Lumsden et al.\ |2009| Luo et al. 201 2 a[). In electron 


overdoped BaFe 2 - x Ni x .As 2 with x = 0.15 ( T c = 14 K) 
and 0.18 (T c = 8 K), spin excitations at E = 8 ± 1 meV 
become weaker and more transversely elongated [Figs. 
19(k) and 19(p)] (Luo et al., 2013a| ). For the nonsuper¬ 
conducting x = 0.3 sample, a large spin gap forms in 
the low-energy excitation spectra [Fig. 19(u)]. Figures 
19(b)-19(e), 19(g)-19(j), 19(q)-19(t), 19(v)-19(y) show 
the evolution of spin excitations at different energies for 
BaFe 2 -x-Ni a; As 2 with x = 0,0.096,0.15,0.18, and 0.30, 
respectively. While electron doping modifies spin excita¬ 
tions at energies below E = 96 ± 10 meV, high energy 
spin excitations remain similar and only soften slightly. 

Figure 20 shows the constant-energy images of spin ex¬ 
citations as a function of hole-doping. For pure KFe 2 As 2 , 
incommensurate spin excitations along the longitudinal 
direction are seen at E = 8 ± 3 meV [Fig. 20(a)] 
and 13 ± 3 meV [Fig. 20(b)] (Lee et al. 2011). How¬ 
ever, spin excitations become much weaker at E = 


53 ± 8 meV (Wang et al. 2013b). For optimally hole- 


doped Bao. 67 Ko. 33 Fe 2 As 2 , the low-energy spin excita¬ 
tions change from longitudinally elongated ellipses at 
E = 5 ± 1 meV [Fig. 20(d)] to transversely elongated 
ellipses at E = 50 ± 2 meV [Fig. 20(f)]. At the neu¬ 
tron spin resonance energy oi E = 15 ± 1 meV, spin 
excitations change from longitudinally elongated ellipses 
above T c (not shown) to isotropic circles below T c in re- 
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x=0.0 x e = 0.096 x„=0.15 x e =0.18 x„=0.30 


(d) 96±10 meV (i) 96±10meV (n) 96±10meV (s) 96±10meV (x) 96±10meV 
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FIG. 19 (Color online) Constant-energy slices through mag¬ 
netic excitations of electron-doped BaFe 2 - :r Ni a ;As 2 iron pnic- 
tides at different energies. The color bars represent the 
vanadium normalized absolute spin excitation intensity in 
the units of mbarn-sr _1 meV _1 f.u. . (a-e) Spin waves in 

BaFe 2 As 2 at excitation energies of E = 9 ± 3, 19 ± 5, 60 ± 10, 
96 ± 10, and 180 ± 10 meV (Harriger et al. 20111. Spin 
waves peak at the AF ordering wave vectors Qaf = (±1,0) 
in the orthorhombic notation. Spin waves are also seen at 
Qaf ~ (0,±1) due to the twin domains of the orthorhom¬ 
bic structure, (f-j) Two-dimensional images of spin excita¬ 
tions for BaFei.904Nio.o96As2 at E = 8 ± 1, 16 ± 2, 60 ± 10, 
96 ± 10, 181 ± 10 meV. Identical slices as that of (f-j) for (k- 
o) BaF ei.85Nio.i5As2 and (p-t) BaFei.g 2 Nio.isAs 2 ( Luo et ai\ 
|20 12a| ). (u-y) Constant-energy slices through magnetic ex¬ 
citations of electron overdoped doped nonsuperconducting 
BaFe 1 . 7 Nio. 3 As 2 at E = 9 ± 3, 19 ±5, 60 ± 10, 96 ± 10, 
181 ± 10 meV ( Wang et al.[|20 13b). The white dashed box 
indicate wave vector integration range at low-energies, while 
the purple dashed boxes in (d-h) mark the integration range 
for high-energy spin excitations. 


ciprocal space [Fig. 20(e)]. For energies above 100 meV, 
spin excitations in hole-doped Bao. 67 Ko. 33 Fe 2 As 2 [Figs. 
20(g)-20(i)] behave similarly to those of electr on-doped 
BaFe 2 - a; Ni a .As 2 (Fig. 19) (Wang et al., 2013b). 

To quantitatively determine the electron and hole¬ 
doping evolution of the spin excitations in iron pnictides, 
one can estimate the energy dependence of the local dy¬ 


namic susceptibility per formula unit \"(w) (Lester et al. 


2010 Liu et al. 2012a). The dashed boxes in Figs. 19 


and 20 show the integration region of the local dynamic 
susceptibility in reciprocal space. At low-energies, we 
only integrate the scattering within the white dashed 
box since it includes all magnetic responses in the Bril- 
louin zone. Approaching to the zone boundary, we inte¬ 
grate the response within the purple dashed box in Fig. 
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FIG. 20 (Color online) Wave vector dependence of spin ex¬ 
citations in hole doped Bai_ a; K I ,Fe 2 As 2 from single crystal 
measurements. Two-dimensional images of spin excitations at 
different energies for hole-doped KFe 2 As 2 at 5 K. (a) E = 8±3 
meV obtained with Et = 20 meV along the c-axis. The right 
side incommensurate peak is obscured by background scat¬ 
tering. (b) 13 ± 3 meV with Ei = 35 meV, and (c) 53 ± 8 
meV with Ei = 80 meV. For Bao.67Ko.33Fe2As2 at T = 45 


K, images of spin excitations at (d) E = 5 ± 1 meV obtained 
with Ei = 20 meV, (e) 15 ± 1 meV with Ei = 35 meV, and 
(f) 50 ± 2 meV obtained with Ei = 80 meV. Spin excitations 
in Ba 0 . 67 Ko. 33 Fe 2 As 2 at energy transfers (g) 115 ± 10 meV; 
(h) 155 ± 10 meV (i) 195 ± 10 meV obtained with Ei = 450 
meV, all at 9 K. Wave vector de pendent background s have 
been subtracted from the images (|Wang et al. (2013b I. 


19 as discussed in Section II A (Fig. 3). The energy 
dependence of the local dynamic susceptibility for hole 
and electron doped iron pnictides are plotted in Figs. 
21(a) and 21(b), respectively. We see that the effect 
of hole-doping near optimal superconductivity is to sup¬ 
press high-energy spin excitations and transfer spectral 
weight to low energies. The intensity changes across T c 
for hole-doped Bao. 67 K 0 . 33 Fe 2 As 2 are much larger than 


that of the electron-doped BaFe1.9Nio.1As2 (Liu et al. 


2012a). As a function of increasing electron-doping, the 


local dynamic susceptibility at low energies decreases 
and finally vanishes for electron-overdoped nonsupercon¬ 


ducting BaFe1.7Nio.3As2 (Luo et al. 2013a Wang et al. 


2013b). 


D. Evolution of spin excitations in iron chalcogenides and 
alkali iron selenides 


Compared with iron pnictides, iron chalcogenide 
(Fei+yTei-xSea;) superconductors have a different static 
AF ordered (bi-collinear instead of collinear) parent com- 


pound (Bao et al. 

2009 Fruclrart et al. 1975] Li et al. 

2009d 

, but a simi 

ar Fermi surface topology (Chen et al. 

2010b 

Nakayama et al. 2010 

Subedi et al. 2008 

). If the 
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FIG. 21 (Color online) Energy and temperature depen¬ 
dence of the local dynamic susceptibility f° r ( a ) 

BaFe 2 As 2 (black solid line), Ba 0 . 67 K 0 . 33 Fe 2 As 2 (purple and 
yellow solid lines for below and above T c , respectively), 
KFe 2 As 2 (blue solid line), and (b) BaFe 2 - :r Ni a; As 2 with x = 
0,0.096,0.15,0.18,0.3, corresponding to black, red, dark red, 
light blue, and green lines respectively. The intensity is in 
absolute units of /^eV _1 f.u. _1 obtained by integrating the 
\"(Q, w) in the dashed regions specified in Figs. 18 and 19 
(Liu et al. 2012a Wang et al. 2013bI. 


resonance originates from the hole and electron Fermi 
surface nesting, one would also expect a resonance at a 
wave vector similar to that of the iron pnictides. The 
neutron scattering experiment s on FeTep. 6 Seo .4 reveal 
that this is indeed the case (Babkevich et al. 2010 


Mook et al. 2010 Qiu et al. 2009). Figure 22(a) shows 
that the resonance energy is weakly temperature depen¬ 


dent and suddenly vanishes above T c (Harriger et al. 


2012 


Qiu et al. 2009). Another interesting aspect of 


Fei+yTei-^Se^ is the presence of transverse incommen- 


and 22(c)[ ( 

Argyriou et al. 

2010 

Lee et al. 

2010 |Li et al. 

2010 Lumsden et al. 

201( 

1). Since the parent compound 


structure, the AF Bragg peaks and associated spin exci¬ 
tations in nonsuperconducting iron chalcogenides stem 
from wave vector positions rotated 45° from those of the 
resonance in reciprocal space. The enhancement of the 
resonance in superconducting Fei+yTei-^Sex occurs at 
the expense of the spin excitations associated with the 


AF nonsuperconducting parent compound (Chi et al. 
2011 Liu et al. 2010 Xu et al. 2010[ ). 


Figure 23 compares the wave vector dependence of spin 
excitations at different energies within the (iL, K) plane 
for nonsuperconducting Fe 1 _ 1 _yTeo. 73 Seo .27 and supercon- 


ducting Fei +y Teo. 5 iSeo .49 (Lumsden et al. 
the tetragonal crystalline lattice unit cel 


2010). Using 
the recipro- 
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FIG. 22 (Color online) Temperature and wave vector depen¬ 
dence of the resonance and low-energy spin excitations in iron 
chalcogenide superconductors, (a) Temperature dependence 
of the resonance energy for optimally doped FeTeo. 6 Seo. 4 . 


The mode energy is essentially temperature independent (Qiu 


et al. 2009). The inset shows the temperature dependence of 


the resonance intensity. The wave vector dependence of the 
spin excitations at different energies along the transverse di¬ 


rection (b) above and (c) below T c for FeTeo. 6 Seo .4 (Argyriou 
et al. |2010| ). 


cal lattice units in Fei+yTei-^Se^ are rotated 45° from 
that for the AF ordered orthorhombic iron pnictides (Fig. 
3). In this notation, spin waves from the bi-collinear 
AF ordered Fei+yTe stem from Qaf = (±0.5, 0) in re- 
ciprocal space while th e resonance occurs at (0.5,0.5) 
(Lumsden et al. 2010). For the nonsuperconducting 
Fe 1 -t-yTeo. 73 Seo. 27 , spin excitations at low energies {E = 
10 ± 1, 22 ± 3 meV) peak at transversely incommensu¬ 
rate positions from (0.5,0.5) [Figs. 23(a) and 23(b)]. 
On increasing the energies to E = 45 ± 5 [Fig. 23(c)] 
and 120 ± 10 meV [Fig. 23(d)], spin excitations be¬ 
come fourfold symmetric and move to positions near 
(±1,0) and (0, ±1) (Lumsden et al. 2010). For super¬ 
conducting Fe 1 - 1 -yTeo. 51 Seo. 49 , the transverse incommen¬ 
surate spin excitations in the nonsuperconducting sample 
at E = 10 ± 1 and 22 ± 3 meV are replaced by the res¬ 
onance and transversely elongated spin excitations near 
(±0.5, ±0.5) [Figs. 23(e) and 23(f)]. Spin excitations at 
E = 45 ± 5 [Fig. 23(g)] and 120 ± 10 meV [Fig. 23(h)] 
are not greatly affected by superconductivity. These re¬ 
sults are similar to spin excitations in electron-doped iron 
pnictides (Liu et al. 2012a), suggesting that supercon- 
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H(r.l.u.) 





120 ±10 meV 


22±3 meV 


FIG. 23 (Color online) Wave vector evolution of the spin 
excitations in FeTei-xSe^ throughout the Brillouin zone. The 
in-plane wave vector dependence of the spin excitations in 
FeTeo. 73 Seo .27 at (a) E = 10 ± 1, (b) 22 ± 3, (c) 45 ± 5, 
(d) 120 ± 10 meV. Identical scans for FeTeo. 51 Seo .49 at (e) 


E = 10± 1, (f) 22±3, (g) 45±5, (h) 120±10 meV (Lumsden 
[eFal) [MTo| . 






Hin(HKO) (r.l.u.) Hin(HKO) (r.I.u.) 


FIG. 24 (Color online) Energy and wave vector dependence 
of the resonance in the superconducting alkali iron selenides. 

(a) Superconductivity-induced neutron scattering intensity 
changes in the superconducting Rb 2 Fe 4 Ses with T c = 32 K. 
A resonance like feature is seen below T c at E = 14 meV. 

(b) The temperature differences in wave vector scans across 
the resonance along the [0.5, I\, 0.5] direction, (c) Possible 
nesting wave vectors connecting the two electron-like Fermi 
surfaces, (d) The temperature difference plot in the in-plane 
reciprocal space reveals the location of the resonance in the 
superconducting alkali iron selenides (|Friemel et al.\ |2012a|b| 


Park et al. 2011b). 


ductivity in iron chalcogenides only affects low-energy 
spin excitations and has commensurate spin excitations 
consistent with the hole and electron Fermi surface nest¬ 


ing (Chi et al. 2011 Liu et al. 2010). 


In iron pnictide and iron chalcogenide superconduc¬ 
tors, the neutron spin resonance is believed to arise from 
quasiparticle excitations between the hole and electron 
Fermi pockets near the T and M points, respectively 
(Hirschfeld et al. |2011 Mazin 20101. Since alkali iron 


selenide superconductors A a ,Fe 2 _ y Se 2 (|Fang et al. 
Guo et al. 2010| do not have hole pockets near the 


2011 


ermi 


energy (Mou et al. 2011 Qian et al. 2011 Zhang et al. 


2011b), it is important to determine if the system also 


has a resonance arising from quasiparticle excitations be¬ 
tween the two electron-like Fermi pockets near the (±1, 0) 
and (0, ±1) positions in reciprocal space (Maier et al., 
2011). From the earlier work on copper oxide super¬ 


conductors, it is generally believed that the resonance 
arises from sign reversed quasiparticle excitations be¬ 


tween two different parts of the Fermi surfaces (Eschrig 


2006). As there are no hole Fermi surfaces in supercon¬ 


ducting 74 a ,Fe2_ySe2, a determination of the location of 
the resonance in reciprocal space will directly test the 
prediction from the RPA and weak coupling calculation 
concerning the nature of the superconducting pairing in¬ 


teraction (Maier et al. 2011). If a resonance is seen ap¬ 
proximately at a wave vector connecting the two electron 


Fermi pockets, one would expect a sign change between 
the two Fermi pockets reminiscent of the d- wave pair¬ 
ing symmetry state of the copper oxide superconductors 


2011a). 


(Das and Balatsky 2011 Maier et al. 2011 Wang et al. 


Experimentally, a neutron spin resonance has been ob¬ 
served at an energy of E r = 14 meV in the superconduct¬ 


ing Rb 2 Fe 4 Se 5 with T c = 32 K [Fig. 24(a)] (Park et al. 
2011b). A complete mapping of the reciprocal space 


within the (H, I\) scattering plane of the system reveals 
that the mode occurs near the wav e vector (0.5,0.25,0.5 ) 
in the tetragonal unit cell notation ( |Friemel et al. 2012b| ). 
Figure 24(a) and 24(b) plots the temperature differ¬ 
ence between 1.5 K(< T c ) and 35 K(> T c ) showing the 
superconductivity-induced resonance in energy and wave 
vector scans, respectively. Figure 24(c) shows the Fermi 
surfaces in the ( H , K , 0) plane corresponding to the dop¬ 
ing level of 0.18 electrons/Fe. The arrows are the in¬ 
plane nesting wave vectors consistent with the resonance 
(Friemel et al. |2012b ). Figure 24(d) plots the difference 
of the RPA calculated dynamic susceptibility between 
the superconducting and normal states at the resonance 
energy (Friemel et al. 2012b). The calculated results 


are in qualitatively agreement with the neutron scatter¬ 
ing experiments, thus suggesting that the mode arises 
from quasiparticle excitations between the electron pock¬ 
ets (Friemel et al. 2012a b). Subsequent neutron scat- 
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tering exper iments on supercon ducting Rbo^FejjjsSe^ 
(T c = 32 K) (|Wang et ed.[|2012b|) and Cs x Fe 2 y Se 2 (Taylor 
et al] 2012 [ ) also found the resonance at wave vector po 
sitions connecting the two electron Fermi surfaces, thus 
confirming this is a general feature of the superconduct¬ 
ing alkali iron selenides. Although the resonance mode 
energy in molecular-intercalated FeSe superconductors 


(Burrard-Lucas et al. 

2013 

Krzton-Maziopa et al. 

2012 

Ying et al. 

2012 

) approximately follows ~ 5 ksT c con- 


et al. 2011 ), its wave vector is better matched to those of 


the superconducting component of A ir Fe 2 _, / Se 2 (Taylor 
et al. 2013). 


E. Impurity effect on spin excitations of iron pnictide and 
chalcogenide superconductors 


As described in the earlier sections, low-energy spin 
excitations in high-T c copper oxide and iron-based su¬ 
perconductors are coupled to superconductivity via the 
opening of a spin gap and re-distributing the weight to 
a neutron spin resonance, both at the AF ordering wave 
vector of their parent compounds (Eschrig 2006). Since 


superconductivity in high-T c superconductors can be al¬ 
tered rather dramatically with impurity doping, it is im¬ 
portant to determine the effect of impurities on spin exci¬ 
tations. In the case of the copper oxide superconductors, 
the resonance and low-energy spin excitations respond 
to magnetic and nonmagnetic impurity doping differently 
(Sidis et al. 2000). When magnetic impurities such as Ni 


are doped into optimally superconducting YBa 2 Cii 207 , 
the resonance peak shifts to lower energy with a pre¬ 
served energy-to-T c ratio (Sidis et al. 2000). In con¬ 


trast, nonmagnetic impurity Zn doping to YBa 2 Cu 207 
restores normal state spin excitations but hardly changes 
the energy of the resonance (Sidis et al. 2000). Similar 


Zn-substitution in the underdoped YBa 2 Cu 2 06.6 induces 
static magnetic order at low temperatures and triggers a 
spectral-weight redistribution from the resonance to the 
low-energy incommensurate spin excitations (Suchaneck 

eToO 12010). 


To see the impurity effect on the resonance and low- 
energy spin excitations in iron-based superconductors, in¬ 
elastic neutron scattering experiments were carried out 
on Ni and Cu-d oped superconducting Fei+yTeo.sSeo.s 
(Xu et al. 20121. Figure 25 shows the temperature de¬ 


pendence of the spin excitations at different energies for 
Fei +!/ _o. 04 Nio. 04 Teo. 5 Seo .5 (Xu et all 2012). In addition 
to reducing the energy of the resonance, the spin exci¬ 
tations at E = 3.5 [Fig. 25(a)], 5 [Fig. 25(b)], and 6.5 
meV [Fig. 25(c)] change from commensurate below T c to 
transversely incommensurate around 100 K. Wave vec¬ 
tor scans at E = 8 [Fig. 25(d)], 11 [Fig. 25(e)], and 
20 meV [Fig. 25(f)] have similar lineshapes on warming 
from 2.8 K to 100 K. Such a dramatic spectral reconstruc¬ 
tion for temperatures up to ~ 3T C is not seen in copper 
oxide and iron pnictide superconductors, and may indi¬ 
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FIG. 25 (Color online) Temperature dependence of the low- 
energy spin excitations in Fei+y-^Ni/Cu^Teo.sSeo.s family 
of iron chalcogenides. Wave vector dependence of the spin ex¬ 
citations along the transverse direction through the AF order¬ 
ing wave vector Qaf for the Ni-doped sample with x = 0.04 
at T — 2.8 K (red circles), 15 K (blue squares), and 100 K 
(green triangles), obtained at (a) E = 3.5, (b) 5, (c) 6.5, (d) 
8, (e) 11, and (f ) 20 meV [which was measured in a higher 
zone, near Q = (1.5, 0.5, 0)]. The low-energy spin excitations 
change from commensurate at low-temperature (T = 2,15 
K) to transversely incommensurate at 100 K. Solid lines are 
guides to the eye (Xu et fli.||20l2) . 


cate the presence of strong electron correlations in iron 
chalcogenide superconductors (Xu et al. 2012). In subse¬ 
quent transport and neutron scattering experiments on 
Cu- doped Feo. 98 - z Cu z Teo. 5 Seo .5 with z = 0,0.02, and 
0.1 (Wen et al. 2013), a metal-insulator transition was 
found for z > 0.02. In addition, low-energy spin exci¬ 
tations of the system are enhanced with increasing Cu- 
doping. These results suggest that the localization of the 
conducting states and electron correlations induced by 


the Cu-doping may play an important role (Wen et al. 


2013). 


While it is well known that hole-doping via K substitu¬ 
tion forBainBaFe 2 As 2 induces high-T c superconductiv¬ 
ity (Rotter et al. 2008a I, substitution of Mn and Cr for 


Fe in BaFe 2 As 2 never induces superconductivity (Sefat 
et al. 2009[ Thaler et al. 2011). In the case of Cr-doping, 


the system adopts a checkerboard G-type AF structure 


for Ba(Fei x Cr x ) 2 As 2 with x > 0.3 [Fig. 6 (a)] (Marty 


et al. 2011). How spin excitations in the parent com¬ 


pound BaFe 2 As 2 are modified by Cr-doping is unclear. 
On the other hand, Mn-doped BaFe 2 As 2 represents a 
more complicated situation: while BaMn 2 As 2 forms a 
simple AF structure with the ordered moment along the 
c axis (Singh et al. 2009), Mn-doping of BaFe 2 As 2 may 
induce a Griffiths regime associated with the suppression 


























































































23 



FIG. 26 (Color online) Spin excitations in Mn impurity 
doped BaFe 2 As 2 and SrCo 2 As 2 pnictides. (a) Spin excitations 
in Ba(Feo. 925 Mno.o 75 ) 2 As 2 with incident beam (Ei = 74.8 
meV) parallel to the crystallographic c axis. Data are dis¬ 
played in the (H + K,H — K) plane and averaged over an en¬ 
ergy transfer of E = 5-15 meV. In addition to the usual spin 
excitations from the collinear AF ordered phase at Qstripe, 
there are spin excitations at QNeei- (b) Spin excitations ema¬ 
nating from QNeei and Qatripe after averaging over the range 
H = 0.50 ± 0.05 (Tucker et al. 2012b). (c) Wave vector de¬ 
pendence of the spin excitations in SrCo 2 As 2 measured with 
incident beam along c axis and Ei = 75 meV at T = 5 K. The 
energy integration range is between E = 15-25 meV, high¬ 
lighting anisotropic spin excitations centered at Qaf- The 
panel (d) plots the same data as in panel (c), but symmetry- 
equivalent quadrants have been averaged together. The wave 


vector anisotropy becomes even more apparent (Jayasekara 


et al., 2013). 


of the collinear AF order in BaFe 2 As 2 by the randomly 
introduced localized Mn moments acting as strong mag¬ 
netic impurities (Inosov et al. 2013). Inelastic neutron 
scattering experiments were carried out on single crys¬ 
tals of Ba(Fei_ a; Mn a ;) 2 As 2 with x = 0.075, which has 
a tetragonal-to-orthorhombic lattice distortion and or¬ 
ders into a collinear AF structure simultaneously below 
T s = Tjy — 80 K (Tucker et al. 2012b). Figure 26(a) 
shows spin excitations of the system measured with the 
crystallographic c axis parallel to the incident neutron 
beam at Ei = 74.8 meV. In addition to spin excitations 
associated with the collinear AF structure denoted as 
Qstripe = Qafj there are excitation s at t h e AF wave vec - 
tor positions of BaMn 2 As 2 (QNeei) (Tucker et al. 2012b). 
At present, it is unclear if this is an intrinsic effect of the 
system or there is real space phase separation between 
Mn and Fe. Figure 26(b) shows energy dependence of the 
scattering at Qstripe and QNeei- While spin excitations 
at Qstripe extend well above 50 meV, they are limited to 
below ~30 meV at QNeei [Fig- 26(b)]. 

In the study of electron-doping evolution of the spin 
excitations in iron pnictides, it was found that electron- 
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FIG. 27 (Color online) Spin excitations in superconduct¬ 
ing and nonsuperconducting LiFeAs without static AF or¬ 
der. (a) Mapping of inelastic neutron scattering intensity 
at E = 5 meV in the ( H, K ) reciprocal space of LiFeAs in 
the tetragonal notation. Note the transverse incommensurate 
peaks away from the Qaf = (0.5, 0.5) position (Qureshi et al. 


2012a). (b) Comparison of incommensurate spin excitations 


for the superconducting (SC) and nonsuperconducting (NSC) 
LiFeAs at E = 5 meV ( |Wang et qZ.||2012a ). (c) Hole and elec¬ 
tron Fermi surfaces in LiFeAs from ARPES (Borisenko et al . I 


2010). The arrow indicates possible nesting condition between 


hole and electron Fermi surfaces. The 8 k v indicates the ex¬ 
pected transverse incommensurability, (d) The experimen¬ 
tally observed transverse incommensurate spin fluctuations 


and its energy dependence for SC and NSC LiFeAs (Qureshi 
| et al . | |2012a| |2014b| |Wang et al. | |2012a| ) . The 5k is the ob 
served transverse incommensurability. 


doping via Co or Ni substitution for Fe in BaFe 2 As 2 in¬ 
duces transversely elongated spin excitations near Qstripe 
due to the mismatched hole and electron Fermi surfaces 
(Figs. 15 and 16) ( |Zhang et al. 2010). If this sce¬ 
nario is correct for all electron-doping levels, one would 
expect transversely elongated spin excitations in heav¬ 
ily Co-doped BaFe 2 As 2 or SrFe 2 As 2 . Figure 26(c) and 
26(d) shows wave vector dependence of spin excitations 
in SrCo 2 As 2 ( Jayasekara et al. 2013). Although spin ex¬ 
citations still appear at the same wave vector positions 
as those of BaFe 2 As 2 , they are longitudinally elongated. 
As SrCo 2 As 2 may have complicated Fermi surfaces like 
that of BaCo 2 As 2 (Xu et al. 2013b), Fermi surface nest¬ 
ing could potentially explain the line-shape of the spin 
excitations. It will be interesting to sort out how spin 
excitations evolve from transversely elongated to longi- 
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(b) m c 



M = sin 2 9-M +cos 2 0-M 

y a c 


FIG. 28 (Color online) Neutron polarization analysis used to 
determine the magnitude of spin excitations along the crys¬ 
tallographic a ( M a ), 6 (Mb), and c ( M c ) axis directions, (a) 
Incident beam neutrons are polarized along the x, y, and z 
directions, corresponding to directions along the momentum 
transfer Q, perpendicular to Q but in the horizontal scatter¬ 
ing plane, and perpendicular to Q and the horizontal scat¬ 
tering plane, respectively. In this geometry, neutron spin-flip 
scattering <rf F ~ M y + M z , where M y and M z are magnitude 
of spin excitations along the y and z directions, respectively. 
Similarly, er jj F ~ M z and uf F ~ M v . If the angle between x di¬ 
rection and the H axis is 9, we have M y = sin 2 9M a + cos 2 9M C 
and M z = Mb- (b) Since M a , Mb, and M c should be the same 
at equivalent wave vectors in reciprocal space except for the 
magnetic form factor, we can conclusively determine M a , Mb, 
and M c by measuring <r® F at two or more equivalent recipro¬ 
cal lattice vectors ( |Luo et al]|2013b| ). 


tudinally elongated in reciprocal space as a function of 
Co-doping for BaFe 2 - x Co x As 2 and SrFe 2 - x Co x As 2 - 
Another iron pnictide worthy of mention is LiFeAs 
(Pitcher et al. 2008[ [Tapp et al. 2008 Wang et al. 2008 1 . 
Although this material has the same crystal structure 
as that of NaFeAs [Fig. 1(d)], it is a superconductor 
without static AF order in stoichiometric LiFeAs and 
Li deficiency tends to suppress superconductivity. Ini¬ 
tial inelastic neutron scattering measurements on pow¬ 
der samples indicate the presence of superconductivity- 
induced resonance near the usual AF ordering wave vec¬ 
tor Qaf (Taylor et al. 2011). Subsequent neutron scat¬ 


tering experiments on single crystals showed that spin 
excitations in this system are transversely incommensu¬ 
rate away from the Qaf for both the superconducting 


LiFeAs [Fig. 27(a)] (Qureshi et al. 

2012a 

and non- 

superconducting Lii_ x FeAs Fig. 27( 

j)] (Wang et al. 


2012a). The absence of the static AF order has been in¬ 


terpreted as due to poor Fermi surface nesting between 
r and M consistent with ARPES measurements [Fig. 
27(c)] (Borisenko et al. 2010[ Brydon et al. 201 1[ > . How¬ 
ever, the quasiparticle scattering between the hole pock¬ 
ets near T and electron pocket M should give rise to 
transverse incommensurate spin fluctuations and this is 


indeed the c ase [Fig . 27(a)] (Qureshi et al. 2012a 2014b 


Wang et al. 2012a). Furthermore, the incommensurate 


spin excitations are weakly energy dependent and only 
broaden slightly for nonsuperconducting Lii_ x FeAs [Fig. 


27(d)] ( |Wang et al . 2012a). These results suggest that 
spin excitations in LiFeAs have the same origin as other 
iron pnictides, and the low-energy spin excitations in the 
system follow the nested Fermi surface scenario. 


F. Neutron polarization analysis of spin excitation 
anisotropy in iron pnictides 


Although most neutron scattering experiments are car¬ 
ried out with unpolarized neutrons, neutron polarization 
analysis can provide some unique information concerning 
the nature of the ordered moment and the anisotropy 
of spin excitations. The neutron polarization analysis 
was first developed at Oak Ridge National Laboratory 
in the pioneering work of Moon, Riste, and Koehler 


(Moon et al. 


1969). This technique was used to un¬ 


ambiguously identify the magnetic nature of the neutron 
spin resonance in optimally doped high-T c copper oxide 
YBa 2 Cu 307 (Mook et al. 1993). In modern polarized 


neutron scattering experiments on high-T c superconduc¬ 
tors at the Institut Laue-Langevin, the Cryopad capa¬ 
bility (Lelievre-Berna et al. 2005) is typically used to 


ensure that the sample is in a strictly zero magnetic field 
environment, thus avoiding errors due to flux inclusion 
and field expulsion in the superconducting phase of the 
sample (Lester et al. 2010). The Cryopad device can 


also be used for neutron polarimetry in which an arbi¬ 
trary incident and scattered neutron beam polarization 
can be measured (Lelievre-Berna et al. 20051. 


Polarized neutrons were produced using a focusing 
Heusler monochromator and analyzed using a focusing 
Heusler analyzer. Polarization analysis can be used to 
separate magnetic (e.g. spin excitations) and nuclear 
(e.g. phonon) scattering because the former has a ten¬ 
dency to flip the spin of the neutron, whereas the latter 
leaves the neutron spin unchanged. More specifically, the 
spin of the neutron is always flipped in a magnetic inter¬ 
action where the neutron polarization is parallel to the 
wave vector transfer Q and the magnetic moment or ex¬ 
citation polarization in the sample is transverse to Q. 
We therefore describe the neutron polarization in a coor¬ 
dinate system where x is parallel to Q. For convenience, 
we then define the other orthogonal directions with y in 
the scattering plane, and z out of plane [see Fig. 28(a)]. 
There are then six independent channels in which the 
instrument can be configured at a specific wave vector 
and energy point: three neutron polarization directions 
x , y, z, each of which can be measured to detect neu¬ 
trons that flip or do not flip their spins when scattering 
at the sample. The measured neutron cross-sections are 
labeled by the experimental configuration in which they 
were measured, and are written <r® F , cr^ SF , where a is 
the neutron polarization direction (x, y or z) and the su¬ 
perscript represents either spin-flip (SF) or non-spin flip 
(NSF) scattering. 

Magnetic neutron scattering only probes the magnetic 
moment perpendicular to Q. The cross-sections can 
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FIG. 29 (Color online) Polarized neutron scattering study 
of spin waves in BaFe 2 As 2 - The AF ordering Brillouin-zone 
centers are at Qaf = (1, 0, L) with L = 1, 3, • • • . The mag¬ 
netic zone boundaries are at Q = (1, 0, L) with L — 0, 2, • • •. 
(a) Inelastic constant-Q scans at the AF zone boundary 
Q = (1,0,0). Here the blue dots are the magnetic scatter¬ 
ing along the y direction, or a y ~ <rf F — cr® F ~ M y , while 
red squares depict <j z ~ crf F — crf F ~ M z . (b,c,d) Similar 
scans at Qaf = (1,0, L) with L = 1,3, and 5, respectively. 
The solid lines are spline-interpolated spin-wave theory cal¬ 


culations folded with the experimental resolution (Qureshi 


et al. 2012b). Using identical experimental setup as that 


of (Qureshi et al. 2012b) but with much more sample mass, 
the energy dependence of M a , Mi,, and M c is determined at 
the magnetic Brillouin (e) zone center and (f) zone boundary. 
The presence of longitudinal spin excitations, or nonzero M a , 
is seen above 20 meV at Qaf = (1,0,1) and above 30 meV 
at Q = (1, 0, 0) ( |Wang et al. 2013a|. 


therefore be written in terms of M y and M z [see Fig. 
28(a)], the two spatial components (perpendicular to Q) 
of the spin direction of the magnetic excitations, and 
the nuclear scattering strength N. However, a mea¬ 
sured cross-section component on an imperfect instru¬ 
ment contains a leakage between SF and NSF channels 
due to imperfect neutron polarization. This leakage can 
be quantified by measuring the nuclear Bragg peak con¬ 
tamination into the spin flip channel, the “instrumen¬ 


tal flipping ratio” R = NSFn/SFn (for an unpolarized 
neutron scattering experiment R= 1, and R —> oo in an 
ideal polarized neutron scattering experiment). The mea¬ 
sured cross-section components can then be written (Lip- 


scombe et al. 

2010 

Moon et al. 

1969) 
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(1) 

where B is a background term to take account of instru¬ 
mental background and the non-magnetic nuclear inco¬ 
herent scattering from the sample, which is assumed to 
be equal in all six cross-sections when measured at the 
same wave vector and energy. NSI is the nuclear spin in¬ 
coherent scattering caused by moments within the nuclei 
of the isotopes in the sample. NSI is independent of Q 
and typically is negligible in magnitude compared with 
the nuclear coherent cross-section N. Furthermore, in 
the case where only SF (or only NSF) cross-section com¬ 
ponents are collected, NSI would be absorbed into the B 
term. 

For SF neutron scattering measurements at Q, one can 
conclusively determine the magnetic components M y and 
M z . If the magnetic components of the system along 
the x, y, and z are M x , M y , and M z , respectively, we 
would have M z = Mb and M y = M a sin 2 9 + M c cos 2 9 
if the sample is aligned in the [ H , 0, L] scattering plane, 
where M a , Mf,, and M c are magnitudes of spin excita¬ 
tions along the orthorhombic a, b, and c-axis directions 
of the lattice, respectively, and 9 is the angle between 
M a and x-axis [Fig. 28(b)] (Lipscombe et al. 2010| ). 
Since there are three unknowns (M a , Mb, M c ) and only 
two equations with known M y and M z , one can only de¬ 
termine the values of M a , Mb, and M c by measuring at 
least two equivalent reciprocal lattice vectors with dif¬ 
ferent 9 angle as illustrated in Fig. 28(b). In the ini¬ 
tial polarized neutron scattering experiments on opti¬ 
mally electron-doped superconducting BaFe 1 .gNio. 1 As 2 , 
spin excitation anisotropy near the resonance energy was 

2010). Similar results were 


observed (Lipscombe et al. 


also found for the resonance in optimal superconduct¬ 
ing Fe(Se,Te) (Babkevich et al.. 2011). For electron- 


overdoped BaFei. 85 Nio.i 5 As 2 , the resonance and spin ex¬ 
citations at all energies probed are isotropic with M y = 
M z (Liu et al. 2012b). 
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FIG. 30 (Color online) Neutron polarization analysis of the 
spin excitation anisotropy in electron-doped iron pnictide su¬ 
perconductors. (a) Energy dependence of the imaginary part 
of the out-of-plane (M c ) and in-plane (Mb) generalized mag¬ 
netic susceptibility at Q = (0.5,0.5,0) in tetragonal nota¬ 
tion for Ba(Feo. 94 Coo.o 6 ) 2 As 2 at T = 2 K. The dashed line 
shows isotropic paramagnetic scattering at T = 30 K. (b) 
Similar data at Qaf = (0.5, 0.5,1), where the mostly c axis 
polarized susceptibility exhibits a peak in the superconduct¬ 
ing phase (Steffens et al .| 2013). (c) Energy dependence of 


M y and M, 


in the normal state at Qaf = (0.5, 0.5, 3) for 
BaFei.904Nio.o96As2 (T c = 19.8 K, Tn ~ T s = 33 ± 2 K). (d) 
Similar data at T = 2 K. The data show clear low-energy spin 
excitation anisotropy in both the normal and superconducting 
states, (e) Temperature dependence of M y and M z , where T* 
marks the temperature of the in-plane resistivity anisotropy, 
(f) Temperature dependence of M a = Muo, Mb = M i_io, 
and M c = Mqoi. Vertical dashed lines mark temperatures for 


T s , Tn, and T c (Luo et al. 2013b). 


Figure 29 summarizes the outcome from polarized neu¬ 
tron scattering experiments on BaFe 2 As 2 - From unpo¬ 
larized neutron scattering experiments, it is well known 
that spin waves in the AF ordered state are gapped be¬ 
low about ~15 meV at the magnetic Brillouin zone center 
Qaf (Matan et al. 2009 Park et al. 2012). Figure 29(a)- 
29(d) shows neutron SF inelastic constant-Q scans at the 
zone center (1,0, L) for different L values and at (1,0,0). 
The blue dots represent the magnitude of magnetic scat¬ 
tering along the y axis direction, or Aly, while red squares 
depict Ms, where cr„ z = (R + l)[crS F — er SF 


My z (Qureshi et al. 2012b). 


yJ/(R - 1 ) 

Consistent with unpolar¬ 


ized measurements fMatan et «L[|2009 Park et al. 2012), 


there are large spin gaps at (1, 0,0) and (1,0,1). How 
ever, the gap value for er? F is significantly larger than 
that for o\y F [Figs. 29(b)-29(d)]. These results indicate 



FIG. 31 (Color online) Polarized neutron diffraction studies 
of the induced magnetization density for different iron pnic¬ 
tide superconductors. The temperature dependence of the 
susceptibility and induced moment for (a) the conventional 
BCS superconductor VsSi, and (b) iron pnictide supercon¬ 
ductor Ba(Feo. 935 Coo.o 65)2 As 2 - The solid lines a re the Yosida 
behavior expected for a singlet order parameter (Le ster et al. \ 
2011). (c) The temperature dependence of the field-induced 
magnetization for superconducting LiFeAs obtained using the 
(1,1,0) and (0,0,2) nuclear Bragg peaks under a 9-T mag¬ 
netic field, (d) The average of the (1,1,0) and (0, 0, 2) shows a 
clear drop below T c , suggesting spin singlet pairing for LiFeAs 
(Brand et al. 1 2014). 


strong single-iron anisotropy within the layer, suggesting 
that it costs more energy to rotate a spin within the or¬ 
thorhombic a-b plane than to rotate it perpendicular to 
the FeAs layers (Qureshi et al. 2012b). In addition, there 


is no evidence for longitudinal spin fluctuations typically 
associated with itinerant electrons and nested Fermi sur¬ 
faces like in the spin-density-wave state of pure chromium 


1994 Qureshi et al. 2012b). 


metal (Fawcett et al. 

In subsequent polarized neutron scattering experi¬ 
ments on BaFe 2 As 2 with greater sample mass (Wang 


et al. 2013a), three distinct spin-excitation components, 
M a , Mb, and M c , with magnetic moments fluctuating 
along the three crystallographic axes are identified at the 
AF Brillouin zone center [Fig. 29(e)] and zone boundary 
[Fig. 29(f)]. The data reveals the presence of finite M a 
at the AF zone center for energies above ~20 meV [Fig. 
29(e)]. Similar measurements at the AF zone bound¬ 
ary suggest nonzero values of M a above ~30 meV [Fig. 
28(f)]. While Mb and M c , the two transverse components 
of spin waves, can be described by a linear spin-wave the¬ 
ory with magnetic anisotropy and inter-layer coupling, 
the presence of M a , the longitudinal component of spin 
waves, is generically incompatible with transverse spin 
waves at low-Q from a local-moment Heisenberg Hamil¬ 
tonian. These results suggest a contribution of itinerant 
electrons to the magnetism that is already in the parent 
compound of this family of Fe-based superconductors. 
This means that one cannot account for spin waves in 
the parents of iron pnictides with a purely local moment 
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picture, and must take the contribution from itinerant 
electrons into account to understand the magnetism in 
these materials (IWang et al. 2013a). 


In a polarized inelastic neutron scattering experiment 
on optimally electron-doped BaFei.8 8 Co 0 .12 {T c = 24 K), 
two resonance-like excitations were found in the super¬ 


conducting state (Steffens et al. 2013). While the high- 


energy mode occurring at E = 8 meV is an isotropic res¬ 
onance with weak dispersion along the c axis, there is a 4 
meV spin excitation that appears only in the c axis polar¬ 
ized channel and whose intensity modulates along the c 
axis similar to spin waves in the undoped BaFe 2 As 2 [Fig. 
30(a) and 30(b)]. These results suggest that spin excita¬ 
tions in undoped and optimally electron doped BaFe 2 As 2 
have similar features, different from what one might ex- 
pect for superco nducting and AF phases of iron pnictides 
(Steffens et al. 2013). 


In a separate polarized inelastic neutron scattering ex¬ 
periment on electron underdoped BaFei.go 4 Nio.o 96 As 2 , 
where the system exhibits an AF order and tetragonal-to- 
orthorhombic lattice distortion temperatures near T/v ss 
T s = 33 ± 2 K, and superconductivity below T c = 19.8 
K, neutron SF cross sections have been measured at var¬ 
ious energies and wave vectors. Figure 30(c) and 30(d) 
shows energy dependence of M y and M z in the normal 
and superconducting states at the AF wave vector (Luo 


et al. 2013b). In addition to confirming that the low- 
energy spin excitations are highly anisotropic below ~5 
meV in the superconducting state [Fig. 30(d)], the mag¬ 
netic scattering appears to be anisotropic in the normal 
state with M z > M y [Fig. 30(c)]. Figure 30(e) shows 
that the magnitudes of M y and M z become different be¬ 
low T*, illustrating that the magnetic anisotropy first 
appears below the temperature where transport mea¬ 
surements on uniaxial strain detwinned samples display 


2010 Fisher 


in-plane resistivity anisotropy (Chu et al. 

2010| ) To quantitatively 


et al. 2011 Tanatar et al. 


determine if the spin excitation anisotropy is indeed 
within the a-b plane, neutron SF cross sections were 
measured at multiple equivalent wave vectors. The out¬ 
come suggests that the presence of in-plane spin excita¬ 
tion anisotropy is associated with resistivity anisotropy 
in strain-induced sample [Fig. 30(f)]. Therefore, spin 
excitation anisotropy in iron pnictides is a direct probe 


ma gnetism may also be coupled to the superconductivity 
(Zhang et al. 2014b). 


Polarized neutron scattering is not only useful for de¬ 
termining the spin excitation anisotropy, it can also be 
used to to measure the susceptibility and induced mag¬ 
netization in the normal and superconducting states of 
a superconductor. The technique of using polarized neu¬ 
tron diffraction to study the magnetization of the para¬ 
magnetic crystal by an externally applied magnetic field 
was developed by Shull and Wedgewood in their study 
of electron sp in pairing of a BCS superconductor VaSi 
(Shull and Wedgewood 1966). Instead of the full neu¬ 


tron polarization analysis as described above, the mag¬ 
netization measurements are performed under a mag¬ 
netic field using a polarized incident beam of neutrons 
(Brown et al. 2010 |Lester et al. 2011). The flipping 
ratio R, defined as the ratio of the neutron scattering 
cross sections with neutrons parallel and antiparallel to 
the applied magnetic field, is associated with the nu¬ 
clear structure factors Fjv(G) and the Fourier transform 
of the real-space magnetization density M(G) via R ss 
1 — 27 roM(G)/[/XB-Piv(G)], where 'yro = 5.36 x 10 -15 


m and G is the reciprocal lattice vector (Lester et al. 


2011). In a conventional BCS superconductor such as 


VaSi, where electrons below T c form singlet Cooper pairs, 
the temperature dependence of the field induced mag¬ 
netization shows the characteristic Yosida drop below 


Wedgewood) 

1966 

. For a spin-triplet superconductor 

such as Sr 2 Ru 04 

Mackenzie and Maeno[ 2003), there 


are no change in the field induced magnetization across 
T c (Duffy et al. 2000). The temperature dependence of 


the field-induced magnetization shows a clear drop below 
T c in nearly optimally electron-doped BaFei . 87 C 00 . 13 AS 2 
[Fig. 31(b)] (Lester et al. 2011), consistent with mea¬ 
surements of the NMR Knight shift in the same com¬ 
pound (Ning et al. 2008[ Oh et al. 2011). The large 
residual contribution of the field induced magnetization 
below T c seen in both VaSi (S hull and We dgewood 1966) 
and BaFei. 87 C 00 . 13 AS 2 (Lester et al. 2011 [ ) has been at¬ 
tributed to the van Vlcck or orbital contribution to the 
susceptibility. 

Similar polarized neutron diffraction experiments have 


also been carried out on superconducting LiFeAs (Brand 


of the spin-orbit coupling in these materials (Luo et al. et al. 20141, which does not have a static AF ordered 


2013b). Recent polarized inelastic neutron scattering ex¬ 


periments on superconducting Ban. 67 Ko 33 Fe 2 As 2 (Zhang 
et al ., 2013b) and Bao. 5 K 0 . 5 Fe 2 As 2 (Qureshi et al. 2014a) 


parent compound 1 

Pitcher et al. 

2008 

Tapp et al. 

2008 

Wang et al. 

2008) and may have triplet electron pair- 


reveal that the low-energy spin excitation anisotropy per¬ 
sists to hole overdoped iron pnictides far away from the 
AF ordered phase. Similar polarized neutron scatter¬ 
ing experiments on underdoped NaFeo.gssCoo.oisAs with 
double resonances [Fig. 17(c)] suggest that the first reso¬ 
nance is highly anisotropic and polarized along the a and 
c axes, while the second mode is isotropic similar to that 
of electron overdoped NaFeo.gasCoo.cnsAs. Since the a 
axis polarized spin excitations of the first resonance ap¬ 
pear below T c , the itinerant electrons contributing to the 


level favoring a ferromagnetic instability (Brydon et al. 


2011). Figure 31(c) shows temperature dependence of the 


field-induced magnetization at wave vectors (1,1, 0) and 
(0, 0, 2). The average of the (1,1, 0) and (0,0, 2) is shown 
in Fig. 31(d). Different from the spin triplet supercon¬ 


ductor Sr 2 Ru 04 (Duffy et al. 2000), the field-induced 
magnetization clearly decreases at the onset of T c , con¬ 


sistent with the spin singlet electron pairing (Brand et al. 


2014). Therefore, the mechanism of superconductivity in 


LiFeAs is likely the same as all other iron-based super- 
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FIG. 32 (Color online) The effect of uniaxial strain on struc¬ 
tural and magnetic phase transitions in the as-grown and an¬ 
nealed BaFe 2 As 2 . (a) Schematic drawing of the device used to 
apply uniaxial strain to detwin single crystals of BaFe 2 As 2 - 
The sample is cut into a square shape with a 0 /b 0 parallel 
to the applied pressure direction, (b) The temperature de¬ 
pendence of the magnetic order parameters under different 
applied unaixial strain for the as-grown BaFe 2 As 2 - The onset 
of the AF ordering temperature increases with increasing uni¬ 
axial pressure (D hital et al.\ 2012]). (c) Magnetic Bragg peak 
intensity at the (1, 0,1) and (0,1,1) positions for the annealed 
BaFe 2 As 2 at zero pressure (green) and P ~ 15 MPa uniax¬ 
ial pressure along the b a axis. No shift in Tjv is seen under 
uniaxial pressure, (d) The blue squares show the tempera¬ 
ture dependence of the (2, —2, 0) nuclear Bragg peak at zero 
pressure. The sharp step at T s is caused by releasing of the 
neutron extinction due to tetragonal to orthorhombic lattice 
distortion. The identical scan under P ~ 15 MPa uniaxial 


pressure is shown as red diamonds (Lu et al. 2014a). 


FIG. 33 (Color online) Temperature dependence of the 
spin excitation anisotropy at wave vectors (1,0,1) and 
(0,1,1) and its comparison with transport measurements for 
BaFe 2 - a; Ni ;r As 2 - (a) Temperature dependence of spin excita¬ 
tions at E = 6 meV for (1,0,1) and (0,1,1) under P ~ 15 
MPa uniaxial pressure. The anisotropy in spin excitations 
vanishes around T = 160 ± 10 K. The inset shows the trans¬ 
port measurement of the in-plane resistivity anisotropy for 
the annealed BaFe 2 As 2 . (b) Temperature dependence of 
E = 6 meV spin excitations at (1,0,1) and (0,1,1) for 
BaFei. 9 i 5 Nio.o 85 As 2 . The data in the top panel were obtained 
by subtracting the background intensity from the peak inten¬ 
sity at every temperature; the data in the bottom panel were 
obtained by fitting the wave vector scans, (c) The electronic 
phase diagram of BaFe 2 - 3: Ni :E As 2 from resistivity anisotropy 
ratio pbo/pao obtained under uniaxial pressure. The spin ex¬ 
citation anisotropy temperatures are marked as T* . The AF 
orthorhombic (Ort), incommensurate AF (IC), paramagnetic 
tetragonal (PM Tet), and superconductivity (SC) phases are 
marked (Lu et al. | 2014a). 


conductors. 


G. Electronic nematic phase and neutron scattering 
experiments under uniaxial strain 

As mentioned in Section F, transport measurements on 
uniaxial strain detwinned electron-doped BaFe 2 As 2 re¬ 
veal clear evidence for the in-plane resistivity anisotropy 
first occurring at a temperature above the zero pressure 
Tjv and T s (Chu et al. 2010 Fisher et al. 2011 Tanatar 


et al. 2010 | ) . As a function of increasing electron-doping, 

the resistivity anisotropy first increases and then vanishes 
near optimal superconductivity (Fisher et al. 2011), con¬ 


sistent with a signature of the spin nematic phase that 
breaks the in-plane fourfold rotational symm etry (C 4 ) of 
the underly i ng tetragonal lattice (Dai et al. 2009[ Fang 
et al. 2008 |Fernandes et al. 2014). NMR experiments 
on 1111 family of materials also indicate the presence of a 
nematic phase below T s (Fu et al. 2012). However, recent 


scanning tunneling microscopy (Allan et al. 2013) and 
transport measurements (Ishid aet al. | 2013[ ) suggest that 
the resistivity anisotropy in Co-doped BaFe 2 As 2 arises 
from Co-impurity scattering and is not an intrinsic prop¬ 
erty of these materials. On the other hand, ARPES mea¬ 


surements on_Co-doped_BaFe 2 As 2 (Yi et al. 2011 ) and 


NaFeAs (Zhang et al. 2012) reveal a splitting in energy 
between two orthogonal bands with dominant d xz and 
d yz character at the temperature of resistivity anisotropy 
in uniaxial strain detwinned samples, thereby suggest¬ 
ing that orbital ordering is also important for the elec¬ 


tronic properties of iron pnictides (Kruger et al. 2009 


Lee et al. 2009 Lv et al. 2009). Finally, since trans¬ 


port measurements were carried out on uniaxial strain 
detwinned samples (Fisher et al. 2011), it is unclear if 


the uniaxial strain can modify the structural and mag¬ 
netic phase transitions in these materials. 

The first neutron scattering experiment carried out un¬ 
der uniaxial strain was on as-grown BaFe 2 As 2 (Dhital 
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et al. 2012). The data show that modest strain fields 


along the in-plane orthorhombic b a axis as shown in Fig. 
32(a) can induce significant changes in the structural and 
magnetic phase behavior simultaneous with the removal 
of structural twinning effects. Both the structural lattice 
distortion and long-range spin ordering occur at temper¬ 
atures far exceeding the strain free phase transition tem¬ 
peratures [Fig. 32(b)], thus suggesting that the resistivity 
anisotropy in transport measurements is a consequence 


of the shift in T)v and T s under uniaxial strain (Dhital 


et al. 2012). In subsequent neutron scattering study of 


the effect of uniaxial pressure on T/v and T s in NaFeAs, 
as-grown and annealed BaFe 2 As 2 (Song et al. 2013a), 
it was found that while the uniaxial strain necessary to 
detwin the sample indeed induces a significant increase 
in T/v and T s for as-grown BaFe 2 As 2 , similar uniaxial 
pressure used to detwin NaFeAs and annealed BaFe 2 As 2 
has a very small effect on their T/v and T s . These re¬ 
sults would suggest that resistivity anisotropy observed 


in transport measurements ( 

Fisher et al. 

2011) 

is an i 

trinsic property of these materials ( 

Song et al., 

2013a) 


In a recent systematic study of magnetic and struc¬ 
tural transitions of the as-grown parent and lightly Co¬ 


doped Ba(Fei_ a ,Co a; ) 2 As 2 under uniaxial pressure (Dhi¬ 


tal et al. 2014), it was found that the uniaxial strain 


induces a thermal shift in the onset of AF order that 
grows as a percentage of T/v as Co-doping is increased 
and the superconducting phase is approached. In addi¬ 
tion, the authors find a decoupling between the onsets 
of the T s and T/v under uniaxial strain for parent and 
lightly-doped Ba(Fei_a,Co^ 2 As 2 _on the f irst order side 
of the tri-critical point (Dhital et al. 2014). 

At around the same time, elastic and inelastic neu¬ 
tron scattering experiments were carried out on the an¬ 
nealed BaFe 2 As 2 (T/v = 138 K), BaFei. 9 i 5 Ni 0 . 085 As 2 
( T c = 16.5 K, T/v = 44 K), and BaFei. 8 sNio.i 2 As 2 
(T c = 18.6 K, tetragonal structure without static AF 
order) to study the temperature dependence of the spin 


excita tion anisotropy at wave vectors (1, 0) and (0,1) (Lu 


et al. 2014a). By comparing the temperature depen¬ 


dence of the magnetic order parameters at wave vectors 
(1, 0,1) and (0,1,1) in zero and ~15 MPa uniaxial pres¬ 
sure on annealed BaFe 2 As 2 [Fig. 32(c)], it was concluded 
that the applied uniaxial strain is sufficient to completely 
detwin the sample, and does not affect T/v- Further¬ 
more, the temperature dependence of the intensity at the 
(2, —2,0) nuclear Bragg reflection for the twinned and de- 
twinned samples both show a dramatic jump at T s = 138 
K arising from the neutron extinction release that occurs 
due to strain and domain formation related to the or¬ 
thorhombic lattice distortion, indicating that the uniaxial 
pressure does not change the tetragonal-to-orthorhombic 
structural transition temperature [Fig. 32(d)]. Since the 
measurable extinction release at temperatures well above 
T s was suggested to arise from significant structural fluc¬ 


tuations related to the orthorhombic distortion (Kreyssig 
et al. 2010), data for the detwinned sample indicates that 


the applied uniaxial pressure pushes structural fluctua- 



FIG. 34 (Color online) The electronic phase diagram of 
Ba(Fei_ 3 ,Co a; ) 2 As 2 as determined from combined /rSR and 
neutron diffraction experiments, (a) Phase diagram of Tjv, as 
determined with /rSR and T c , obtained from resistivity and 
magnetic susceptibility measurements as well as from the spe¬ 
cific heat data, (b) The Co dependence of the normalized T/v 
and T c , and the normalized values of the average magnetic 
field at the muon site, B M and of its relative spread, A B^. 
The open symbols show the magnetic properties in the spa¬ 
tially inhomogeneous magnetic state near optimum doping 
(Bernhard et al 20121. 


tions to a temperature similar to that at which resistiv¬ 
ity anisotropy emerges. These results are different from 
those of Ref. (Dhital et al. 2014) carried out on as- 


grown BaFe 2 As 2 - It remains to be seen how the different 
results in these experiments can be reconciled (Dhital 


et al. 

2014 

2012 

Lu et al. 

2014a 

Song et al. 

2013a 


In addition to determining the effect of uniaxial strain 
on structural and magnetic phase transitions in annealed 
BaFe 2 As 2 , inelastic neutron scattering experiments on 
BaFe 2 - a; Ni a ;As 2 also reveal that low-energy spin excita¬ 
tions in these materials change from fourfold symmet¬ 
ric to twofold symmetric in the uniaxial-strained tetrag¬ 
onal phase at temperatures corresp onding to the on¬ 
set of in-plane resistivity anisotropy (Lu et al. 2014a). 


The inset in Figure 33(a) shows the in-plane resistiv¬ 
ity anisotropy on annealed BaFe 2 As 2 under uniaxial 
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strain. The temperature dependence of the E = 6 
meV spin excitations (signal above background scatter¬ 
ing) at (1,0,1) and (0,1,1) is shown in Fig. 33(a). In 
the AF ordered state, there are only spin waves at the 
AF wave vector Qaf = (1,0,1). On warming to the 
paramagnetic tetragonal state above TV and T s , we see 
clear differences between (1,0,1) and (0,1,1) that van¬ 
ish above ~160 K, the same temperature below which 
anisotropy is observed in the in-plane resistivity [inset 
in Fig. 33(a)]. Similar measurements on underdoped 
BaFei. 9 i 5 Ni 0 . 085 As 2 reveal that the E = 6 meV spin 
excitations at the (1,0,1) and (0,1,1) wave vectors be¬ 
comes anisotropic below ~80 K [Fig. 33(b)], again consis¬ 
tent with the in-plane resistivity anisotropy from uniax¬ 


ial strain detwinned BaFe 2 _ x Ni E As 2 [Fig. 33(c)] (Fisher 


et al. 2011). Finally, uniaxial strain on electron over¬ 


doped BaFei. 88 Nio.i 2 As 2 induces neither spin excitations 
nor in-plane resistivity anisotropy at all temperatures 
(Fisher et al. 2011 Lu et al. 2014a). Therefore, resis¬ 
tivity and spin excitation anisotropies both vanish near 
optimal superconductivity and are likely intimately con¬ 
nected, consistent with spin nematic phase induced elec¬ 
tronic anisotropy (Fernandes et al. 2014). 


H. Comparison of magnetic order and spin excitations 
determined from neutron scattering with those from fiSR, 
NMR, and RIXS measurements 


In the discussions above, we have summarized recent 
progress in the AF order and spin dynamics in iron-based 
high temperature superconductors determined from elas¬ 
tic and inelastic neutron scattering. However, neutron 
scattering is not the only technique to study magnetism 
in these materials. Other probes such as /iSR, NMR, 
and RIXS experiments have also been used to investi¬ 
gate their magnetic properties. Compared with neutron 
scattering, which is a global probe determining the av¬ 
erage magnetic properties of the bulk solids, the /iSR 
technique is a local probe measuring the magnitude of 
the static random field experienced by muons that re¬ 
side on interstitial lattice sites of the studied material. 
In addition to being able to detect static AF order and 
superfluid density, it can determine the volume fractions 
of the magnetic and supercon duct ing phases and their 
temperature dependence (Uemura 2009; ?). Similarly, 
NMR is also a local probe that can detect magnetic and 
sup ercon d ucting properties of the studied materials (Al- 
loul et al., [2009). In this section, we briefly summarize 


recent results from these techniques and compare them 
with results obtained from neutron scattering. 

We begin by comparing the electronic phase diagrams 



the /xSR results find full volume AF ordered phase co- 



FIG. 35 (Color online) NMR determination of the phase 
diagram for electron-doped iron pnictides. (a) The AF or¬ 
der TV and T c in Ba(Fei_a,Coj;) 2 As 2 as determined from the 
NMR measurements. The left and right insets are schematic 
representations of the Fermi surfaces in unfolded first Bril- 
louin zone for the x < 0.15 and overdoped nonsupercon¬ 
ducting samples, respectively. The absence of electron-hole 
Fermi surface nesting is believed to be responsible for the 
suppression of superconductivity in the overdoped regime, 
(b) The Co-doping dependence of the strength of the para¬ 
magnetic spin excitations as measured by 1/T\T at 25 K 
(> T c ). (c) Weiss temperature 6 obtained from fitting the 
interband (electron-hole pocket excitations) AF spin excita¬ 
tions (1/ Ti T) j n ter with a Curie-Weiss term (l/TiT)j n ter = 
C/{T + 9). Here 1/TiT = (1/TiT) inter + (l/TiT) intra , where 
the intraband scattering can be fitted with a phenomenolog¬ 


ical form (l/TiT)i n tra = a + /3exp(—A/fesT) (Ning et al. 
[2010). (d) The electronic phase diagram of BaFe 2 _ x Ni a ,As 2 


as determined from NMR and transport measurements (Zhou 
et al[ |2013b). 


existing and competing with superconductivity. This 
is consistent with neutron diffraction results indicating 
static commensurate AF order coexists and competes 
with su perconduct ivity in the underdoped regime (Fig. 
13) (Christianson et al. 2009 Pratt et al. 2009). For 
Ba(Fei_ a ,Co s ) 2 As 2 near optimal superconductivity with 
x > 0.05, magnetic order only develops in parts of the 
sample volume and the normalized values of the average 
magnetic field at the muon site, are reduced dramat¬ 
ically [Fig. 34(b)]. In addition, the AF Neel temperature 
of the system saturates to a value near or slightly above 
T c before vanishing in a first order fashion with increas¬ 
ing Co-doping [Fig. 34(b)] (Bernhard et al. 2012). This 
region of Co-doping is consistent with the appearance of 
the transverse incommensurate AF order seen in neutron 
diffraction experiments (Pratt et al. 2011). However, in¬ 
stead of a uniform incommensurate spin-density-wave or¬ 
dered phase, the /iSR data indicates a spatially inhomo¬ 
geneous magnetic state for which the volume fraction of 
the ordered phas e decr eases with increasing Co-doping 
(Bernhard et al. 2012). The NMR measurements on 
Ba(Fei_ a; Co z ,) 2 As 2 samples near optimal superconduc¬ 
tivity reveal that the system is in the spin-glass state 
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FIG. 36 (Color online) The elastic properties of the iron- 
based superconductor Ba(Fei_ a .Co :E ) 2 As 2 as determined from 
ultrasonic spectroscopy measurements. The temperature 
dependence of the anisotropic elastic stiffness Cee for (a) 
BaFe 2 As 2 and (b) BaFei. 84 Coo.i 6 As 2 as determined from res¬ 
onant ultrasonic spectroscopy. The solid lines are fits from 
a model considering spin nematic phase ((Fernandes et al. 
2010). (c) Temperature dependence of the inverse of the 

elastic stiffness 566 = 1/C66 = S° 6 + 566, Cr, where Sg 6 is 


the normal (background) contribution, for Ba(Fei_ :E Co :E ) 2 As 2 
with different x. (d) Structural and magnetic phase diagram 
of Ba(Fei_ a; Co :! ;) 2 As 2 as determined from ultrasonic spec¬ 
troscopy measurements. A structural quantum critical point 


is identified near optimal superconductivity (Yoshizawa et al. 


2012 ). 


which competes with superconductivity (Dioguardi et al. 


2013). Instead of being a consequence of Fermi surface 
nesting ( |Pratt et al. 2011), the incommensurate mag¬ 
netic order in Co-doped BaFe2As2 may arise from inho¬ 
mogeneous short-range magnetic order similar to the mo¬ 
ment modulating cluster spin glass in Ni-doped BaFe2As2 
(Lu et al. [2014b). Therefore, the electronic phase dia¬ 
grams of the Co- and Ni-doped BaFe2As2 are similar, 
showing a homogeneous commensurate long-range static 
AF ordered phase competing with superconductivity in 
the underdoped region and a short-range spin gl ass phase 
near optimal superconductivity [Fig. 5(b)] (Bernhard 


et al. , 20T2j|Dl oguard i et al. l |20T3l|WTZ||2013[|2014~b[ 


Luo et al. 2012b). There is no evidence for a conven¬ 


tional magnetic quantum critical point near optimal su¬ 
perconductivity, and the AF order disappears in the su- 


first order fashion ( 

Bernhard et al. 

2012 

Dioguardi et al. 

2013 

Lu et al. 

2013 

2014b 

Luo et al. 

2012 b 

)• 


The electron phase diagrams of the Co- and Ni- 
doped BaFe 2 As 2 systems have also been mapped out 
by NMR m e asurements [Figs. 35(a) and 35(d)] ( |Ning 
et al.. 2009 2010 Zhou et al. 2013bI. In the case of 
Ba(Fei_ a: Co a ;)2As2, the NMR experiments on electron 
overdoped samples suggest that the absence of quasi¬ 
particle excitations with momentum transfer Qaf be¬ 
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FIG. 37 (Color online) RIXS measurements of the high- 
energy spin excitations in optimally hole-doped BaFe 2 As 2 . 
(a) Schematic view of the reciprocal space covered by Fe L 3 
RIXS is shaded by a yellow circle. T, B and C are the recip¬ 
rocal space positions at which RIXS spectra were collected. 
Black (blue) squares represent the tetragonal (orthorhombic) 
Brillouin zone. All RIXS spectra use the orthorhombic Bril- 
louin zone convention for defining relative momentum trans¬ 
fer values, r point is the structural zone center, while F-M is 
the AF ordering wave vector, (b) Three typical RIXS spec¬ 
tra of BFe 2 As 2 collected at 15 K with n polarized incoming 
light, at F, B, and C positions in reciprocal space, (c) Disper¬ 
sion of spin excitations of BFe 2 As 2 (BFA) in the AF phase, 
and Bao. 6 Ko. 4 Fe 2 As 2 (BKFA) in the superconducting phase, 
(d) Half Width Half Maximum (HWHM, damping) and inte¬ 
grated intensity of spin excitations of BFA and BKFA. The 
horizontal dotted line in upper panel marks the HWHM of 
the total instrumental resolution of the RIXS experiment (40 
meV) (Zhou et al. 2013aI. 


tween the hole and electron Fermi surfaces results in 
complete suppression of the low-energy spin fluctuations 
for x > 0.15. The insets in Fig. 35(a) show that the 
hole bands sink below the Fermi surface for x > 0.15, 
disallowing inter-band quasiparticle transitions between 
the hole and electron Fermi surfaces (Ning et al. 2010). 
These results are consistent with neutron scattering data 
on electron-doped BaFe2_ a ,Ni ir As2, where there is a large 
(~50 meV) s pin gap in th e nonsu perconducting sample 
with x = 0.3 (Wang et al. 20131?). 


In addition to determining the electronic phase dia¬ 
gram of Ba(Fei_ 2; Co a; )2As2 as shown in Fig. 35(a), the 
1/T]T obtained by NMR measurements is also related 
to the wave vector integral of the low-energy spin dy¬ 
namic susceptibility x"(Q, /) via 1/TiT = A/(T — 9) ~ 
|A(Q)| 2 x"(Q, /)//, where 9 is the Curie-Weiss tem¬ 
perature (the temperature at which a plot of the re¬ 
ciprocal molar magnetic susceptibility against the ab- 
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solute temperature T intersects the T-axis), / is the interpreted as due to the spin nematic phase (Fernan- 


NMR frequency, a is the lattice constant, |A(Q)| = des et aL [2010). Figure 36(c)_shows the temperature 

\Acos(Q x a/2) cos(Q y a/2)\ 2 is the form factor of trans¬ 


ferred hyperfme coupling at the 75 As sites, and the wave 
vector summation of Q is taken over the entire first Bril- 
louin zone (Ning et al. 2009). By measuring the Co¬ 


doping dependence of 1/TiT, one can fit the data with 
Curie-Weiss law and obtain the electron doping depen¬ 
dence of 6. Figure 35(b) and 35(c) shows the Co-doping 
dependence of the 1 /T{T in the normal state (T > T c ) 
and 0(K), respectively ( |Ning et al. \ 2010). The enhance¬ 
ment of 1/TiT and the negative to positive crossing of 
9(K) near vanishing AF order suggest the presence of 
a magnetic quantum critical point at x ~ 0.07 (Ning 


et al. 2010). Similar NMR data and transport mea¬ 


surements on BaFe2-£cNi K As2 suggest the presence of 
two quantum critical points associated with AF order 
and the tetragonal-to-orthorhombic lattice distortions, 
respectively [Fig. 35(d)] (Zhou et al. 2013b). These 
results are in direct contrast to the conclusions of neu¬ 
tron scattering and fi SR experiments as described earlier 
in the session, which suggests a weakly first order phase 
transition from static AF order to superconductivity in 
electron-doped BaFe 2 As 2 . We note that these NMR mea¬ 
surements have not been carried out for the Co and Ni- 
doped samples with incommensurate magnetic order and 


near optimal superconductivity (Ning et al. 
et al. 2013b). 


2010 Zhou 


Although determining whether a conventional mag¬ 
netic quantum critical point exists in the iron pnic- 
tides is important, it is equally important to under¬ 
stand what happens to the tetragonal-to-orthorhombic 
lattice distortion in the Co and Ni-doped BaFe 2 As2 phase 
diagram near optimal superconductivity. From initial 
high-resolution X-ray diffraction experiments, it was sug¬ 
gested that Ba(Fei_ a; Co a ;)2As2 has reentrant behavior 
near optimal superconductivity, exhibiting a tetragonal- 
orthorhombic transition above T r and orthorhombic- 


tetragonal structural transition below T c (Nandi et al. 
2010). While the overall trends of these results were 


confirmed by later high-resolution X-ray diffraction ex¬ 
periments on BaFe2- a; Ni a ,As2 near optimal superconduc¬ 
tivity, the presence of low-temperature incommensurate 
AF order suggests that the system is still in the or¬ 
thorhombic phase (Lu et al. 2013). Using resonant 


ultrasonic spectroscopy, one can measure the tempera¬ 
ture dependence of the anisotropic elastic stiffness C\j 
associated with the tetragonal-orthorhombic lattice dis¬ 


tortion in Ba(Fei_ a ,Co a; )2As2 (Fernandes et al. 2010 


Yoshizawa et al.\ 2012). Figure 36(a) and 36(b) shows 


the temperature dependence of the resonant ultrasonic 
spectroscopy measured squared resonant frequency / 2 
(red points), and of the calculated elastic stiffness (shear 
modulus) C s = Cqq of the tetragonal phase (solid lines) 


for BaFe 2 As2 and BaFei.84Coo.i6As2, respectively (Fer¬ 


nandes et al. 2010). The dramatic softening of the Cqq 


shear modulus at temperatures above the tetragonal- 
orthorhombic lattice distortion temperature T s has been 


dependence of the inverse of C 6e for Co-doping levels 
of x = 0,0.037,0.060,0.0 84,0.098,0.116,0.161, and 0.245 
(Yoshizawa et al. 2012). The resulting phase diagram 
shown in Fig. 36(d) suggests the presence of a structural 
quantum critical point near x « 0.07. While this is con¬ 
sistent with the NMR results of a magnetic quantum crit¬ 
ical point (Ning et al. 2009 20101, it differs from the first 


order nature of the AF order to superconducting phase 
transition in electron-doped pnictides determined from 


neutron and /xSR experiments 

Bernhard et al. 

Dioguardi et al. 

2013 

Lu et al. 

[2013 

2014b 

)• 


2012 


In the past, the only way to explore the energy and 
wave vector dependence of the spin excitations in mate¬ 
rials was via inelastic neutron scattering, which in prin¬ 
ciple can map out the spin excitations in absolute units 
from low energy to high energy throughout the Brillouin 
zone as described in section III of this review. However, 
such technique suffers from the need for large amounts 
of single crystals, which may not be available. Recent 
advances in RIXS provide an alternative method to look 


for high-energy spin excitations in copper (Tacon et al. 


2011) and iron (Zhou et al. 2013aI based high-T c super¬ 


conductors, although the precise RIXS cross section is 
difficult to calculate and includes couplings to orbital and 
electronic excitations in addition to magnetic excitations 
(Ament et al. 2011). For the hole-doped YBa 2 Cu306+a; 


family of materials, RIXS measurements showed the ex¬ 
istence of damped but well-defined dispersive magnetic 
excitations, deep in the Stoner continuum of hole-doped 
cuprates with doping beyond the optimal level (Tacon 


et al. 2011). The high-energy spin excitation spectral 


weights are found to be similar to those of spin waves in 
the undoped, anti ferromagnetically ordered parent mate¬ 
rial (Tacon et al. 2011). So far, these measurements on 
hole-doped copper oxide superconductors have not been 
independently confirmed by inelastic neutron scattering 
experiments, which mostly probe spin excitations near 
the AF ordering wave vector instead of near the origin as 
in RIXS experiments due to the small neutron scattering 
cross section near T (Fig. 4) (Fujita et al. 2012 Tran- 
quada et al. 2014). Given that the RIXS cross section 
is not well-known (Ament et al.. 20111, it is extremely 


important to compare inelastic neutron scattering and 
RIXS experiments on similar samples. Recently, such a 
comparison has been made for high-energy spin excita¬ 
tions of Sr 2 Cu02Cl2 (Plumb et al.\ 2014). 


Figure 37(a) shows schematic view of the reciprocal 
space that can be covered by Fe L 3 RIXS shaded by a yel¬ 
low circle. Black (blue) squares represent the tetragonal 
(orthorhombic) Brillouin zone (Zhou et al. 2013a). Sim¬ 


ilar to the case of copper oxide superconductors, RIXS 
by Fe L 3 edge will not be able to cover the same region 
of the reciprocal space as that of inelastic neutron scat¬ 
tering [Fig. 3(e)]. Figure 37(b) plots Fe L 3 edge RIXS 
spectrum on BaFe 2 As2 collected at 15 K with n polar¬ 
ized incoming light, at ( Q Xl Q y ) = (0,0) (T), (0.5,0) 
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FIG. 38 (Color online) Local moments of various iron- 
based superconductors in their paramagnetic states deter¬ 
mined by X-ray emission spectroscopy and core level pho¬ 
toelectron spectroscopy, (a) The IAD values derived from 
the XES spectra for various samples. The room-temperature 
data are shown in circles, and the low-temperature IAD 
values at T = 15 K are shown in triangles for Fei.i 2 Te 
and Rb 2 Fe 4 Se 5 ( |Gretarsson et al\ |2011| ). The right hand 
axis is scaled to the total (static+fluctuating) moments of 
Rbo. 89 Fei. 58 Se 2 from ( |Wang et al. 2011c I (Y.-J. Kim, pri¬ 
vate communication), (b) Temperature dependence of the 
IAD values for Cao. 9 2 Ndo.Q 8 Fe 2 As 2 , Cao. 7 8 L ao. 22 Fe 2 As 2 , and 
Cao. 85 Ndo.i 5 Fe 2 As 2 (Grctarsson et al. 2013). (c) Estimation 
of the spin moment on the Fe sites from the multiplet en¬ 
ergy separation AF 3( g obtained from PES measurement. The 
continuous line is the extrapolation of the linear fit of the 
AE 3 S values plotted against (2 Sv + 1) for the Fe ionic com¬ 
pounds FeF 3 , FeF 2 , FeO, for which Sv is known to be 5/2 
(FeF 3 ) and 2 (FeF 2 , FeO). Here SFA, SFA-10%, CFAO, and 
CFAO-11% are SrFe 2 As 2 , SrFei. 76 Coo. 24 As 2 , CeFeAsO, and 
CeFeAsOo.89Feo.11, respectively. The linear fit results in the 
relation AF 35 = 0.94+1.01(2SV-|-l) (Vilmercati et al. 20121. 


(. B ) and (0.35, 0.35) (C) using orthorhombic recipro¬ 
cal lattice units [see Fig. 37(a) for the T, B , and C 
positions in reciprocal space]. In addition to the in¬ 
tense Fe 3d fluorescence at around —2 eV energy trans¬ 
fer, there are momentum dispersive excitations centered 
around 200 meV near the quasi-elastic peak at zero en¬ 
ergy [Fig. 37(b)]. Since the dispersions of these excita¬ 
tions are identical to those of spin waves in BaFe2As2 
determined from inelas tic neutron scattering [Fig. 37(c)] 
(Harriger et al. 2011), they are believed to arise from 
the same spin waves in BaFe2As2 measured by RIXS 


(Zhou et al. 2013a). Similar RIXS measurements on 
optimally hole-doped superconducting Bao.6Ko.4Fe2As2 
(BKFA) indicate that while the dispersions of spin ex¬ 
citations are softened [Fig. 37(c)], the excitation inten¬ 
sity and widths are unchanged from its undoped parent 
BaFe2As2 [Fig. 37(d)] (Zhou et al.\ |201 3a). This is a 
surprising result, different from the hole-doping evolu¬ 
tion of the high-energy spin excitations in copper oxides 
(Tacon et al. 2011). While the dispersion of the spin 
excitations in BKFA determined by RIXS has been con¬ 
firmed by later inelastic neutron scattering experiments, 
the absolute spin excitation intensity in BKFA measured 
by neutrons is much smaller than that of the spin waves 
in BaFe2As2 (Wang et al. [2013b ). At present, it is un¬ 
clear how to reconcile the RIXS and neutron scattering 
results, although we note that region of the reciprocal 
space probed by RIXS is different from that probed by 
inelastic neutron scattering. 

As discussed in the introduction, neutron scattering 
can measure the overall strength of the magnetic excita¬ 
tions through the local fluctuating moment (m 2 ). For 
a local moment Heisenberg model with spin S , the inte¬ 
grated spectral weight when integrated over all wave vec¬ 
tor and energy space should yield (m 2 ) = ( gfj.B) 2 S(S + 
1). Therefore, one can in principle determine the mag¬ 
nitude of S by measuring elastic and inelastic magnetic 
scattering of the system in absolute units throughout the 
Brillouin zone, and the outcome should be the same as 
the local moment S determined from the X-ray emission 


spectroscopy (XES) (Rueff and Shukla 

2010) and core 

level photo-emission spectra (PES) (Vi 

mercati et al. 


2012). However, the current unpolarized time-of-flight 


neutron scattering technology can only measure corre¬ 
lated spin excitations, and thus will underestimate the 
size of the effective spin S when spin excitations become 
diffusive and broad as in the case of most iron pnictides 
[except for th e true local moment system Rbo.89Fei.5sSe2 
(Wang et al. 2011c)]. By measuring the local dynamic 
susceptibility for electron-doped BaFe2- a; Ni :E As2, we see 
that the magnitude of (m 2 ) decreases from (m 2 ) ss 3.5 
for x = 0 to 2.7 ± 0.1 /z^/Fe, correspo nding to S w 1/2 
(Luo et al. 2013a Wang et~al\\20 13b ). 

Figure 38(a) shows the local moment S for various 
iron-based superconductors in the paramagnetic phase 
as determined from XES (Gretarsson et al. 20111. Here 
the material variation of the local moment of Fe can be 
extracted using the overall shape of the Fe K(3 emis¬ 
sion spectra by applying the integrated absolute differ¬ 
ence (IAD) analysis. The absolute values of S were ob¬ 
tained by scaling the results to neutron scattering data. 
Figure 38(b) plots the temperature dependence of the 
IAD for Cao.92Ndo. o8Fe2As 2 , Cao.7 8 Lao.22F e2As2, and 
Cao.85Ndo.i5Fe 2 As2 (Gretarsson et al. 2013). The local 
moments in the Nd- and Pr-doped samples disappear in 
the cT phase (Goldman et al. 2009), indicating that the 
Fe 2+ ions go through a spin-state transition by taking 
on the low-spin state in the cT phase. Inelastic neutron 
scattering experiments on CaFe2As2 also reveal vanish- 
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FIG. 39 (Color online) The dispersions of spin excitations in 
various copper oxide superconductors and the energy scale of 
the spin excitations, (a) Magnetic dispersion relation along 
the (0.5 + h, 0.5, 0)tetra in various cuprates, corresponding to 
wave vectors parallel to the Cu-O bonds. The energy is scaled 
by J ~ 110 meV for the AF parent compound La 2 CuC> 4 . (b) 
Large symbols: estimates of doping-dependent energy scale at 
which magnetic spectral weight falls below half that of the AF 
state based on inelastic neutron scattering studies of various 
cuprates (?). Small gray symbols: pse udogap en ergy f rom 
various electronic spectroscopies (Tranquada et al. 20141. 


ing spin excitations in the cT phase (Soh et al. 


similar to the XES results (Gretarsson et al. 20131. Fig- 


2013), 


ure 38(c) shows the values of local Fe moment for var¬ 
ious iron pnictides estimated from PES measurements 
(Vilmercati et al. 2012). We see that the values of ~ 2 
Hb for BaFe 2 As 2 and SrFe 2 As 2 family of iron pnictides 
are much larger than those obtained by inelastic neutron 
scattering, and decrease by about 40% in the optimally 
doped compound, also different fr om doping depen dence 
results from neutron scattering (Mannella 2014). As 
emphasized in a recent Review (Mannella 2014), the 
PES experiments sample spin excitations on time scales 
shorter than 10 _15 -10 -16 s, while inelastic neutron scat¬ 
tering measurements probing spin excitations up to ^300 
meV correspond to time scales on the order of 10-15 fs, at 
least one or two orders of magnitude slower than those 
in PES experiments. Furthermore, PES measurements 
determine the total uncorrelated spin while present in¬ 
elastic neutron scattering measures correlated spin exci¬ 
tations near the AF ordering wave vector within the first 
Brillouin zone. 


I. Comparison of spin excitations in iron-based 
superconductors with those in copper oxide and heavy 
fermion superconductors 


With the discovery of I^-xBa^CuC^ family of cop¬ 


per oxide superconductors in 1986 (Bednorz and Muller 


1986), the field of high-T c superconductivity was born. 


Although research in high-T c copper oxide superconduc¬ 
tors is still active with many exciting new results almost 


30 years later (Fradkin and Kivelson 2012 Ghiringhelli 


et al. 20126, the discovery of iron-based superconduc¬ 


tors in 2008 provided an entire new class of materials 
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FIG. 40 (Color online) Neutron spin resonance in heavy 
fermion superconductor CeCoIns and its magnetic field de¬ 
pendence. (a) The imaginary part of the dynamic suscepti¬ 
bility at Qaf = (1/2,1/2,1/2) in the normal (3 K) and super¬ 
conducting (1.3 K) states. A background scattering taken at 
Q = (0.3, 0.3, 0.5) and (0.7, 0.7, 0.5) was subtracted. The hor¬ 
izontal bar is the instrumental resolution width (Stock et al. 


2008). (b) The peak positions of the field spit resonance as 


a function of applied field within the a-b plane. The solid 
lines are fits to E± = Tiu>o ± g^Bf-oH/2 with g = 1.92 ± 0.10. 
(c) Resonance at zero applied magnetic field, (d) Identical 
scan under an applied field of 3-T. A 10 K background was 
subtracted from the scans (Stock et al. 2012). 


where high-T c superconductivity occurs (Kamihara et al. 


2008). Since high-T/ copper oxide and iron-based super¬ 


conductors, as well as heavy fermion superconductors are 
close to AF instability, magnetism may be a common 


thread for unconventional superconductivity (Scalapino 


2012). If this is indeed the case, it will be interesting to 


determine the similarities and differences in spin excita¬ 
tions of these superconductors. 


As discussed in recent review articles (Fujita et al. 


2012 Tranquada et al. 2014), an important feature of 


the spin excitations in copper oxide superconductors is 
the neutron spin resonance and hourglass like dispersion 
of the spin excitation spectra. An hourglass magnetic 
spectrum is also seen in an insulati ng, hole-dop ed an- 
tiferromagnet La 5 / 3 Sr 1 / 3 Co 04 (|Boothroyd et al. 20111. 
Figure 39(a) shows the spin excitation dispersions of dif- 
ferent copper oxide superconductors away from the AF 
ordering wave vector Qaf = (0.5, 0.5), where J is the 
nearest neighbor magnetic exchange coupling (J ss 120 
meV) ( |Fujita et al. 2012 Tranquada et al. 2014). While 
the hourglass dispersion of spin excitations appears to be 
a ubiquitous feature of different families of hole-doped 
copper oxide superconductors, they are clearly absent in 
spin excitation spectra of electron-doped iron pnictides. 
For optimally hole-doped iron pnictides, spin excitations 
change from longitudinally to transversely elongated el¬ 
lipses centered at Qaf = (1, 0) on moving from below the 
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resonance energy to above it (Wang et al. 2013b). This is 
different from the hourglass dispersion seen in hole-doped 
copper oxides. Another important feature of the spin ex¬ 
citation spectra is the electron and hole-doping evolution 
of the local dynamic susceptibility. For hole-doped cop¬ 
per oxide superconductors, the strength of the high en¬ 
ergy magnetic response near Qaf = (0.5, 0.5) decreases 
with increasing doping level. To quantify this behavior, 
Stock et al. evaluated the hole doping dependence of the 
energy at which the local dynamic susceptibility 
falls below half of that for an undoped AF parent com¬ 
pound (Stock et al. 2010). The outcome suggests that 
the hole-dependence of the magnetic energy scale corre¬ 
sponds very well with the pseudogap energy determined 


f rom electronic spectroscopies [Fig. 39(b)] (Hiifner et al. 


2008). These results indicate that AF spin excitations in 


copper oxides decreases dramatically with increasing hole 
doping. Although there are only limited data a vailable 
for hole-doped iron pnictides (Wang et al. 2013b), they 
show a similar trend as that of hole-doped copper oxides. 
More inelastic neutron scattering experiments on hole- 
doped iron pnictides are necessary in order to make a 
detailed comparison between spin excitations in iron and 
copper based high-T c superconductors. 



As discussed in Section H, recent advances in RIXS 
have allowed a direct study of spin excitations in cop¬ 
per and iron-based superconductors. However, the re¬ 
sults obtained by RIXS for hole-doped iron pnictides dis¬ 
agree with those obtained by inelastic neutron scattering 
(Wang et al. 2013b Zhou et al. 2013aI. RIXS measure¬ 


ments on copper oxide superconductors also find excita¬ 
tions that soften no more than 10% with doping, with 
negligible change in integrated intensity. While these 
results are very similar to those obtained on iron pnic¬ 
tides (Zhou et al. 2013a), they are again in conflict with 


those of inelastic neutron scattering probing spin exci¬ 
tations near Qaf = (0.5,0.5) ( |Tranquada et aZ7[ |2014| ). 
As the RIX S cross section i nclud es both charge and spin 
excitations ( Ament et al.\ 2011), it is not known how 
to compare directly the RIXS intensity with that of the 
well-established magnetic cross section obtained from in¬ 
elastic neutron scattering. 

In the case of electron-doped copper oxide supercon¬ 
ductors, recent RIXS experiments have confirmed that 
magnetic excitations harden across the antiferromag¬ 
netism and superconductivity boundary seen originally 


with neutron scattering experiments ( 

Fujita et al. 2006 

Lee et al. 

2013b 

Wilson et al. 2006b 

). However, RIXS 


finds a negligible change in the integrated magnetic in¬ 
tensity as a function of electron doping at high energies, 


FIG. 41 (Color online) Superconductivity-induced changes 
in spin dynamic susceptibility and magnetic contribution 
to the superconducting condensation energy in CeCu 2 Si 2 . 
(a) Energy scans in S'-type CeCu 2 Si 2 at Q = Qaf = 
(0.215, 0.215,1.458) in the superconducting state at zero field 
and normal state at B = 2 T (T = 0.07 K). (b) Imagi¬ 
nary part of the dynamic susceptibility at Qaf in the normal 
(x(v(Qaf> and superconducting (Xs(Qaf, hui)) states. 
The blue area (marked with a +) leads to an increase in mag¬ 
netic exchange energy A E x , whereas the green area (marked 
with a —) leads to a decrease in A E x (Stockert et al ., 2011). 


assuming that the excitations measured by RIXS near 
the r (Q = 0) point are the same as those near Qaf 
obtained by neutron scattering for a doped metallic sam¬ 
ple (Tranquada et al. 2014). Whereas the neutron scat¬ 


tering cross section is well understood, the RIXS cross 
section is complicated and the significance of the RIXS 
measurements is unclear at present. Only future RIXS 
and neutron scattering measurements performed on the 
same sample at the same region of the reciprocal space 
will shed new light to our understanding of the spin and 
electronic excitations as revealed by RIXS and determine 
their significance to high-T c superconductivity. 

In addition to copper and iron based supercon¬ 
ductors, unconventional superconductivity also includes 


measured in neutron scattering experiments with RIXS 
in the AF ordered phase of copper oxide and iron-based 
materials due to symmetry of the equivalent Brillouin 
zones (see dashed box near T and magenta shaded box 
near Qaf in Fig. 3), there is no physical justification for 


again in conflict with neutron scattering results on simi- 

heavy fermion superconductors ( 

Gegenwart et al. 

2008 

lar materials 

(Fujita et al. 

2006 

Lee et al. 

2013b 

Wilson 


v. Lohneysen et al., 2007] Steglich et al., 1979; Stewart 

et al. 

2006b) 

While one can directly compare spi 

n waves 

2001 

). Compared with iron and copper-based supercon- 


ductors, the parent compounds of heavy fermion super¬ 
conductors are also long-range ordered antiferromagnets 
but with a magnetic exchange coupling energy scale much 
lower than that of the AF parents of copper oxide and 
iron-based superconductors. For example, in a recent 






































































































36 


study of spin waves in AF ordered CeRhlns, the par¬ 
ent compound of the CeCoIns family of heavy Fermion 
superconductors (Thompson and Fisk 2012), the dom¬ 


inant in-plane nearest neighbor magnetic exchange cou¬ 
pling is SJ 0 = 0.74 nreV (Das et al. 2014), much smaller 


than that of the AF parents of iron based superconduc¬ 
tors (see Table II). In spite of the dramatically reduced 
energy scale, heavy Fermion superconductors still have 
some interesting features also seen in copper and iron- 
based superconductors. A case in point is the neutron 
spin resonance in heavy Fermion superconductors such 


as UPd 2 Al 3 (Sato et al. 2001) and CeCoIn 5 (Stock et al. 

20081 ). 


Figure 40(a) shows energy dependence of the spin ex¬ 
citations of CeCoIns below and above T c at Qaf = 
(0.5, 0.5, 0.5) ( |Stock efal. 2008). While the normal state 
spin excitations are featureless within the probed energy 
range, superconductivity induces a sharp peak reminis¬ 
cent of the resonance in iron and copper based supercon¬ 
ductors. From neutron polarization analysis of the reso¬ 
nance discussed in Section F, we know that the mode in 
iron-based superconductors is isotropic, consistent with 
the singlet-to-triplet excitation of electron Cooper pairs 
(Eschrig, 2006). If this is indeed the case, application of 
a magnetic field should Zeeman split the resonance into 
three peaks a rising from a single ground state to a triplet 
excited state (Dai et al. 2000), where the central field in¬ 


dependent peak is longitudinally polarized while the field 


dependent peaks are transversely polarized (Ismer et al. 


2007). Unfortunately, magnetic field experiments on cop- 


per oxide (Dai et al. 

2000) and 

iron-based superconduc- 

tors 

Li et al., 2011a 

Wen et al. 

2010 

Zhao et al. 

2010) 


have been unable to determine the ground state of the 
resonance. Surprisingly, the application of a magnetic 
field on heavy Fermion superconductor CeCoIn 5 splits 
the resonance into two peaks as shown in Figs. 40(b)- 


40(d) (Stock et al. 2012). This would suggest that the 


resonance in CeCoIns is a doublet instead of a singlet-to- 
triplet excitation. In separate polarized inelastic neutron 
scattering experiments under a magnetic field, the dou¬ 
ble peak nature of the resonance was confirmed under a 2 


T applied field (Raymond et al. 2012). In addition, the 
resonance line shape is found to depend on the neutron 
polarization, suggesting that the resonance is a degener¬ 
ate mode with three fluctuation channels: a Zeeman split 


contribution and an additional longitudinal mode (Ray¬ 


mond et al. 2012|. While these results on CeCoIn 5 are 


interesting, they still have not established conclusively 
the doublet or singlet-to-triplet nature of the resonance. 

Using established models to calculate the magnetic 
contributions to the superconducting condensation en- 


ergy in copper oxide superconductors 

Dernier and 

Zhang 

1998 

Scalapino and White 

1998 

, one can in 


principle estimate the lowering of the magnetic exchange 
energy in Yf^CuaOg+z family of m aterial s using spin 
excitation spectra above and below T c (|Dahm et al. 2009 


Dai et al. 1999 Woo et al. 2006). However, the large 


ide superconductors means that it is difficult to obtain 
the overall spin excitation spectra in the low-temperature 
normal state using a magnetic field to eliminate super¬ 


conductivity (Dahm et al.. 2009 Dai et al. 1999 Woo 


et al. 2006). This makes the estimation of the lowering 


of the magnetic exchange energy in the low-temperature 
state ambiguous (Dernier and Zhang 1998 Scalapino and 


White 1998). Since heavy Fermion superconductors gen¬ 


erally have a much smaller energy scale, one can easily 
obtain the low-temperature normal state spin excitation 
spectra using a magnetic field to suppress superconduc¬ 
tivity. Figure 41(a) shows the energy dependence of the 
spin excitation spectra in the low-temperature normal 


and superconducting states for CeCi^SU (Stockert et al. 


2011), which exhibits superconductivity below T c « 0.6 
K (Steglich et al. 1979). One can see the opening of a 


spin gap and clear enhancement in spin excitations above 
the gap. Using the energy dependence of the dynamic 
susceptibility in absolute units in the low-temperature 
normal and superconducting states [Fig. 41(b)], Stock- 
ert and co-workers estimated the lowering of the mag¬ 
netic exchange energy in CeCi^SU and found it to be 
considerably larger than the superconducting condensa¬ 
tion energy (Stockert et al. 2011). Although the lowering 


of the magnetic exchange energy is also found to be much 
larger than the superconducting condensation energy in 
the superconducti ng iron pnictides (Wang et al. 2013b) 
and chalcogenides ( |Leiner et al. 2014| ), these experiments 
also suffer the problem of not being able to determine the 
overall spin excitation spectrum in the low-temperature 
normal state. 


IV. THEORETICAL DESCRIPTIONS OF STATIC AF 
ORDER AND SPIN EXCITATIONS IN IRON-BASED 
SUPERCONDUCTORS 

Shortly after the discovery of iron-based superconduc¬ 
tors, band structure calculations predicted that the Fermi 
surfaces of parent compounds consist of two quasi-two- 
dimensional near-circular hole pockets centered around 
the zone center F, and two quasi-two-dimensional elliptic 
electron pockets center around the (1,0) and (0,1) points 
in the orthorhombic unfolded Brillouin zone (Fig. 14) 
(Hirschfeld et al.. 2011 Mazin 2010). The sign reversed 


quasiparticle excitations between the hole and electron 
pockets (nesting) can give rise to the spin-density-wave 


order observed experimentally (de la Cruz et al. 


2008 


Dong et al. [2008). In addition, the same quasiparti¬ 


cle excitations in doped superconductors are expected 
to induce a neutron spin resonance in the superconduct- 


III. B) (Korshunov and Eremin 

00 

o 

o 

CM 

Maier et al. 

2009 

Maier and Scalapino 2008 

). Within this weak-coupling 


energy scale of the spin excitations in high-T c copper ox¬ 


analysis, iron-based superconductors and their parents 
are assumed to be good metals made of itinerant elec¬ 
trons with spin-density-wave type AF order. Spin waves 
and spin excitations can then be calculated using RPA in 
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a multiband Hubbard model with appropriate Fermi sur¬ 
faces for hole and electron pockets ( |Kno lle and Eremin 


2013). In this approach, the large in-plane effective mag¬ 


netic exchange coupling anisotropy (Table II) in the spin 


waves of iron pnictides (Harriger et al. 2011 Zhao et al. 


2009) can be understood as due to the ellipticity of the 


electron pockets (Fig. 14), which induces frustration be¬ 
tween the (1, 0) and (0, 1) wave vectors connecting the 


hole and electron pockets (Kaneshita and Tohyama 2010 


Knolle et al.. 

2010|. In the 

underdoped regime where the 
d competes with superconduc- 

static AF order coexists an 

tivity 

(Christianson et al. 

2009 

Pratt et al. 

2009 

Wang 

et al. 

2010 

), spin excitations at (1,0) are determined by 


the presence of the AF order and associated spin waves, 
while the excitations at wave vector (0,1) are dominated 
by the superconductivity and formation of the resonance 
in the single domain sample (Knolle et al. |2011). How¬ 


ever, since most neutron scattering experiments in under¬ 
doped iron pnictides were carried out on twinned sam¬ 
ples that cannot distinguish the wave vector (1, 0) from 
(0,1) (Christianson et al. 2009 Pratt et al.. 2009[ Wang 


et al. 


2010), it is unclear how the resonance associated 


with superconductivity interacts with itinerant electrons 
contributing to the spin waves. In a systematic study of 
spin excitations in BaFe2-r C Ni a; As2, the electron-doping 
evolution of the low-energy spin excitations was found to 
qualitatively agree with RPA calculations of the nested 
Fermi surfaces (Luo et al. 2012b). However, the high- 


energy spin excitations are weakly electron-doping inde¬ 
pendent, and have values much different from that found 


by RPA calculations (Liu et al. 2012a Luo et al. 2013a 


Wang et al. 2013b). These results suggest that the weak- 
coupling analysis based on purely itinerant electrons is 
insufficient to explain the entire spin excitation spectrum 
and its electron or hole doping evolution. 

Although the weak-coupling approach using a Fermi 
surface nesting picture provides a nice framework to un¬ 
derstand static AF order, spin excitations, and their con¬ 
nection with superconductivity in iron-based materials 
( Hirschfeld et al.\ 2011 Mazin |2010| ), calculations show 
that the ordered moment in the AF iron pnictides is 
around 2 /rs/Fe (Ma and Lu, |2008 ), much larger than 
the largest experimental value (~0.9 hb/ Fe, see Table I). 
Furthermore, the ordering wave vectors of the bi-collinear 
AF structure in iron chalcogenides shown in Fig. 2(c) 
do not match the nesting wave vector of the Fermi sur¬ 
faces (Subedi et al. 2008). In the strong coupling limit, 


all unpaired electrons, not just itinerant electrons near 
the Fermi surface, participate in forming magnetic order, 
much like the magnetic moment of Cu 2+ in the insulat- 


large in-plane magnetic exchange coupling anisotropy in 
the parent compounds of iron pnictides is understood in 
terms the presence of the biquadratic exchange coupling 
K between the nearest spins in the AF ordered states, 
which can be mapped onto the Ji a -Jib model with a 
specif i c relationship between J\ a -Jib and J\-K (Wysocki 
et al . | 201 1[ Yu et al. 2012). This means that the usual 
local moment Heisenberg Hamiltonian will be modified to 
H = Ji Ei,i Si-S t , s + J 2 Ei,s S,-S t+J —A^^S-rS.j) 2 , 
where J\ and J2 are the nearest and next nearest neigh¬ 
bor exchan ge c ou plings, respect ivel y, an d S, is the spin 
at site i (Wysocki et al.. 2011 Yu et al. 2012| ). The 
calculated dynamical structure factor S(Q,uj) forms el¬ 
lipses in momentum space around Qaf = (±1,0) at low 
excitation energies. With increasing energy, the elliptic 
features expand towards the zone boundary and gradu¬ 
ally split into two parts, forming a pattern around the 


experiments (Harriger et al. 

2011 

Wysocki et al. 

2011 

Yu et al. 

2012 Zhao et al. 

2009) 



By comparing the AF exchange interactions deter¬ 
mined from neutron scattering in the parent compounds 
of iron-based superconductors with those of the super¬ 
conducting gap functions in the corresponding supercon¬ 
ducting materials observed by ARPES, it was argued 
that the key ingredients in determining the high-T c su¬ 
perconductivity of iron-based materials are the next near¬ 
est neighbor local AF exchange interactions in real space 
and a Fermi surface topology in reciprocal space that 
matches to the pairing form factor provided by AF in¬ 
teractions (Hu and Ding 2012). From analysis of spin 


excitation spectra in hole and electron-doped iron pnic¬ 
tides and magnetic contributions to the superconduct¬ 
ing condensation energy (Scalapino 2012), we find that 


high-T c superconductivity is associated with materials 
having large magnetic exchange coupling and strong itin¬ 


erant electron-spin excitation interactions (Wang et al. 


2013b), similar to the large Debye temperatures and 


strong electron-phonon interactions in high-T c BCS su¬ 
perconductors (Bardeen et al. 1957). Finally, the strong¬ 


coupling approach predicts the presence of a quantum 
critical point separating a paramagnetic metallic phase 
from an AF phase phase in phosphorus doped iron pnic¬ 


tides ([D ai et al.[ |2009). Indeed, neutron scattering exper¬ 
iments on powder samples of CeFeAsi-^P^O suggest the 
presence of a quantum critical point near x = 0.4 con¬ 
trolled by the pnictogen height away from the Fe-plane 
[Fig. 7(a)] ( |de la Cruz et al. 2010). Similarly, thermody¬ 
namic and transport measurements on BaFe2(Asi_ x P a; )2 
have identified a quantum critical point near optimal su- 


ing copper oxides (Fang et al. 2008 Si and Abrahams perconductivity at x = 0.3 [Fig. 6(d)] (Shibauchi et al. 


2008 

Xu et al. 

2008 

). Here, the AF ordered state of 

2014 


iron-based superconductors can be described by a local 
moment Heisenberg Hamiltonian with the nearest neigh¬ 
bors (Jia-Jib) and next nearest neighbor (J2) exchange 


interactions (Fig. 4) (Dai et al. 

2009 Fang et al. 

2008 

Han et al.\ 2009 [ |Moreo et al. 

2009( Xu et al. 

2008 

Yildirim 

2008] [Yu et al. 

2013). In this picture, the 


surements on BaFe 2 (Asi_ 2 ,P a; )2 suggest that structural 
quantum criticality cannot exist at compositions higher 
than x = 0.2 8 due to the actual phase stability range 
(Allred et al.\ 2014). It is therefore desirable to carry 


out neutron scattering and SR experiments on single 
crystals of BaFe2(Asi_ a: P a ;)2 near x = 0.3 to establish 
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the nature of the AF phase transition with increasing 
P-doping. 

While there are many phenomenological reasons for us¬ 
ing the strong coupling approach to understand the elec¬ 
trical transport, spin and charge dynamical properties of 
iron pnictides and chalcogenides (|Dai et al. 2009 Fang 
eTalj [20081 |Xu~ et al. ] |2U08l |Yildirim[ |2008| |Yu et al.\ 


2013), such an approach is incompatible with some spin 


dynamical properties. For example, in a strict local mo¬ 
ment Heisenberg Hamiltonian, single particle excitations 
or spin waves should only have transverse components 
and would not support longitudinal spin excitations in 
the AF ordered phase of iron pnictides as seen in polar- 
ized neutron scattering experiments [Fig. 29(f)] (|Wang 


et al. 2013a). Furthermore, the electron and hole-doping 


evolution of the low-energy spin excitations are consis¬ 
tent with the Fermi surface nesting predictions, but it is 
unclear whether the data is also compatible with a pure 
local moment Heisenberg Hamiltonian (see Section III. 
C). Finally, spin waves of iron chalcogenides are heav¬ 
ily damped at high-energies near the zone boundary and 
exhibit a number of anomalous properties difficult to un¬ 


derstand within a local moment model (Lipscombe et al. 
2011 Zaliznyak et al. 2011). 


Instead of a strong or weak coupling approach, the iron 
pnictides may be Hund’s metals where the interaction be¬ 
tween the electrons is not strong enough to fully localize 
them to form a Mott insulator, but is sufficient so that 
the low energy quasiparticles have much enhanced mass 
(Haule and Kotliar 2009). Here the electron correlation 


strength would be primarily controlled by the Hund’s 
coupling Jh , which depends on the pnictogen heights 
and tends to align spins of all the electrons on a given 
Fe-atom, and hence enhances spin excitations without 


appreciably affecting the charge excitations (Yin et al. 


2011). This is different from the effect of large Coulomb 
repu sion U in a Mott insulator, which hampers charge 
excitations in order to enhance spin fluctuations (Lee 


et al. 2006). The electronic excitations in iron-based 


superconductors are neither fully itinerant nor fully lo¬ 
calized, but have a dual nature that can be realistically 


described by a combination of DFT and DMFT (Kotliar 


et al. 2006). This idea is similar to the picture where 


single electron spectral function is composed of coherent 
and incoherent parts representing electrons near (itiner¬ 
ant electrons) and far away (local moments) from the 
Fermi surface (Abrahams and Si 2011). 


Using the combined DFT and DMFT method, one 
can estimate the size of the ordered moment for differ¬ 
ent iron pnictides and find them to be close to the ob¬ 
served value (Yin et al. |2011 ). The same method has also 
been used to calculate the spin wave spectra in BaFe2As2 
and good agr eement was found with neutron scat tering 
experiments ( |Park et al. 2011a Yin et al. 20141. Fi¬ 
nally, the combined DFT and DMFT method has been 
used to calculate the electron and hole-doping depen¬ 
dence of the spin excitation spectrum in absolute units 


for Bai_ a ,K a; Fe2As2 and BaFe2- a; Ni a: As2 (Wang et al. 


2013b). Compared with RPA calculations based on the 


weak-coupling approach (Hirschfeld et al. 2011), com¬ 


bining DFT and DMFT gives a more realistic estimation 
of the absolute intensity of the local dynamic susceptibil¬ 
ity, and can quantitatively model the electron and hole¬ 
doping evolution of the spin excitations in absolute units 
(Wang et al. 2013b). Furthermore, it can account for 


the pni ctoge n hei ght de pendence of the spin-wave spec¬ 
tra (Zhang et al.. 2014aI. 


The static AF order and spin excitations in iron-based 
materials can also be understood by hybrid models con¬ 
sisting of local moments on each Fe site and itinerant 
electrons from the degenerate d xz and d yz orbitals (Kou 


et al. 2009 Yin et al. 2010). In this picture, the io- 


cal moments interact with each other via J\ and J2 
Heisenberg exchanges, and they are coupled to the itin¬ 
erant electrons via Hunds rule coupling. Since itinerant 
electrons are only associated with d xz and d yz orbitals 
that break the C4 rotational symmetry of the under¬ 
lying x-y lattice plane due to their different occupan¬ 
cies, these orbitals can form a Hamiltonian that drive 
the in-plane magnetic anisotropy, producing unfrustrated 
collinear AF order and lifting the degeneracy of the (1,0) 


and (0,1) magnetic states (Chen et al.. 2010a 


2009 

Kruger et al., 2 

2010 

Yin et al. 

2010) 


Lee et al., 2009 


Kou et al. 


Lv et al. 


is due to purely electronic ferro-orbital order that spon¬ 
taneously breaks the rotational symmetry of the under¬ 
lying lattice and drives the observed magnetic and struc¬ 
tural transitions without Fermi surface nesting or mag¬ 


netic frustration (Kruger et al. 2009 Lee et al. 2009). 


Using a fermionic representation of the local moments 
and a generalized RPA framework, one can calculate the 
spin wave spectra of BaFe2As2 and find that the out¬ 
come is consistent with spin excitations in the AF ordered 
and paramagnetic states obtained from inelastic neutron 
scattering (Leong et al. 2014 Yang et al. 2010). In 


addition, the global phase diagram for the AF and su¬ 
perconducting states calculated from the hybrid model 
on the mean-field level is qualitatively consistent with 
experiments (You et al. 2011). At high characteristic 


temperatures, electrons in more localized orbitals of the 
multiband system may first form short-ranged AF or¬ 
der. Upon cooling to lower temperatures, the electrons 
in more itinerant orbitals can be driven into a true static 
AF ordered or superconducting state via Hund’s coupling 
to the preformed localized AF state. This is analogous to 
the orbital-selective Mott transition, where itinerant and 
localized electrons in different orbitals may separate as 
independent degrees of freedom (Kou et al. 2009 Leong 


et aZ-l |2014| |Lv et al.\ |2010[ [Yang et aL[|20l6||You et al. \ 

20TT|. 


V. SUMMARY AND OUTLOOK 

In this paper, we have reviewed recent progress of 
neutron scattering studies of the static AF order and 
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spin dynamics in iron-based high temperature supercon¬ 
ductors. Soon after the discovery of these materials in 
2008 (Kamihara et al.\ 2008), neutron diffraction mea¬ 


surements at the NIST center for neutron research and 
high flux isotope reactor at Oak Ridge National Lab¬ 
oratory have determined the AF order and crystalline 
structures of the parent and superconducti ng compounds 
(de la Cruz et al. 2008 |Huang et al. 2008a I. These mea¬ 


surements have established the basis that superconduc¬ 
tivity in iron-based materials arises from the suppression 
of static long-range ordered antiferromagnets, much like 
copper oxide superconductors (Tranquada et al. 20141. 


When single crystals of the 122 family of iron pnictides 
and iron chalcogenides have become available, the ad¬ 
vanced time-of-flight neutron spectrometers at spallation 
neutron sources at Oak Ridge National Laboratory and 
ISIS at Rutherford-Appleton Laboratory have allowed 
detailed mapping of th e sp in wav e spectra t hroughout the 
Brillouin zone (Diallo et al. 2009 Zhao et al. 2009). This 


first occurred only slightly more than one year after the 
discovery of the iron pnicti de superco nd uctors . In cop¬ 
per oxide superconductors (Bednorz and Muller, 1986), 
the first complete spin wave spectrum was measured 15 


years after its discovery (Coidea et al. 2001). Using the 


overall spin wave spectra in the AF ordered iron pnic¬ 
tides, one can fit the dispersion curves with a Heisenberg 
Hamiltonian, revealing the anisotropic in-plane effective 
magnetic exchange couplings. This has inspired much 
discussion on the microscopic origin of the effective mag¬ 
netic anisotropy as described in Section IV. 

Since high quality single crystals of electron and hole- 
doped BaFe2As2 are available, most of the elastic and 
inelastic neutron scattering experiments have been car¬ 
ried out on these materials. With elastic neutron scatter¬ 
ing experiments, one can map out the electron and hole 
doping evolution of the structural and magnetic phase 
diagrams. For electron-doped materials obtained via Co 
and Ni substitution of Fe, the tetragonal-to-orthorhombic 
lattice distortion precedes the AF phase transition, and 
the static long-range order coexists/competes with super¬ 
conductivity in the underdoped regime. However, the AF 
order becomes incommensurate with a short-range corre¬ 
lation length near optimal superconductivity, indicating 
that it is a spin-glass phase in the matrix of the supercon¬ 
ducting phase, coexisting and competing with supercon- 


ductivity 

(Bernhard et al. 

2012 

Dioguardi et al. 2013 

Lu et al. 

2014b). NMR measurements on 1111 family of 


materials also suggest nanascale electronic inhomogene¬ 
ity (Lang et al. 2010). For hole-doped Ba 1 _ a ,K a ,Fe2As 2 


and Bai_a,Na x Fe2As2, the structural and magnetic phase 


transitions are coupled at at all doping levels (Avci et al. 


2013 2014 2012). Near optimal superconductivity, there 


is a new magnetic phase within the tetragonal structure, 
possibly associated with a spin nematic phase, and the 
transition from AF order to superconductivity may also 
occur in the first order fashion (Avci et al. 2013 2014). 


Finally, although transport and NMR measurements sug¬ 
gest the presence of a quantum critical point near x = 0.3 


for isoelectr onically doped BaFe2(Asi_ ;r P ;E )2 (,Shibauchi 


et al. 2014), neutron diffraction measurements have only 


mapped out the magnetic and structural pha se diagram 
in the underdoped regime (Allred et al. 20141, and much 


is not known about the evolution of AF order near opti¬ 
mal superconductivity. 

Similar to neutron diffraction work, most of the in¬ 
elastic neutron scattering studies of spin excitations 
in iron pnictides have been focused on electron-doped 
BaFe2As2- Compared with the undoped parent com¬ 
pounds, electron-doping appears to modify spin excita¬ 
tions below ~80 meV while leaving high energy spin ex¬ 
citations mostly unchanged (Luo et al. 2013a). How¬ 


ever, hole-doping suppresses high-energy spin excitations 
and transfers the spectral weight to low-energies (Wang 


et al. 2013b). In addition, the wave vector dependence 


of the low-energy spin excitations in iron pnictides ap¬ 
pears to be controlled by the quasiparticle nesting be¬ 
tween the hole and electron Fermi surfaces. These re¬ 
sults are consistent with the notion that spin excitations 
in iron pnictides have both local and itinerant character 
with the electron correlations controlled by the pnictogen 


Kotliar 

2009 Kou et al , 

2009; |Yin et al ., 

2010) 


The 

availability of large single crystals of iron chalcogenides 


terials have been carefully mapped out ( 

Lipscombe et al. 

2011 Lumsden et al. 

2010 

Zaliznyak et al. 

2011 

). In 


particular, application of the sum rules of neutron scat¬ 
tering indicate that the integrated spin excitation in¬ 
tensity of Fei+3,Te is inconsistent with an S = 1 Fe 2+ 
ground state expected in t he pre s ence o f a strong crys¬ 
talline electric field (Stock et al. 2014), suggesting the 


importance of itinerant electrons even for the iron chalco¬ 
genides, which exhibit strong electron correlations and 
localized moments (Yin et al. 2011). 


Compared with electron-doped BaFe2As2, spin excita¬ 
tions in hole-doped Bai_ a; K a ;Fe2As2 and isoelectronically 
doped BaFe 2 As2-xP :E iron pnictides have been much less 


studied ( 

Lee et al. 

Zhang et al. 

2011a 


2011 2013a Wang et al. 2013b 


possible spin nematic phase in the tetragonal phase (Avci 


et al. 2014), it will be interesting to study the evolution 


of the overall spin excitations in hole and isoelectroni¬ 
cally doped iron pnictides. In particular, since the elec¬ 
tron pairing symmetry of the heavily hole-doped super¬ 


conducting Bai_ a ,K x Fe2As2 is still unclear (Ota et al. 


2014 Tafti et al. 2013), it is important to carry out 


temperature dependent measurements to study the ef¬ 
fect of superconductivity on low-energy spin excitations. 
A determination of the wave vector and energy of the 
superconductivity-induced neutron spin resonance will 
put considerable constraint on the nature of the super¬ 
conducting pairing state. 

Although most neutron scattering work has focused on 
the 122 family of iron pnictides and iron chalcogenides, 
the 111 family, including Co-doped NaFeAs and LiFeAs, 
is equally interesting since these materials may be more 



































































































































































40 


correlated than the 122 family (Yin et al.\ |201 1). At 
present, spin waves throughout the Brillouin zone have 
been mapped out for NaFeAs (Zhang et al. 2014aI. It is 
important to determine how electron-doping affects the 
spin excitations and to compare the outcome with the 
pure LiFeAs and Co-doped LiFeAs. Similarly, it is impor¬ 
tant to study temperature and doping dependent spin ex¬ 
citations in Se-overdoped Fei+j,Tei_ x Sea; and pure FeSe. 
The case of pure FeSe is particularly interesting as this 
is the system where the structural phase transition hap¬ 


pens without static AF order (Johnston 2010 Stewart 


2011). A complete understanding of this material may 


reveal spin or orbital driven electronic nematic phase. 
Detailed experiments on other iron-based superconduc¬ 
tors and associated materials are necessary to establish 
the common features of the magnetism in various mate¬ 
rials and their connection to high-T c superconductivity. 
Neutron scattering, together with RIXS, /rSR, and NMR, 
can play a unique role in our quest to find the microscopic 
origin of high-T c superconductivity. 
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